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Abstract 

The  thermal  neutron  capture  cross  section  of  the  rare  earth  metal  isotope  157Gd 
is  255,000  barns  and  is  the  largest  of  all  known  natural  isotopes,  which  distinguishes 
the  material  as  a  logical  candidate  for  neutron  detection.  Semiconductors  that  in¬ 
corporate  large  neutron  absorption  cross  section  elements,  such  as  gadolinium,  are 
gaining  attention.  Despite  the  relative  immaturity  of  the  development  of  such  semi¬ 
conductors,  gallium  nitride  is  one  such  semiconductor  material  that  has  attained 
considerable  recognition.  In  spite  of  study  over  the  past  two  decades,  there  still 
exists  an  incomplete  understanding  of  the  physical  and  electronic  properties  of  GaN 
materials,  particularly  whether  even  low  concentrations  of  a  rare  earth  in  a  GaN  host 
can  alter  the  surface  electronic  structure.  To  address  this  concern,  investigations  of 
the  surface  electronic  structure  and  interface  properties  of  GaN  thin  films  doped 
with  three  distinct  rare  earths  (Yb,  Er,  Gd)  were  undertaken  using  photoemission 
spectroscopy. 

The  effective  Debye  temperatures  of  ytterbium  and  gallium  in  GaN:Yb  thin  films 
were  obtained  using  X-ray  photoemission  spectroscopy.  The  vibrational  motion  nor¬ 
mal  to  the  surface  resulted  in  a  dimunition  of  photoemission  intensities  from  which 
the  effective  Debye  temperatures  of  221  ±  30  K  and  308  ±  30  K  for  Yb  and  Ga, 
respectively,  were  estimated.  The  similarity  between  the  Yb  and  Ga  Debye  tempera¬ 
tures  is  indicative  of  a  substitutional  occupation  of  a  Ga  site  by  a  Yb  ion.  The  slightly 
smaller  effective  surface  Debye  temperature  for  Yb  correlates  to  a  soft,  strained  sur¬ 
face,  possibly  due  to  an  increased  Yb-N  bond  length  as  compared  to  the  Ga-N  bond 
length. 

Bonding  with  the  GaN  and  rare  earth  4/  hybridization  were  also  examined.  The 
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4 d  — >  4/  super  Coster- Kronig  resonances  for  various  rare  earth  doped  GaN  thin  films 
(RE  =  Yb,  Er,  Gd)  were  investigated  using  synchrotron  photoemission  spectroscopy. 
Resonant  photoemission  Fano  profiles  showed  that  the  major  Gd  and  Er  rare  earth 
4/  weight  was  at  about  5-6  eV  below  the  valence  band  maximum,  similar  to  the  4/ 
weights  in  the  valence  band  of  many  other  rare  earth  doped  semiconductors.  For 
Yb,  there  was  very  little  resonant  enhancement  of  the  valence  band  of  Yb  doped 
GaN,  consistent  with  a  largely  4/14  occupancy.  The  placement  of  the  occupied 
Gd,  Er  and  Yb  4/  levels,  deep  within  the  valence  band,  suggested  that  the  intra- 
atomic  /-/  transitions  may  be  more  ‘blue’  than  predicted  by  theoretical  models.  The 
Schottky  barriers  formed  at  the  interface  between  gold  and  the  various  rare  earth 
doped  GaN  thin  films  were  also  investigated  in  situ  using  synchrotron  photoemission 
spectroscopy.  We  found  compelling  evidence  that  thin  layers  of  gold  do  not  wet  and 
uniformly  cover  the  GaN  surface,  even  with  rare  earth  doping  of  the  GaN.  The 
resultant  Schottky  barrier  heights  were  measured  as  1.68  ±0.1  eV  (GaN:Yb),  1.64 
±  0.1  eV  (GaN:Er),  and  1.33  ±  0.1  eV  (GaN:Gd),  which  were  25-55%  larger  than 
those  reported  at  the  gold  to  undoped  GaN  interface.  The  trend  of  the  Schottky 
barrier  heights  followed  the  trend  of  the  rare  earth  metal  work  function. 

The  utility  of  gadolinium  as  a  neutron  detection  material  in  a  hypothetical  direct 
conversion  homojunction  was  examined  via  fundamental  nuclear  and  semiconductor 
physics.  In  the  single  neutron  detection  limit,  the  charge  produced  in  a  Gd-based 
semiconductor  by  an  internal  conversion  electron  created  from  15 'Gd  neutron  capture 
was  found  to  approach  the  intrinsic  noise  level  of  commercial  preamplifiers  required 
for  pulse  height  analysis.  The  fraction  of  captured  neutrons,  and  hence  efficiency, 
was  maximized  with  Gd-layer  thickness,  and  the  range  of  internal  conversion  elec¬ 
trons  created  from  neutron  capture  by  15 'Gd  was  estimated  to  be  ~  20-25  /im.  Full 
depletion  of  a  Gd-based  homojunction  to  the  25  /im  range  of  the  internal  conversion 
electrons  was  seen  to  be  contrary  to  low  power  operation. 
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THE  EFFECTS  OF  RARE  EARTH  DOPING  ON  GALLIUM  NITRIDE  THIN 


FILMS 


I.  Introduction 


1.1  Motivation 

Nuclear  radiation  detection  has  long  been  one  of  the  most  troublesome  homeland 
security  challenges,  and  is  not  exclusively  a  post-9/11  concern.  Months  after  the 
United  States  used  the  Erst  nuclear  weapons  against  Japan,  apprehensive  lawmakers 
asked  Robert  Oppenheimer,  widely  regarded  as  “The  Father  of  the  Atomic  Bomb,” 
how  authorities  would  detect  a  nuclear  weapon  in  an  incoming  shipping  crate.  He 
replied,  “with  a  screwdriver,”  meaning  the  only  device  that  would  do  the  job  was 
a  screwdriver  to  open  every  incoming  box  and  crate  [1].  Many  of  the  conclusions 
from  the  aptly  dubbed  and  still  classified  “Screwdriver  Report,”  co-authored  by 
top  physicists  and  Manhattan  Project  veterans  Robert  Hofstadter  and  Wolfgang 
Panofsky,  remain  largely  unchanged  today.  Consequently,  in  2007,  while  serving  as 
the  director  emeritus  of  the  Stanford  Linear  Accelerator  Center,  Panofsky  defended 
the  validity  of  the  report’s  assertions,  stating  that  “the  laws  of  physics  have  not 
changed  one  bit,”  and  that  nuclear  devices  remain  virtually  undetectable  “unless 
you  are,  say,  ten  feet  away  from  it-and  even  then  it  can  be  quite  easily  shielded”  [2], 

While  the  central  point  of  Panofsky’s  report  is  still  valid,  in  that  radioactive  ma¬ 
terial  emits  less  radiation  with  distance  and  shielding  from  the  source,  not  everyone 
completely  agrees  with  him;  Department  of  Homeland  Security  (DHS)  officials  claim 
that  state-of-the-art  detectors  offer  far  better  performance.  Nevertheless,  DHS  con¬ 
tinues  its  vigilant  pursuit  of  improved  radiation  detector  design,  performance,  and 
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efficiency. 

The  creation  of  DHS  provided  the  impetus  to  invigorate  a  radiation  detector 
market  considered  by  many  to  be  saturated  or,  at  best,  growing  slowly  because 
neither  incentives  to  invest  in  new  technology  nor  demands  for  new  applications 
existed.  In  2002,  DHS  enlisted  the  assistance  of  the  American  National  Standards 
Institute  (ANSI),  which  oversees  the  development  of  voluntary  consensus  standards 
for  products,  services,  processes,  systems,  and  personnel  in  the  United  States.  DHS 
and  ANSI  collaborated  to  develop  and  adopt,  within  an  unprecedented  two-year 
time  frame,  four  American  national  standards  designed  to  assist  federal  agencies,  as 
well  as  state  and  local  officials,  in  procurement  decisions  related  to  radiological  and 
nuclear  detection  equipment  [3]. 

These  new  standards,  coupled  with  DHS  funding,  stimulated  demand  for  innova¬ 
tive  ideas,  technology,  and  instruments,  thus  improving  miniaturization,  efficiency, 
accuracy,  and  operation  of  radiation  detectors  for  homeland  security  users.  As  one 
example,  Instrument  Business  Outlook  projected  14.5%  and  13.7%  global  market 
growth  rates  for  personal  radiation  detectors,  radioactive  isotope  identification  de¬ 
vices,  and  radiation  portal  monitors  for  the  years  2006  and  2007,  respectively,  with 
anticipated  increases  in  subsequent  years  [4].  Furthermore,  DHS’s  Radiation  Portal 
Monitor  Program  (RPMP)  seeks  to  deploy  a  suite  of  stationary,  mobile,  and  hand¬ 
held  radiation  detection  systems  at  all  U.S.  Customs  and  Border  Protection  (CBP) 
ports  of  entry.  The  objective  is  to  scan  100%  of  all  traffic  and  cargo  for  nuclear  and 
radioactive  materials,  while  maintaining  the  efficient  flow  of  commerce.  In  November 
2009,  DHS  Secretary  Janet  Napolitano  announced  the  RPMP’s  full  deployment  of 
non-intrusive  scanning  equipment  to  detect  radiation  emanating  from  nuclear  device 
material  across  all  northern  border  land  points  of  entry  into  the  United  States  [5]. 

While  these  efforts  are  noteworthy,  we  cannot  underestimate  our  enemies’  ambi¬ 
tions  to  circumvent  the  capabilities  of  our  current  detection  systems.  The  RPMP 
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directs  the  preponderance  of  efforts  toward  the  proactive  detection  of  radioactive 
material  at  high  volume  entry  points;  hence,  its  assets  are  mostly  fixed,  or  vehicle- 
mobile,  and  these  assets  are  concentrated  at  roadway  entry  points.  Accordingly,  we 
must  continue  to  develop  novel  radiation  detection  technologies,  which  accurately 
and  efficiently  identify  nuclear  material  via  proactive  cargo  screening.  We  must  also 
provide  our  first  responders  (firefighters,  law  enforcement,  emergency  medical)  with 
reaction  capability  during  crisis  situations,  such  as  a  hostile  dispersion  of  radiological 
material  on  American  soil. 

In  the  succeeding  sections,  analysis  of  the  radioactive  decay  modes  of  nuclear 
weapon  material  isotopes  will  show  that  neutron  detection  potentially  offers  the 
most  accurate  method  of  detecting  Special  Nuclear  Material  (SNM).  Radiation  de¬ 
tector  design  considerations  will  be  discussed  briefly,  and  the  ever-present  problem  of 
detecting  nuclear  weapon  material  in  an  incoming  shipping  crate  will  be  addressed. 

1.2  Special  Nuclear  Material  Detection 

1.2.1  Detection  Challenges. 

In  response  to  recent  worldwide  events,  our  nation’s  ability  to  monitor  and  inter¬ 
cept  special  nuclear  material  has  become  a  top  national  security  priority  [6].  Cargo 
screening  at  ports  of  entry  requires  standoff  detection  at  distances  of  several  meters 
or  more,  but  the  low  intrinsic  activity  of  special  nuclear  material,  coupled  with  the 
limited  range  of  SNM-emitted  charged  particles  and  gamma  rays,  combine  to  present 
challenges  to  standoff  detection  capability. 

Neutron  background  also  complicates  the  ability  to  discern  a  legitimate  signal 
from  persistent  noise.  Specifically,  the  cosmic  ray-induced  spallation  neutron  back¬ 
ground  of  ~  120  —  150nm~2s_1  [7,8],  which  is  due  to  cosmic  ray  interactions  with 
the  air  and  ground,  is  competitive  with,  and  sometimes  dominates,  the  signal  arising 
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from  special  nuclear  material.  The  creation  of  neutron  background  depends  upon 
atomic  number,  and  this  further  complicates  detection  of  a  SNM  source,  as  the  back¬ 
ground  due  to  cosmic  ray  interactions  at  an  air/iron  interface  (shipping  container) 
is  ~  770nm_2s_1  [8].  When  compared  to  an  unmoderated,  unshielded  hypothetical 
mass  of  5  kg  of  94%  enriched  239Pu  (weapons  grade)  at  a  distance  of  3  m,  which 
would  have  an  isotropic  neutron  flux  of  ~  3700nm~2s_1,  the  neutron  background 
clearly  affects  the  detection  process  for  nuclear  material  at  US  ports  of  entry. 

Most  radiation  detectors  operate  using  the  same  basic  principles;  radiation  enters 
a  detection  region  and  interacts  with  the  detector  material,  which  ultimately  liberates 
electrons  that  are  collected  and  measured  as  a  current  or  light  output  (scintillation). 
Therefore,  the  challenge  of  detecting  radiation  from  special  nuclear  material  requires 
an  understanding  of  the  nuclear  decay  modes  of  the  fissile  isotopes  used  for  nuclear 
weapons,  as  well  as  how  radiation,  produced  via  these  decay  modes,  interacts  with 
matter. 

Uranium-235  (235U)  and  plutoniunr-239  (239Pu)  are  the  commonly  used  isotopes 
for  nuclear  weapons.  Uranium  is  a  metal  that  is  obtained  by  mining  from  natu¬ 
ral  deposits,  and  contains  approximately  99.3%  Uranium-238  (238U)  and  0.7%  235U. 
Following  mining,  the  uranium  is  enriched  to  approximately  80%  235U,  through  a  se¬ 
ries  of  chemical  and  mechanical  processes,  at  which  point  the  material  is  considered 
suitable  for  weapon  production. 

Unlike  uranium,  only  trace  amounts  of  plutonium  exist  in  nature.  Therefore, 
239Pu  is  synthesized  in  a  nuclear  reactor  using  uranium  and  neutrons  (n),  via  beta 
decay  (/%)  with  neptunium  (Np)  as  an  intermediate,  in  the  reaction 
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which  inevitably  produces  240Pu  and  241  Pu  via  the  reactions 
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Table  1  displays  several  properties  of  the  SNM  isotopes  from  reactions  (1),  (2), 
and  (3),  and  highlights  several  concerns  regarding  their  timely  and  accurate  detec¬ 
tion.  First,  a  long  half  life,  or  low  intrinsic  radioactivity,  requires  extended  collection 


Table  1.  Selected  properties  of  common  special  nuclear  material. 


Isotope 

Half  Life 

[years] 

Primary 
Decay  Mode 

Spontaneous  Fission 
Neutron  Emission  Rate4 
[n/kg-sec] 

235  pj 

7.0  x  108 

a 

1.04  x  10"2 

238  pc 

4.5  x  109 

a 

12.6 

239  pu 

2.4  x  104 

a 

19.9 

24°  pu 

6.5  x  103 

a 

1.38  x  106 

241  Pu 

14.4 

a 

2.20  x  102 

f  Based  upon  average  neutrons  per  fission  of  2.44  and  2.89  for  uranium 
and  plutonium,  respectively  [9]. 


times  to  successfully  observe  the  infrequent  decay  events  for  the  isotopes  of  concern. 
Second,  the  dominant  decay  mode  for  these  isotopes  is  alpha  (a)  decay,  described 
by  the  emission  of  a  doubly-ionized  helium  atom  (|He2+) ,  or  alpha  particle,  as 

zX  — >■  it Y  +  *He2+  +  Q,  (4) 

where  A  represents  the  atomic  mass  or  the  sum  of  protons  and  neutrons,  Z  is  the 
atomic  number  or  the  number  of  protons,  X  and  Y  are  the  reactant  and  product 
isotopes,  respectively,  and  Q  is  the  reaction  energy.  In  theory,  Q  is  calculated  from 
the  ground  state  mass-energies  of  the  reactant,  product,  and  alpha  particle  and  by 
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conservation  of  energy  and  momentum.  The  reaction  energy  is  distributed  between 
the  product  isotope  and  the  alpha  particle,  with  the  alpha  particle  receiving  nearly  all 
of  the  reaction  energy.  In  practice,  however,  a  distribution  of  alpha  particle  energies 
exists,  due  to  the  fact  that  the  reactant  isotope  or  the  product  isotope  may  be  in 
an  excited  or  meta-stable  state  prior  to  or  immediately  following  the  nuclear  decay. 
These  distributions  are  well-known  and  tabulated  for  easy  reference  in  the  Chart  of 
the  Nuclides  [10]  or  online  via  the  National  Nuclear  Data  Center  (NNDC)  [11],  as 
two  examples.  For  the  SNM  of  greatest  concern,  the  most  probable  alpha  particle 
energies  for  235 U  and  239Pu  are  4.39  MeV  and  5.15  MeV,  respectively  [11]. 

It  should  be  noted  that  the  term  “special  nuclear  material”  in  these  discussions 
implies  that  the  mass  of  material  is  composed  of  isotopes  from  Table  1  and  contains  at 
least  a  small  percentage  of  the  isotope  240Pu.  This  assumption  is  required  because,  of 
the  isotopes  listed  in  Table  1,  only  240Pu  has  a  spontaneous  fission  neutron  emission 
rate  that  exceeds  the  background  neutron  emission  from  cosmic  ray  interactions.  In 
fact,  although  not  explicitly  stated  earlier,  the  unmoderated,  unshielded  hypothetical 
mass  of  5  kg  of  weapons  grade  material  used  for  comparison  with  neutron  background 
contained  5%  or  250  g  of  the  isotope  240Pu. 

In  addition  to  alpha  particles,  meta-stable  state  decay  and  spontaneous  fission  re¬ 
actions  yield  gamma  rays,  electrons,  and  neutrons.  The  material  response  to  charged 
particles  and  photons  classify  both  as  ionizing  radiation,  due  to  their  interaction  with 
atomic  electrons.  On  the  other  hand,  neutrons  are  classified  as  non-ionizing  radi¬ 
ation,  due  to  the  fact  that  their  interaction  is  with  the  nucleus  of  the  absorbing 
material.  Thus,  an  understanding  of  both  ionizing  radiation  and  non-ionizing  ra¬ 
diation  interactions  with  matter  is  required  in  order  to  determine  the  feasibility  to 
identify  special  nuclear  material  via  the  radiation  interaction  mechanisms. 
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1.2.2  Ionizing  Radiation  Interactions  with  Matter. 


Heavy  charged  particles,  such  as  the  alpha  particle,  interact  with  matter  primarily 
through  Coulomb  forces  between  their  positive  charge  and  the  negative  charge  of  the 
orbital  electrons  within  the  absorber  atoms.  In  any  single  encounter,  the  electron 
experiences  an  impulse  from  the  attractive  Coulomb  force  as  the  particle  passes  in  its 
vicinity.  Depending  on  the  proximity  of  the  interaction,  this  impulse  may  be  sufficient 
either  to  raise  the  electron  to  a  higher-lying  shell  within  the  absorber  (excitation) 
or  to  completely  remove  the  electron  from  the  atom  (ionization).  In  either  case,  the 
energy  that  is  transferred  to  the  electron  must  come  at  the  expense  of  the  charged 
particle,  and  its  velocity  is  therefore  decreased  as  a  result  of  the  encounter.  At 
any  given  time,  a  charged  particle  interacts  with  many  electrons,  so  the  net  effect 
is  to  decrease  its  velocity  continuously  until  the  particle  is  stopped  at  its  depth  of 
penetration  in  the  material,  or  range. 

The  alpha  particle  energies  mentioned  in  the  previous  paragraph  have  approxi¬ 
mate  ranges,  in  air,  of  less  than  4  cm  [12].  For  materials  with  higher  atomic  numbers, 
that  value  is  significantly  reduced;  for  example,  the  same  particles  would  travel  less 
than  100  fin i  in  silicon  [12].  Therefore,  the  primary  decay  mechanism  for  common 
nuclear  weapon  isotopes  is  easily  shielded  and  does  not  facilitate  successful  identifi¬ 
cation  of  special  nuclear  material. 

As  mentioned  earlier,  a  product  isotope  may  be  in  an  excited  or  meta-stable 
state  following  nuclear  decay,  and  subsequent  decay  to  a  stable  ground  state  is  often 
accompanied  by  gamma  ray  emission.  As  with  alpha  particle  decay,  distributions 
of  characteristic  gamma  ray  energies  are  well-known  and  tabulated.  For  example, 
the  approximate  ranges  of  characteristic  gamma  ray  energies  associated  with  235U 
and  239Pu  are  110-280  keV  and  40-400  keV,  respectively  [11],  Given  these  energy 
ranges,  the  dominant  photon  interaction  mechanisms  are  the  photoelectric  effect  and 
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Compton  scattering,  whereby  electrons  receive  the  partial  or  complete  transfer  of  the 
gamma  ray  photon  energy.  The  resulting  interactions  exponentially  attenuate  the 
gamma  ray  photons,  as  functions  of  photon  energy,  interaction  material  density,  and 
interaction  material  thickness.  Specifically,  a  narrow  beam  of  monoenergetic  photons 
with  an  incident  intensity  Jo,  penetrating  a  layer  of  material  with  thickness  x  and 
density  p,  will  emerge  with  intensity  /  given  by  the  exponential  attenuation  law 

I(x)  =  I0e-(p)pxy  (5) 

where  p/ p  is  the  photon  mass  attenuation  coefficient,  which  is  a  function  of  incident 
photon  energy  and  absorber  material  density.  These  experimentally  determined  co¬ 
efficients  are  well-known  and  tabulated.  As  an  example,  the  National  Institute  of 
Standards  and  Technology  (NIST)  provides  a  comprehensive,  electronic  listing  [13]. 

Using  (5)  and  appropriate  values  from  the  NIST  database,  the  highest  gamma 
ray  energies  associated  with  SNM  isotopes  listed  in  Table  1  require  less  than  a  3  cm 
thickness  of  lead  shielding  to  attenuate  99.9%  of  emitted  gamma  ray  photons;  there¬ 
fore,  gamma  rays  are  easily  shielded.  Furthermore,  distinguishing  SNM  emitted 
gammas  from  the  background  continuum  of  gamma  rays  and  X-rays  from  naturally 
occuring  terrestrial  and  cosmic  radiation  sources  compounds  the  problem.  Hence, 
identification  of  special  nuclear  material  from  either  characteristic  alpha  decay  or 
characteristic  gamma  ray  emission  is  problematic. 

Electrons  and  positrons  are  primarily  created  following  the  SNM  fission  product 
beta  (/?)  decay  mechanisms  of 


n 


P  +  P  +V, 


(6) 


and 


p — >n  +  (3+  +  v,  (7) 

where  u  and  v  are  the  antineutrino  and  neutrino  particles.  The  beta  decay-produced 
charged  particles  offer  little  potential  for  identifying  SNM.  Analysis  of  the  particles’ 
interactions  with  matter  shows  that  collisional  and  radiative  energy  losses  yield  beta 
particle  ranges  that  are  approximately  100  times  larger  than  those  for  an  energeti¬ 
cally  equivalent  alpha  particle,  but  are  still  easily  attenuated  by  minimal  shielding 
material.  Thus,  neutron  emission  detection  is  the  remaining  choice  to  identify  SNM 
via  radiation  detection. 

1.2.3  Non-ionizing  Radiation  Interactions  with  Matter. 

Due  to  their  charge  neutrality,  neutrons  interact  with  the  nucleus,  as  opposed 
to  with  the  atomic  electrons,  and  the  mechanisms  that  govern  their  attenuation  in 
a  material  differ  slightly  from  those  associated  with  ionizing  radiation.  Typically,  a 
neutron  will  lose  energy  in  material  through  a  series  of,  for  the  most  part,  elastic 
scattering  events.  As  neutron  energy  decreases,  scattering  continues,  but  the  likeli¬ 
hood  of  its  capture  by  a  nucleus  generally  increases.  First,  we  shall  consider  neutron 
production  by  special  nuclear  material,  followed  by  a  more  detailed  discussion  of 
neutron  interaction  mechanisms. 

Recall  that  Table  1  listed  spontaneous  fission  neutron  emission  rates  for  common 
nuclear  weapon  materials.  Spontaneous  fission  is  a  form  of  radioactive  decay  char¬ 
acteristic  of  very  heavy  isotopes.  As  the  atomic  number,  Z,  and  atomic  mass,  A, 
of  an  isotope  increase,  the  sensitive  balance  between  the  nuclear  binding  energy  and 
the  Coulombic  repulsion  between  the  protons  within  the  nucleus  may  be  disrupted, 
at  which  point  it  becomes  energetically  possible  for  the  nucleus  to  spontaneously 
divide  into  two  smaller  fragments  (fission).  A  quantum  mechanical  analysis  of  the 
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semi-empirical  mass  formula  (SEMF)  for  nuclear  binding  energy  dictates  that  this 
phenomenon  occurs,  approximately,  for  isotopes  in  the  limiting  condition  [14] 


Z2 

A 


2  ds 
CLc 


49, 


(8) 


where  as  and  ac  are  the  surface  and  Coulomb  SEMF  correction  factors  to  the  liquid 
drop  binding  energy  model.  The  spontaneous  fission  neutron  emission  rate,  hsF, 
listed  in  Table  1,  is  calculated  as 


nsF 


f  mNAV\ 
V  Maw  ) 


'P SF  Tl  fission  j 


(9) 


where  ti/2  is  the  isotope  half-life,  m  is  the  isotope  mass,  NAv  is  Avogadro’s  number, 
Maw  is  the  isotope  molar  mass,  Vsf  is  the  isotope  probability  of  spontaneous  fission, 
and  h fiSSion  is  the  average  neutrons  emitted  per  fission  event.  These  data  are  accessible 
from  NNDC,  with  the  exception  of  the  n fission  data,  which  are  indicated  in  the  Table  1 
footnote  via  [9]. 

In  addition  to  spontaneous  fission  neutrons,  neutrons  may  be  liberated  by  fission 
that  is  induced  from  neutrons  originating  from  terrestrial,  cosmic,  or  other  radioac¬ 
tive  sources.  In  general,  the  prompt  fission  responses  of  the  SNM  of  interest  follow 
the  reactions 


zX  +  Qii  — >■  A+zX*  — >■  Fission  Fragments  (2)  +  gii  (2  or  3)  +  7  +  Q .  (10) 

As  with  alpha  particle  and  gamma  ray  decay,  the  reaction  energies  and  their  dis¬ 
tributions  are  tabulated  for  the  nuclear  weapon  material  isotopes,  although  explicit 
calculation  of  the  reaction  energies  is  accomplished  using  the  mass  energies  of  the 
product  and  reactants,  as  done  in  (4). 

It  should  be  noted  that  the  fission  fragment  pair  produced  in  reaction  (10)  does 
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not  always  consist  of  the  same  nuclides,  even  for  the  case  of  induced  fission  with  iden¬ 
tical  target  nuclei  and  identical  bombarding  neutron  energy.  Instead,  a  range  of  pairs 
is  possible,  subject  to  the  requirement  that  the  total  number  of  nucleons  is  conserved, 
although  experimentally  determined  maxima  exist  at  A  =  95  and  A  =  140  for  235U. 
Further  experimental  analysis  indicates  that  the  average  value  of  prompt  fission  en¬ 
ergy  for  all  possible  fission  fragment  pairs  is  180  MeV  for  235U  and  185  MeV  for  239Pu 
and  that  these  energies  are  distributed  among  the  products  according  to  Table  2. 
Therefore,  the  average  number  of  neutrons  emitted  per  fission  event  (2-3)  from  re- 


Table  2.  Prompt  fission  energy  distribution  [9]. 


Energy  Form 

235-(j 

239Pu 

Fission  Fragment  Kinetic  Energy 

168 

172 

Neutron  Kinetic  Energy 

5 

6 

Prompt  Gamma  Energy 

7 

7 

Total  Prompt  Energy  (MeV) 

180 

185 

action  (10)  and  the  neutron  kinetic  energy  allocation  from  Table  2  indicates  that 
the  average  neutron  energy  resulting  from  235U  and  239Pu  fission  is  approximately 
2  MeV;  however,  the  most  probable  neutron  energy  is  approximately  1.5  MeV  [9]. 

The  interaction  of  these  neutrons  with  matter  and  the  likelihood  of  identifying 
SNM  via  neutron  detection  must  be  examined  next.  A  neutron’s  ability  to  generate 
an  electronic  signal  in  a  detector  requires  conversion  of  the  neutron  into  a  charged 
particle  through  a  secondary  reaction,  and  the  likelihood,  or  probability,  of  this 
conversion  is  governed  by  nuclear  cross  sections. 

The  concept  of  a  microscopic  nuclear  cross  section  cr  is  that  the  effective  size  of  a 
nucleus  should  be  proportional  to  the  probability  that  an  incident  particle  would  react 
with  it.  This  probability  is  expressed  conventionally  in  terms  of  the  barn  (10_28m2) 
for  each  type  of  interaction.  For  example,  each  isotope  will  have  an  elastic  scattering 
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cross  section,  a  radiative  capture  cross  section,  and  so  on,  each  of  which  will  be  a 
function  of  the  neutron  energy.  The  hierarchy  of  microscopic  cross  sections  used  to 
determine  neutron  interaction  rates  is  shown  in  Figure  1,  courtesy  of  [15]. 
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Figure  1.  Microscopic  cross  section  hierarchy  used  to  calculate  neutron  interaction 
rates.  Adapted  from  [15]. 


When  multiplied  by  the  the  number  of  nuclei  per  unit  volume  N ,  the  microscopic 
cross  section  is  converted  to  a  macroscopic  cross  section,  E,  then  as 

E  =  N  a  ,  (11) 

which  has  dimensions  of  inverse  length,  and  represents  physically  the  probability  per 
unit  path  length  for  the  specific  process  described  by  the  microscopic  cross  section. 
Combining  the  processes  for  each  interaction 

E  total  Escalt,e|.  T  Erad.  capture  T 

gives  the  resulting  total  macroscopic  cross  section,  or  the  probability  per  unit  path 
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length  that  any  interaction  type  will  occur.  Consequently,  the  neutron  flux-dependent 
volumetric  reaction  rate,  or  reaction  rate  density  Rvo\  is  directly  proportional  to  £ 
and  is  calculated  as 

Rv oi  =  Sd> ,  (13) 

where  <f>  is  the  neutron  flux. 

As  with  photons,  and  in  a  simplistic,  one-dimensional  analysis,  if  a  narrow  beam 
neutron  experiment  were  conducted,  the  entering  neutron  flux  <h0  would  attenuate 
exponentially  [12, 15]  with  absorber  thickness  x  as 

<h(a;)  =  $0e's,ota11 .  (14) 

Rigorous  neutron  attenuation  analysis  requires  replacing  the  simple  position-dependent 
neutron  flux  $(a;)  considered  in  (14)  with  a  multivariable  counterpart  described  by 
the  Boltzmann  neutron  transport  equation.  However,  full  analysis  of  the  Boltzmann 
equation  is  unwarranted  here,  as  a  simplistic  treatment  of  neutron  flux  yields  a  similar 
conclusion.  That  is,  analysis  of  (14)  with  appropriate  cross  section  data  would  show 
that  neutrons  produced  by  fission  reactions  possess  sufficient  energy  to  escape  most 
typical  shielding  materials  with  minimal  attenuation  such  that  the  neutron  flux  value 
resulting  from  the  assumptions  of  an  unmoderated,  unshielded  SNM  mass,  discussed 
in  Section  1.2.1,  is  a  reasonable  initial  estimate. 

Therefore,  neutrons  offer  the  most  accurate  and  efficient  method  to  identify  spe¬ 
cial  nuclear  material.  However,  while  SNM  identification  via  neutron  detection  ap¬ 
pears  favorable  from  a  physical  standpoint,  several  engineering  design  considerations 
must  be  addressed. 
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1.3  General  Neutron  Detector  Design  Considerations 


A  neutron’s  ability  to  generate  an  electronic  signal  in  a  detector  requires  conver¬ 
sion  of  the  neutron  into  a  charged  particle  through  a  secondary  reaction,  as  shown 
below  [12]. 


Target  Nucleus  +  neutron 


Recoil  Nucleus 
Fission  Fragments 
<  Alpha  Particle 
Proton 
Electron 


However,  these  neutron  induced  secondary  reactions  occur  only  with  appreciable 
probability  (cross  section)  in  the  thermal  neutron  energy  range,  about  0.025  eV  at 
room  temperature.  Therefore,  accurate  and  efficient  neutron  detection  requires  a 
careful  combination  of  moderator,  conversion,  and  detector  materials,  as  shown  hy¬ 
pothetically  in  Figure  2,  that  successfully  reduces  the  energy  of  a  fast  neutron  to 


Figure  2.  General  schematic  (top  view)  of  a  hypothetical  solid  state  semiconducting 
device  capable  of  neutron  detection. 


the  thermal  range,  converts  thermal  neutrons  to  charged  particles,  and  responds 
favorably  to  this  ionizing  radiation  by  emitting  and  collecting  charged  particles  or 
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photons. 

Neutron  opacity,  or  complete  neutron  capture,  within  the  conversion  volume  or 
layer  presents  one  design  challenge.  If  the  moderator  and  conversion  layers  are  too 
thin,  neutrons  are  not  captured  and  converted  to  charged  particles;  if  the  layers  are 
too  thick,  the  reaction  products  fail  to  reach  the  active  layer  of  the  detector.  Simi¬ 
larly,  the  distance  traveled  by  the  charged  particle  reaction  products  has  important 
consequences  within  the  detector.  Specifically,  in  order  to  capture  the  full  kinetic 
energy  of  the  products,  the  detector  must  be  designed  with  an  active  volume  that  is 
large  enough  to  fully  stop  the  products. 

Discrimination  between  neutron-induced  and  gamma-induced  reactions  in  a  de¬ 
tector  is  another  consideration.  If  the  active  layer  of  the  detector  is  much  smaller 
than  the  mean  free  path  of  gamma  rays  in  the  material,  then  the  detector  would 
be  insensitive  to  these  photons.  However,  if  the  gamma  mean  free  path,  which  is  a 
function  of  material  density  and  Z,  or  the  number  of  electrons  per  atom,  approaches 
the  active  layer  width  of  the  detector,  the  neutron  and  gamma  signals  may  overlap, 
and  detection  via  gross  neutron  counting  can  be  problematic. 

Of  course,  efficient  neutron  capture  and  conversion  are  inadequate  if  the  neutron 
conversion  reaction  products  cannot  produce  electron-hole  pairs  for  efficient  collec¬ 
tion.  In  general,  the  charge  created  from  a  single  ionizing  radiation  interaction  is 
proportional  to  the  energy  lost  by  an  ionizing  agent  Eion  in  the  material,  at  the 
cost  of  approximately  3 Eg  per  electron-hole  pair  [16],  where  Eg  is  the  bandgap  en¬ 
ergy.  Assuming  full  energy  deposition  of  the  ionizing  radiation,  the  maximum  charge 
created  in  a  single  event  is  approximated  as 

Qmax  =  X  (1.6  X  10-19)  [C]  .  (15) 

SEg 

The  delicate  balance  of  material  composition  and  dimensions  must  be  coupled 
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with  a  device  design  that  maximizes  the  electronic  signal,  such  that  the  generated 
charge  exceeds  the  intrinsic  noise  level  of  the  charge-sensitive  commercial  preamplifier- 
amplifier  combination.  In  short,  the  cross  section  for  the  neutron  conversion  reaction 
must  be  as  large  as  possible  so  that  efficient  detectors  can  be  built  with  small  di¬ 
mensions.  If  the  detection  medium  is  a  solid,  this  requirement  is  typically  achieved 
because  the  ranges  of  the  ionizing  radiation  reaction  products  in  solid  materials  are 
on  the  order  of  103  less  than  products  of  equivalent  energy  in  a  gas  detection  medium. 
Hence,  the  ranges  of  these  products  in  gas-based  detection  systems  can  be  significant 
compared  to  the  detector  dimensions  and  can  prevent  full  collection  of  the  products. 

Common  neutron  detection  materials  include  boron-10  (10B),  lithium-6  (6Li),  and 
helium-3  (3He),  which  are  configured  in  designs  that  include  proportional  counters, 
scintillators,  thermoluminescent  dosimeters,  solid  state  conversion  layers,  and  neu¬ 
tron  absorbing  solid  state  layers.  Proportional  counters  employ  gas  as  the  detection 
medium  and,  by  design,  tend  to  be  large  compared  to  the  remaining  designs,  which 
use  solid  detection  materials.  Furthermore,  limited  production  of  3He,  caused  in 
combination  by  the  1988  closure  of  the  Savannah  River  production  facility  and  the 
post-9/ 11  growth  in  security- related  neutron  detectors  [17],  has  elevated  alternatives 
to  3He  as  a  top  priority  for  the  Department  of  Homeland  Security  [17]. 

Currently,  most  commercially  available  solid  state  neutron  detectors  are  based  on 
lithium  iodide  scintillation.  Scintillators  and  thermoluminescent  dosimeters  (TLD), 
shown  in  Figures  3(a)  and  3(b),  respectively,  consist  of  neutron  absorbing  materials 
that  decay  into  ionizing  radiation  reaction  products  and  create  electron-hole  pairs. 
Scintillation  detectors  and  thermoluminescent  dosimeters  are  similar  in  that  both 
use  solid  state  materials  that  incorporate  dopants,  which  capture  the  electrons  and 
holes  and  then  measure  the  light  output,  via  a  photomultiplier  tube  (PMT),  when 
the  electrons  and  holes  recombine.  However,  while  scintillators  facilitate  rapid  re¬ 
combination  of  electrons  and  holes,  thermoluminescent  dosimeters  require  heating  in 
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order  to  de-trap  the  electron-hole  pairs. 


Conduction  band 


Conduction  band 
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(a)  Scintillator. 


(b)  Thermoluminescent  dosimeter. 


Figure  3.  Scintillation  detector  schematics. 


Solid  state  neutron  detectors  may  also  take  the  form  of  p  —  n  junction  or  Schot- 
tky  heterostructures,  and  are  separated  into  two  subclasses.  Indirect  conversion, 
or  thin-film-coated  [18,19],  neutron  detectors,  as  depicted  in  Figure  4(a),  use  neu¬ 
tron  absorbing  material  to  generate  ionizing  radiation  reaction  products  which  in¬ 
teract  with  adjacent  conventional  semiconductor  devices.  Clever  planar  and  stack 
designs  are  used  to  maximize  the  surface  area  between  conversion  materials  and  space 
charge.  Direct  conversion,  or  solid-form  [18, 19]  neutron  detectors,  shown  in  Figure 
4(b)  are  fabricated  from  semiconductors  doped  with  neutron  conversion  material. 
Unconverted  reaction  products  are  minimized,  and  higher  theoretical  efficiencies  are 
possible.  While  clever  designs  are  not  required,  maturity  of  materials  and  an  under¬ 
standing  of  the  physical  and  electronic  properties  of  the  heterostructure  materials 
are  the  limiting  factors  in  direct  conversion  designs  [18-20] . 

Figure  5  shows  the  detector  efficiencies  for  several  common,  thin-hlm  semicon¬ 
ductor  diodes,  and  indicates  that  the  most  efficient  diode,  containing  an  enriched 
10B  layer,  requires  an  active  volume  thickness  of  approximately  40  /irn  for  optimal 
performance,  whereas  the  other  designs  require  thicknesses  between  several  hundred 
microns  and  over  one  centimeter  to  achieve  comparable  performance.  Device  de¬ 
signs  that  minimize  the  conversion  layer  thickness  are  important  for  several  reasons. 
First,  the  ranges  and  mean  free  paths  of  charged  particles  and  photons  created  from 
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(a)  Indirect  conversion.  (b)  Direct  conversion. 

Figure  4.  Solid  state  neutron  conversion  detector  schematics.  Adapted  from  [19]. 


neutron  conversion  decrease  with  increased  thickness.  Thus,  charge  collection  effi¬ 
ciency  is  potentially  reduced  and  the  ability  to  discriminate  between  neutron  and 
gamma  interactions  is  complicated.  Second,  thicker  hints  are  more  costly  and  poten¬ 
tially  more  difficult  to  grow  and  fabricate  into  devices,  as  compared  to  thinner  hints. 
Thin-hlm-coated  semiconductor  devices  have  been  investigated  for  thermal  neutron 
detection  by  various  research  groups  [18,21-40],  all  of  which  have  used  10B,  6Li,  6LiF, 
Cd,  U,  or  Gd  coatings  as  the  neutron  reactive  layer,  although  the  most  successful 
conversion  layer  devices  are  based  on  10B  enriched  boron  coatings. 


Detector  Thickness  (cm) 


Figure  5.  Neutron  absorption  (detection  and  capture  efficiencies)  as  a  function  of 
detector  thickness. 


As  shown  in  Figure  6,  gadolinium-157  (15'Gd)  possesses  a  thermal  neutron  cap- 
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ture  cross  section  that  is  nearly  two  orders  of  magnitude  larger  than  that  of  10B,  and 
is  the  highest  of  all  known  isotopes,  thus  making  it  a  potentially  favorable  choice  for 
neutron  detection.  Moreover,  due  to  the  unique  energy  dependence  of  its  cross  sec¬ 
tion,  157Gd  maintains  its  utility  for  thermal  neutron  capture  and  detection  at  neutron 
energies  in  the  range  of  200-300  meV,  as  compared  with  approximately  25  meV  for 
the  common  neutron  absorbers.  Finally,  the  secondary  conversion  reactions  result 
in  the  emission  of  ionizing  radiation,  to  include  gamma-ray  and  characteristic  X-ray 
photons,  as  well  as  Auger  and  internal  conversion  electrons.  When  absorbed  by  a 
detector  material,  these  ionizing  radiation  photons  and  charged  particles  can  ulti¬ 
mately  provide  reliable  indication  of  the  initial  neutron  absorption.  From  a  neutron 


Figure  6.  Energy  dependent  neutron  capture  cross  sections  for  various  materials. 


detector  device  standpoint,  these  favorable  properties  indicate  that  neutron  opacity 
in  Gd  occurs  for  smaller  film  thicknesses  than  in  the  other  materials  shown  in  Figure 
6.  ft  follows  that  the  thinner  Gd  films  are  transparent  to  photons  at  a  lower  en¬ 
ergy  threshold  than  the  remaining  materials.  This  favorable  combination  of  neutron 
opacity  and  gamma  discrimination  for  Gd  films  offers  the  likelihood  that  a  signal 
generated  from  a  Gd-based  neutron  detector  could  be  definitively  correlated  with 
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neutron  emission  from  a  radioactive  source,  as  opposed  to  gamma  decay.  Several 
researchers  [21-23]  have  incorporated  gadolinium  in  semiconductor  device  neutron 
detectors  with  varying  results.  However,  there  remains  great  opportunity  for  contin¬ 
ued  exploration  and  development  of  high-quality  Gd  doped  semiconductor  materials 
and  neutron  detector  device  architectures. 

1.4  Research  Objectives  and  Overview 

Semiconductors  that  incorporate  large  neutron  absorption  cross  section  materi¬ 
als,  such  as  157Gd,  are  gaining  attention,  despite  the  relative  immaturity  of  their 
development.  Gallium  nitride  (GaN)  is  one  semiconductor  material  which,  by  it¬ 
self,  has  attained  considerable  recognition  over  the  past  two  decades.  Its  favorable 
properties  as  a  device  material  include  a  large  bandgap  (~  3.4  eV),  thermal  and 
chemical  stability,  and  high  carrier  mobility,  ensuring  that  GaN-based  designs  are 
robust  in  the  presence  of  fluctuating  temperature  and  ionizing  radiation  environ¬ 
ments.  In  the  past,  successful  incorporation  of  rare  earth  (RE)  dopants,  to  include 
Gd,  into  GaN  structures  had  been  limited,  but  during  the  past  decade,  rare  earth 
doped  semiconductors  have  generated  considerable  attention  for  their  application  in 
new  optoelectronic  devices  [41-44],  Despite  the  increased  attention,  incomplete  un¬ 
derstanding  of  the  physical  and  electronic  properties  of  these  materials  classifies  rare 
earth  doped  GaN  heterostructures  as  non-conventional,  as  compared  to  silicon-based 
devices.  It  is  therefore  of  considerable  interest  to  know  whether  even  low  concentra¬ 
tions  of  a  rare  earth  in  a  GaN  host  can  alter  the  surface  electronic  structure.  With 
these  considerations  in  mind,  investigations  of  the  surface  electronic  structure  and 
interface  properties  of  RE^Gai-^N  (RE  =  Gd,  Er,  Yb)  semiconductors  are  neces¬ 
sary  in  order  to  advance  the  scientific  community’s  knowledge,  and  to  determine  the 
viability  of  these  materials  in  neutron  detection  heterostructures. 

To  that  end,  systematic  studies  of  these  rare  earth  doped  GaN  thin  films  were 
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conducted  over  a  period  of  20  months  in  this  research.  Specifically,  the  research  was 
designed  to  answer  several  questions,  to  include: 

(1)  Can  rare  earth  dopants  be  incorporated  into  a  GaN  lattice  using  various  epi¬ 
taxial  growth  techniques? 

(2)  Is  the  location  of  the  rare  earth  dopant  in  the  GaN  lattice  classified  as  substi¬ 
tutional  or  interstitial? 

(3)  Where  are  the  rare  earth  dopant  4/  electronic  levels  located  within  the  valence 
band  of  rare  earth  doped  GaN  thin  films? 

(4)  How  does  the  Schottky  barrier  height  at  the  metal  to  rare  earth  doped  GaN 
thin  film  interface  compare  with  undoped  GaN  interfaces? 

(5)  What  does  fundamental  nuclear  and  semiconductor  physics  predict  about  the 
viability  of  gadolinium  as  a  neutron  detection  material? 

First,  temperature  dependent  X-ray  photoemission  spectroscopy  (XPS)  studies 
of  Yb  doped  GaN  (GaN:Yb)  were  conducted  in  order  to  determine  the  effective 
surface  Debye  temperature,  and  concluded  that  the  rare  earth  dopant  occupied  a 
substitutional  Ga  site  within  the  GaN  lattice,  as  had  been  previously  reported  [41,42], 
Armed  with  confidence  that  the  epitaxial  growth  process  produced  thin  him  sam¬ 
ples  with  rare  dopants  incorporated  successfully  into  the  GaN  lattice,  resonant  pho¬ 
toemission  studies  of  several  rare  earth  dopants  (Gd,  Er,  Yb)  in  GaN  were  conducted 
next.  The  key  value  of  resonant  photoemission  is  to  probe  what  valence  bands  of 
the  semiconductors  have  strong  4/  and/or  rare  earth  weight  [35,45,46].  The  value 
in  studying  the  resonant  photoemission  process  of  a  semiconductor  with  several  dif¬ 
ferent  rare  earths  is  that  while  each  is  likely  to  dope  the  semiconductor  in  a  similar 
fashion,  the  resonant  enhancement  of  the  valence  band  will  occur  at  different  pho¬ 
ton  energies.  These  studies  concluded  that  the  placement  of  the  occupied  Gd,  El¬ 
and  Yb  4/  levels  were  more  at  the  bottom  of  the  valence  band  than  predicted  by 
theory.  This  suggests  that  photons  emitted  from  intra-atomic  /-/  transitions  may 
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be  more  energetic  or  ‘blue’  than  predicted  by  many  theoretical  models,  which  has 
implications  for  the  physics,  device  fabrication,  and  engineering  communities. 

Next,  device  fabrication  and  neutron  detector  applications  were  explored  more 
carefully  by  measuring  the  Schottky  barrier  heights  at  the  gold  to  rare  earth  doped 
GaN  interface  for  the  Gd,  Er,  and  Yb  dopants  via  ultraviolet  photoemission  spec¬ 
troscopy.  The  value  in  studying  Schottky  barrier  heights  is  as  follows.  Heterostruc¬ 
ture  diodes  that  function  as  radiation  detectors  are  reverse  biased  to  minimize  current 
flow  through  the  device;  thus,  current  is  limited  to  leakage  only.  Under  reverse  bias, 
the  depletion  region,  which  correlates  directly  with  the  radiation-induced  charge 
collection  volume,  widens  with  increasing  reverse  voltage.  Despite  efforts  to  limit 
current  flow  due  to  leakage,  and  maximize  the  charge  collection  volume,  the  gener¬ 
ated  signal  due  to  radiation- induced  currents  can  prove  difficult  to  detect  above  noise 
due  to  leakage  currents.  To  alleviate  this  obstacle,  blocking  or  Schottky  contacts  are 
often  used  to  reduce  leakage  currents  and  enhance  the  signal-to-noise  ratio  in  the 
radiation  detector  device  architecture. 

The  completed  Schottky  barrier  studies  discovered  larger  (25-55%)  barrier  heights 
at  the  metal  to  rare  earth  doped  GaN  interface  than  at  interfaces  with  undoped 
GaN.  These  results  were  uniform  for  all  of  the  rare  earths  studied,  and  imply  further 
reduction  in  leakage  currents  as  compared  with  undoped  GaN  devices.  Moreover, 
the  larger  barrier  height  at  the  metal  to  rare  earth  doped  GaN  interface  means  that 
radiation-induced  charge  collection  volumes  are  increased,  and  thus  more  efficient, 
over  what  would  be  the  case  in  pure  GaN  devices. 

Finally,  the  utility  of  gadolinium  as  a  neutron  detection  material  in  a  hypothet¬ 
ical  direct-conversion  p  —  n  homojunction  diode  was  examined  using  fundamental 
nuclear  and  semiconductor  physics.  The  charge  produced  in  a  Gd-based  semicon¬ 
ductor  by  an  internal  conversion  electron  created  from  157Gd  neutron  capture  was 
found  to  approach  the  intrinsic  noise  level  of  commercial  preamplifiers  required  for 


22 


pulse  height  analysis.  Neutron  capture  efficiency  was  shown  to  be  maximized  with 
Gd-layer  thickness,  and  the  bias-dependent  depletion  width  of  the  Gd-based  diode 
to  the  full  range  of  the  internal  conversion  electrons  (~  25  /mi)  was  considered  for 
low  power  operation. 

These  preliminary  results  are  promising,  but  more  work  is  needed  to  fully  under¬ 
stand  the  material  interfaces  and  their  implications  for  full  device  fabrication. 
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II.  Theory 


This  chapter  discusses  theory  pertaining  to  the  experimental  techniques  used  dur¬ 
ing  research,  beginning  with  a  general  overview  of  photoelectron  spectroscopy.  Subse¬ 
quent  sections  present  applications  of  the  theory,  including  temperature-dependent 
photoemission,  resonant  photoemission,  and  the  photoemission  response  to  metal 
overlayer  deposition.  Each  technique  was  used  during  the  research  owing  to  its  valu¬ 
able  insight  to  surface  electronic  structure.  As  devices  continue  to  scale  smaller,  and 
as  rare  earth  doped  GaN  optoelectronic  devices  continue  to  gain  popularity,  com¬ 
prehensive  understanding  of  their  material  surface  properties  becomes  increasingly 
valuable  to  efficient  device  design  and  engineering. 

2.1  Photoelectron  Spectroscopy  Models 

In  a  general  sense,  photoelectron  spectroscopy  (PES)  provides  information  about 
occupied  states  in  an  atom,  molecule,  or  solid  in  that  the  spectroscopic  techniques 
directly  measure  a  system’s  orbital  energies.  The  theoretical  basis  behind  these  tech¬ 
niques  is  photoemission,  described  by  the  photoelectric  effect,  which  was  first  identi¬ 
fied  in  1887  by  Hertz  [2]  and  subsequently  explained  by  Einstein  [3].  Photoemission, 
depicted  for  an  arbitrary  solid  in  Figure  7,  results  from  incident  photon  annihilation, 
which  excites  an  electron  from  a  bound  state  of  a  system  (atom,  molecule,  solid) 
to  a  continuum  state,  such  as  the  vacuum.  For  the  process  to  occur,  the  photon 
energy  hv  must  be  sufficient  to  provide  the  electron  with  enough  energy  to  overcome 
the  work  function  cj)  and  the  binding  energy  Eb  of  the  bound  state,  such  that,  in  a 
vacuum,  the  kinetic  energy  Ekin  of  a  photoelectron  is  measured  as 

-Ekin  =  hv  -  \Eb\  -  (/>■ 
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(16) 


Figure  7.  The  photoelectron  energy  distribution  produced  by  incoming  photons  and 
measured  as  a  function  of  the  kinetic  energy  -Ekin  of  the  photoelectrons  (right)  is  more 
conveniently  expressed  in  terms  of  the  binding  energy  Eb  (left)  when  one  refers  to 
the  density  of  states  inside  the  solid  (Eb  =  0  at  Ep.)  Adapted  from  [1]. 


Prior  to  discussing  the  photoemission  techniques  relevant  to  this  research,  an 
overview  of  the  two  generally  accepted  models  of  photoemission  is  required.  The 
one-step  model  and  three-step  model  explain  the  phenomenon  in  either  a  single  or  a 
three  step  process,  as  their  names  indicate. 

The  three-step  model,  developed  by  Berglund  and  Spicer  [4],  adequately  and  ped- 
agogically  explains  the  photoemission  process  as  three  independent  and  sequential 
processes.  In  the  first  step,  the  photon  is  absorbed  locally  and  an  electron  is  excited. 
In  the  second  step,  this  electron  travels  through  the  sample  to  the  surface.  Finally, 
in  the  third  step  the  electron  escapes  through  the  surface  into  the  vacuum  where  it 
is  detected.  As  expected,  the  one-step  model  treats  the  entire  process  as  a  single 
coherent  step,  but  this  model  is  mathematically  rigorous  and  less  conducive  to  a 
qualitative  understanding  of  the  photoemission  process  than  its  multi-step  counter¬ 
part.  Furthermore,  although  the  separation  of  events  is  somewhat  artificial  in  that 
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the  entire  process  should  be  treated  as  a  single  step,  the  division  of  the  process  into 
three  distinct  steps  yields  results  that  are  not  very  different  from  the  conceptually 
unwieldy  one-step  model  [1],  Thus,  given  the  adequacy  of  a  three-step  model  analysis 
and  the  unnecessary  rigor  of  a  one-step  model  examination,  a  detailed  discussion  of 
the  one-step  model  is  omitted  in  the  interests  of  clarity  and  simplicity. 

In  addition  to  the  three-step  model’s  treatment  of  photoemission  as  sequential, 
independent  steps,  the  model  also  assumes  that  the  energy  distribution  of  pho¬ 
toelectrons  N(E,hu)  results  from  a  distribution  of  primary  unscattered  electrons 
Npr[(E,  hv)  and  a  distribution  of  background  secondary  electrons  Nsec(E,  hv)  that 
result  from  one  or  more  scattering  events,  or 

N(E,  hv)  =  Npri(E ,  hv )  +  Nsec{E,  hv) .  (17) 

Furthermore,  the  primary  photoelectron  energy  distribution  results  from  the  convo¬ 
lution  of  three  independent  functions,  namely  the  transmission  function,  the  escape 
function,  and  the  internal  energy  distribution  of  photoelectrons  N-mt  [5] ,  each  of  which 
will  be  discussed  in  the  following  sections. 

2.1.1  Optical  Excitation  of  the  Electron  in  the  Solid. 

The  following  three-step  model  theory  is  summarized  from  Hiifner  [1],  The  pho¬ 
tocurrent  produced  in  a  PES  experiment  results  from  the  excitation  of  electrons  from 
an  initial  energy  eigenstate  with  wave  function  'It,  to  a  continuum  of  final  states  with 
wave  function  T f  by  a  photon  field  having  the  vector  potential  A,  and  may  be  ex¬ 
plained  using  time-dependent  perturbation  theory.  Assuming  a  small  perturbation 
H' ,  the  transition  probability  per  unit  time  w  between  the  initial  eigenstate  dq  and 
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the  final  electron  state  T/  is  given  by  Fermi’s  golden  rule  as 

w<x^\(*f\H,\*i)\26(Ef-Ei-hv)  (18) 

where  h  is  Planck’s  constant  divided  by  2n,  | (T /  \H'\  'F*)|2  is  the  square  of  the  tran¬ 
sition  matrix  element  of  the  time-dependent  perturbing  Hamiltonian  operator  H1, 
and  the  Kronecker  delta  function  S(Ef  —  E{  —  hu)  dictates  that  the  initial  and  final 
states  of  the  system  must  be  energetically  identical  to  ensure  a  non-zero  transition 
probability.  It  should  be  noted  that  w  must  vary  slowly  enough  to  determine  time 
intervals  that  satisfy  the  condition  of  time  dependence,  yet  small  enough  to  classify 
H'  as  a  small  perturbation  [6].  Thus,  an  appropriate  expression  for  the  perturbing 
Hamiltonian  H'  is  required  to  determine  the  transition  rate  and,  hence,  the  internal 
energy  distribution  of  photoexcited  electrons. 

In  the  most  general  form,  the  Hamiltonian  H'  of  an  electron  in  an  arbitrary 
electromagnetic  field  is  given  by  [6] 

H>  =  2^  '  P  +  P  '  A)  +  G)  (A  '  A^]  “  61P  ’ 

where  e  and  m  are  the  electron  mass  and  charge,  respectively,  c  is  the  speed  of  light 
in  a  vacuum,  A  is  the  vector  potential,  p  is  the  momentum  operator  p  =  ihV,  and 
( p  is  the  scalar  potential;  the  scalar  potential  ip  is  not  to  be  confused  with  the  work 
function  (f).  Using  the  commutation  relation 

[p,  A]  =  (p  •  A)  -  (A  •  p)  =  — i h(V  ■  A) ,  (20) 
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the  first  two  terms  in  equation  (19)  reduce  as 


(A  ■  p)  +  (p  •  A)  =  (p  •  A)  -  (A  •  p)  +  2(A  •  p) 

=  -ift(V  •  A)  +  2(A  •  p) , 


and  H'  can  be  re-written  as 


H' 


e 

2  me 


2(A  p) 


-  ift(V  ■  A)  + 


(A -A) 


(21) 

(22) 


(23) 


The  perturbing  Hamiltonian  H'  in  (23)  can  be  simplified  further  by  selecting 
the  Coulomb  gauge  such  that  V  •  A  =  0,  which  is  a  valid  assumption  because  of 
translational  invariance  in  the  solid  [1].  A  more  physically  intuitive  justification 
would  be  that  A  is  constant  over  the  atomic  dimensions  of  the  solid.  For  photons  of 
energy  hv  =  10  eV  one  has  wavelength  A  =  103  A  and  can  therefore  assume,  at  least 
for  a  wide  range  of  experiments,  that  the  wavelength  is  large  compared  to  atomic 
distances;  thus,  A  can  be  assumed  constant.  If  non-linear,  two-photon  processes  are 
neglected,  then  the  term  A  •  A  may  also  be  ignored,  and  (23)  reduces  to 


H'=—  (A-p).  (24) 

me 

The  electronic  states  in  a  solid  are  described  by  bands,  which  are  typically  plotted 
in  the  reduced  zone  scheme.  Within  this  scheme,  the  bands  outside  of  the  first 
Brillouin  zone  are  folded  back  into  the  first  zone  by  adding  the  appropriate  reciprocal 
lattice  vector  G.  Thus,  in  the  reduced  zone  scheme,  optical  transitions  between  initial 
i  and  final  /  states  are  “vertical”  such  that  k f  =  k;,  where  k  is  the  wave  vector 
of  the  electronic  state  in  the  crystal,  whereas  in  the  extended  zone  scheme,  k f  = 
kj  +  G.  Neglecting  the  momentum  of  the  photon,  the  internal  energy  distribution 
of  photoexcited  electrons  N-mt(E,  hu)  where  E  is  the  final  electron  kinetic  energy  is 
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determined  via  summation  over  the  initial  and  final  states.  Thus,  if  the  the  square  of 
the  transition  matrix  element  of  the  perturbing  Hamiltonian  operator  given  in  (24) 
is  re-written  as  |M/j(k*,  k/)|2,  then  Nint(E,hv)  is  determined  by  [1] 

Nint(E,hu)  oc^  |M/i(k  i,kf)\2  5(Ef{kf)  -  Ei{  kt)  -  hv) 

U 

xS(E-[Ef(kf) -(/>]),  (25) 

where  Ef(kf)  and  Ej(kj)  denote  the  energies  of  the  final  bands  | /,  kj)  and  the  initial 
bands  |i,kj),  respectively.  The  first  ^-function  imposes  energy  conservation  during 
excitation  and  the  additional  5- function  term  8(E  —  [Ef(kf)  —  0])  accounts  for  the 
fact  that  only  photoemitted  electrons  with  an  energy  above  the  vacuum  level  0  can 
be  experimentally  detected.  In  other  words,  the  kinetic  energy  measured  outside  the 
sample  equals  the  final  state  energy  inside  minus  the  work  function,  as  originally 
shown  in  (16). 

2.1.2  Transport  of  the  Electron  to  the  Surface. 

Following  optical  electronic  excitation,  scattering  mechanisms,  dominated  by 
electron-electron  interactions,  reduce  the  number  of  photoexcited  electrons  that  reach 
the  sample  surface  with  energy  Ef{kf).  Photoelectrons  that  avoid  these  interactions 
during  transport  to  the  surface  retain  their  momentum  and  energy  and  increase  the 
likelihood  of  escape.  If  the  scattering  frequency  1/r,  where  r  is  the  electron  lifetime, 
is  isotropic  and  dependent  on  E  only,  then  the  electron  inelastic  mean  free  path  A 
is  directly  proportional  to  the  probability  of  successful  transport  to  the  surface  such 
that 

A  {E,k)=rv,=  Q^,  (26) 

where  vg  is  the  group  velocity  in  the  final  state. 
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Berglund  and  Spicer  [4]  formulated  a  classical  treatment  of  electron  transport 
described  by  the  transmission  coefficient  d(E,  k )  as 

(27) 

where  a  is  the  optical  absorption  coefficient  of  the  photon  and  has  units  of  inverse 
length.  Qualitatively,  a-1  describes  the  penetration  depth  of  the  photon,  A  is  the 
average  distance  traveled  by  the  photoexcited  electron,  and  d(E,  k )  represents  the 
fraction  of  photoelectrons  created  within  one  mean  free  path  from  the  surface.  For 
the  energy  range  of  most  PES  experiments,  a-1  «  100  —  1000  A  and  A  «  10  —  20  A, 
which  implies  that  a  A  <§:  1  [1].  In  this  limit,  in  which  the  mean  free  path  of  the 
electron  is  smaller  than  the  penetration  depth  of  the  incident  photon,  inspection 
of  (27)  shows  that  d(E,  k )  — >  a\,  which  corresponds  to  a  range  of  approximately 
0.1  —  0.01  for  d(E,k).  Using  a  typical  density  of  states  (~  1023cm-3)  and  the 
dimensions  of  sample  surfaces  in  PES  experiments,  the  range  of  d(E,  k )  would  not 
be  expected  to  significantly  inhibit  electron  transport  to  the  surface.  Therefore,  it 
is  likely  that  a  sufficient  quantity  of  electrons  will  escape  from  the  surface  to  be 
captured  by  the  analyzer. 


2.1.3  Escape  of  the  Electron  into  Vacuum. 

In  a  free  electron  model,  escaping  photoelectrons  must  have  sufficient  kinetic 
energy,  normal  to  the  surface,  to  overcome  the  surface  potential  barrier 


E2Kl 

2m 


A  Evac  Eq 


(28) 


where  Kj^  is  the  perpendicular  component,  in  the  direction  of  the  surface  normal,  of 
the  photoelectron  wave  vector  K  inside  the  crystal,  and  Evac  and  E0  are  the  energy 
levels  of  the  vacuum  and  at  the  bottom  of  the  valence  band,  respectively.  It  follows 
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that  the  minimum  value  of  Kj^  required  to  escape  an  assumed  perfect  crystal  surface 

is 

Kl  min  =  (£vac  ~  ^0)V2  •  (29) 

It  is  a  key  assumption  of  the  three-step  model  that  electron  transmission  through 
the  surface-vacuum  interface  leaves  the  parallel  component  of  the  wave  vector  con¬ 
served  such  that 

P||/^  —  K||  =  k||  +  G|| ,  (30) 

where  p||  is  the  photoelectron  momentum  in  vacuum,  K||  is  the  parallel  component 
of  the  wave  vector  K  of  the  photoexcited  electron  inside  the  crystal,  k||  is  the  parallel 
wave  vector  component  of  the  electron  initial  state,  and  G||  is  the  parallel  component 
of  the  reciprocal  lattice.  Equation  (30)  for  the  free  electron  model  represents  Snell’s 
law,  such  that 

/  2m  \  ’^2  r2m  1 x/2 

k/n  =  sind  (—Ekinj  =  sin  d'  —(Ef-E0)  ,  (31) 

where  d  and  i)'  are  the  angles,  with  respect  to  the  surface  normal,  outside  and  inside 
the  crystal,  respectively. 

The  photoelectron  kinetic  energy  outside  the  crystal  could  be  determined  via 
knowledge  of  the  parallel  and  perpendicular  photoelectron  momentum  components 
in  vacuum,  py  and  pj_,  respectively,  by 

£ki„  =  ^  [(p,| /hf  +  (p  Jhf]  =  Ef{  k)  -  Evac  ,  (32) 

but  pj_  cannot  be  determined  explicitly  from  the  above  equations.  Alternatively,  in 
a  free  electron  solid,  the  crystal  potential  depth  Evac  —  E0  could  be  determined,  from 
which  Kj_,  the  full  wave  vector  of  the  crystal  k,  and  mapping  of  the  band  structure 
E(k)  are  determined.  However,  the  free  electron  model  is  not  accurate  enough  to 
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permit  full  determination  of  Ki.  Specifically,  the  model  fails  owing  to  the  fact  that 
the  wave  function  'h  /-(k)  of  a  state  at  Ef(  k)  in  a  real  solid  is  not  a  single  plane  wave, 
but  a  Bloch  wave  containing  plane  wave  contributions  from  various  reciprocal  lattice 
vectors.  Furthermore,  the  individual  components  can  be  attributed  to  photoelectrons 
emerging  from  the  crystal  in  many  possible  directions,  dictated  by  (31)  and  (32). 

Plane  wave  components  of  an  individual  Bloch  wave  of  final  energy  state  E  with 
identical  values  of  ky  +  G||  leave  the  crystal  in  the  same  direction  and  are  added 
coherently.  For  such  a  k||  +  G||,  the  total  transmission  factor  \T(Ef,  K| | ) |2  of  final 
energy  state  Ef  is  expressed  as  the  sum  of  reduced  transmission  factors  \t(Ef,K.\\)\ 
for  each  plane  wave  w/(G,k)  as 


2 


lr(-E/.kM)|2 


|i(£/,K„)|2 


£  «/(G,k) 

(fc+Gx)>0 


(33) 


where  the  summation  includes  only  those  components  propagating  toward  the  surface 
(k  +  GjJ  >  0.  An  expression  for  pj _/h  can  be  found  from  (32) 


pi 


p 


=  Ef(k)  -  Evac  -  -dh  , 


2  m 


2m 


(34) 


and  inserting  p||  =  k||  +  G  from  (30)  yields  a  “classical”  estimate  for  the  reduced 
transmission  factor  tf(Ef,  K||)  as 


!*(£/.  K||)|2 


1  if  Ef( k)  -  >  R2(kN  +  G,|)2/2m 

0  if  Ef( k)  -  £voc  <  «2(k||  +  G|,)2/2m 


(35) 


Although  this  is  an  overestimate  as  electrons  experience  inelastic  scattering  processes 
due  to  surface  plasmons,  meaning  tf(Ef,  K||)  <  1,  analysis  of  experimental  data  has 
shown  that  the  three-step  model  is  a  useful  and  accurate  approximation.  Assembling 
the  pieces  gives  a  final  expression  for  the  angle  resolved  energy  spectrum  of  primary 
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photoelectrons 


JV(£,KM,Hpri  |M/,(k,.k/)|2<i(£/,k/)|T(£/,Kll)|2 

f,i 

x  S(Ef( kf)  -  Efc)  -  hu)6(E  -  [Ef( kf)  -  0]) 
x  <5(ki  +  G  -  K)5(K||  -  pM  (tf,  ip) /h) .  (36) 

The  last  5  function  ensures  that  the  component  of  momentum  parallel  to  the  crystal 
surface  is  conserved  inside  and  outside  the  crystal. 

2.2  Resonant  Photoemission  Spectroscopy 

Resonant  photoemission  spectroscopy  (RPES)  is  an  interesting  application  of  PES 
that  provides  information  about  occupied  electronic  states  in  a  system.  However,  the 
key  value  of  resonant  photoemission  to  the  study  of  rare  earth  doped  GaN  thin  films  is 
to  probe  which  valence  bands  of  the  semiconductors  have  strong  4/  and/or  rare  earth 
weight  [7-9].  Resonance  occurs  when  one  stimulates  photoemission  using  photons 
of  energy  hu  near  an  absorption  threshold  of  a  core  level  electronic  state  of  binding 
energy  huy  For  the  rare  earth  dopants  studied  in  this  research,  resonance  results 
from  a  signal  overlap  between  direct  emission  of  photoelectrons  from  the  4/  state 

4d104fN  +  hu  — >  4d104/Ar_1  +  e~  ,  (37) 

and  Auger-like  electrons  emitted  in  a  super  Coster-Kronig  process  [10,11] 

4dw4fN  +  hu  — >  [4d94fN+1]*  — »  id1HfN~1  +  e~  ,  (38) 

where  [  ]*  denotes  an  excited  state.  The  final  states  for  both  the  direct  and  recom¬ 
bination  processes  are  identical. 
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Figure  8  depicts  a  general  schematic  of  the  process  for  a  hypothetical  solid,  using 
a  photon  energy  approaching  the  4 d  core  absorption  threshold,  and  is  explained  as 
follows.  Direct  photoemission  occurs  from  the  4/  state,  as  shown  in  Figure  8(a) 

•e' 


(b)  Valence  band  Coster-Kronig  process. 


Figure  8.  Schematic  molecular-orbital  diagram  of  the  resonant  photoemission  process 
for  a  hypothetical  solid. 


and  the  emission  process  (37).  Additionally,  electrons  from  the  4 d  core  state  absorb 
photons  and  are  excited  to  a  bound,  unoccupied  4/  valence  band  state.  The  4 d  hole  is 
subsequently  occupied  by  a  4/  electron  (de-excitation),  and,  if  the  energy  difference 
between  the  4/  and  4 d  states  is  sufficient,  4/  electron  emission  occurs  as  depicted  in 
Figure  8(b)  and  process  (38).  The  result  is  an  enhanced  photoemission  signal  due  to 
the  overlap  between  direct  emission  of  the  valence  band  electrons  and  the  Auger-like 
electron  emission  following  excitation  and  decay  between  the  core  electronic  level 
and  the  bound  valence  band  state;  both  processes  result  in  an  identical  final  state  of 
the  system. 

The  transitions  of  Figure  8(b)  represent  a  special  case  of  the  Auger  process, 
referred  to  as  a  Coster-Kronig  transition  [10, 11],  in  which  the  vacancy  is  filled  by  a 
higher  electron  subshell  of  the  same  electron  shell.  When  the  emitted  electron  also 
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belongs  to  the  same  shell,  as  shown  in  Figure  8(b),  the  process  is  called  a  super 
Coster-Kronig  transition.  Thus,  the  schematic  shows  both  resonant  photoemission 
and  a  super  Coster-Kronig  transition. 

The  intensity  of  the  valence  band  features  for  which  the  direct  and  super  Coster- 
Kronig  processes  overlap  is  described  as  a  function  of  photon  energy  hu  by  the  Fano 
lineshape  [11,12]  as 


N(hu)  2* 


(e  +  g )2 

e2  +  1 


hu  —  huj 

f 


(39) 


where  huj  is  the  photon  energy  equal  to  the  specified  core  level  absorption  threshold 
binding  energy,  q  is  the  line  profile  index  fitting  parameter  for  the  core  level  and  is 
related  to  the  Coulomb  and  dipole  matrix  elements  [12],  and  T  =  A(huj) /2  where 
A  (huj)  is  the  width  (FWHM)  of  the  core  level.  Thus,  fitting  (39)  with  experimental 
data  is  a  common  analytical  technique. 

Figure  9  displays  the  functional  form  of  (39)  and  highlights  the  importance  of 
the  fitting  parameter  q.  Inspection  of  (39)  shows  that  as  hu  — >  huj ,  N(E ,  hu)  — >  q2. 
It  follows  in  Figure  9  that  the  difference  between  N(hu  «  140)  and  N{hu  «  140) 
is  approximately  q 2  =  4,  such  that  the  fitting  parameter  q  provides  an  estimate 
of  the  relative  photoemission  signal  between  values  of  hu  nearby  and  far  from  huy 
Therefore,  when  analyzing  the  data,  special  attention  should  be  given  to  the  fitting 
parameter  q,  as  an  indicator  of  the  resonant  photoemission  strength. 

Resonant  photoemission  can  be  used  to  confirm  or  deny  the  presence  of  intended 
or  unintended  impurity  ions  in  a  solid.  When  coupled  with  other  PES  measure¬ 
ments,  such  as  angle-resolved  ultraviolet  photoemission  spectroscopy  (ARUPS)  or 
temperature-dependent  X-ray  photoemission  spectroscopy,  inferences  can  be  made 
about  the  surface  concentration  of  impurities,  as  well  as  the  location  of  impurities  in 
a  crystal  structure  as  either  interstitial  or  substitutional  ions,  as  will  be  discussed  in 
the  next  section. 
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Figure  9.  Plot  of  Fano  function  from  (39)  using  huj  —  140  eV,  A hvj  —  3  eV,  and 
<7  =  2. 

2.3  Temperature-Dependent  X-ray  Photoemission  Spectroscopy 

PES  experiments  are  often  performed  at  room  temperature,  which  is  sufficient  to 
vibrationally  excite  atoms  or  molecules  in  a  crystal  structure.  These  vibrations  can 
couple  to  atomic  electrons  via  electron-phonon  interactions  and  influence  the  PES 
data.  Specifically,  if  PES  is  considered,  phenomenologically,  as  an  electron  scattering 
process  by  the  ions  of  a  crystal,  then  detected  electrons  must  reflect  the  effect  of  the 
lattice  vibrations. 

A  free  electron  model  treatment  of  valence  electrons  creates  a  minor  pedagogical 
problem  as  free  electrons,  being  only  slightly  perturbed  by  the  crystal  potential, 
are  hardly  influenced  by  lattice  vibrations.  However,  if  a  tight-binding  approach  is 
considered,  such  that  the  electrons  are  rigidly  attached  to  the  vibrating  ion  cores, 
the  effect  of  lattice  vibrations  is  introduced  and  the  pedagogical  problem  becomes 
somewhat  transparent.  Under  this  approximation,  the  transition  matrix  element  of 
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the  perturbing  Hamiltonian  from  the  Golden  Rule  (18)  is  calculated  [1]  as 

(/|H,|i)  oc  crpe  S ,  (40) 

where  ape  is  the  photoabsorption  cross  section  and  S  is  the  scattering  factor 

S  =  y^exp[iAk  •  (RM  —  R„)];  Ak  =  k/  — k/.  (41) 

/J,,U 

Indices  ji  and  v  span  the  atomic  sites,  not  the  electron  initial  and  final  states,  and 
represents  the  coordinate  of  the  atom.  Thus,  lattice  vibrations  are  accounted 
for  as 

R =  RM,„  +  u p,v(T) ,  (42) 

where  uM)i,  is  the  displacement  of  the  respective  ion.  The  temperature-dependent 
photocurrent  is  then  determined  as 


N(E,  hu,  T )  oc|crpe|2  <  exp[— A/c2Mq(T)]  5(Ak  —  G)  +  N(  1  —  exp[— Ak2ul(T)]) 

{  G 

direct  PE  current  indirect  PE  current 

(43) 

The  direct  current  term  represents  the  temperature-dependent  attenuated  photoe¬ 
mission  signal,  N  is  the  number  of  atoms  that  contribute  to  the  signal,  and  the 
indirect  current  results  from  phonon  assisted  momentum  transitions  that  occur  at  a 
finite  temperature. 

Qualitatively,  as  the  temperature  of  a  solid  is  increased,  lattice  vibrations  in¬ 
crease,  which  reduce  an  electron’s  ability  to  traverse  through  the  material,  escape, 
and  get  detected.  Quantitatively,  the  intensity  of  an  emitted  or  scattered  electron 
beam  decays  exponentially,  clue  to  increased  thermal  vibrations  from  increased  tem- 
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perature  as  [1, 13] 


I  =  I0e~2W  ,  (44) 

where  Jo  is  the  reference  intensity,  which  is  typically  the  lowest  temperature  spec¬ 
trum.  The  Debye- Waller  factor  W  is  expressed  as 

2W  =  Ak2(u0)2  ,  (45) 


where  A k  is  the  wave  vector  transfer  and  (uq)2  is  the  mean  square  displacement  of 
the  atoms.  Within  the  Debye  model  of  thermal  vibrations,  in  the  case  of  isotropic 
vibrations,  W  is  described  as 


3h2(Ak)2T 
m  ks  Q2d 


(46) 


In  (46),  h(Ak)  is  the  electron  momentum  transfer  due  to  phonon  coupling,  m  is 
the  mass  of  the  scattering  ion,  ks  is  the  Boltzmann  constant,  and  ©d  is  the  Debye 
temperature.  A k  is  calculated  from  the  experimentally  measured  electron  kinetic 
energy  Ekin  as 


(A  k)2 


2  me  -Ekin 
K2 


(47) 


where  me  is  the  electron  mass.  From  equations  (44)  and  (46), 


In 


-2  W 

3H2(Ak)2T  T 

m  kB  h‘i,  Q2d  ' 


(48) 


Therefore,  when  plotted,  the  slope  \S\  of  the  logarithm  of  relative  intensity  I/Iq 
as  a  function  of  temperature  is  inversely  proportional  to  the  square  of  the  Debye 
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temperature,  as 


3/t2(A  kf 
m  ks  Q2d 

As  mentioned,  resonant  photoemission  can  be  combined  with  temperature-dependent 
XPS  to  draw  conclusions  about  the  presence  and  location  of  impurity  ions  in  a  crys¬ 
tal  structure.  As  an  example,  consider  a  ternary  compound  semiconductor  AxBi_iTC 
where  ion  A  is  intended  to  be  substitutional  for  ion  B  in  x%  of  B  sites  within  the 
lattice  structure  BC.  Resonant  photoemission  can  effectively  confirm  the  presence  of 
ion  A  by  scanning  through  photon  energies  near  an  absorption  threshold  of  a  core 
level  electronic  state  of  A.  By  graphing  the  intensities  of  the  resulting  spectral  fea¬ 
tures  and  analyzing  the  energetically  allowed  electronic  transitions,  one  can  make 
inferences  about  the  occupied  states  of  the  ternary  compound.  If  this  analysis  shows 
an  enhancement  of  spectral  features  located  at  electron  binding  energies  that  are 
consistent  with  ion  A,  then  one  may  conclude  that  the  ion  was  successfully  incor¬ 
porated  into  the  compound,  and  inferences  about  band  hybridiztion  can  be  made. 
However,  the  precise  location  of  ion  A  within  within  the  crystal  structure  would  be 
undetermined.  Temperature-dependent  XPS  analysis  and  calculation  of  the  Debye 
temperatures  for  ions  A  and  B  permits  inferences  about  the  location  of  ion  A.  Specif¬ 
ically,  if  ©d  for  A  and  B  are  equal,  within  experimental  uncertainty,  then  one  may 
conclude  that  the  electron-phonon  interactions  for  both  ions  are  identical,  and  one 
may  classify  ion  A  as  substitutional  for  ion  B  in  the  structure.  On  the  other  hand, 
if  the  ions’  Debye  temperatures  are  not  equal,  then  one  may  conclude  that  ion  A  is 
vibrating  from  an  interstitial  location  within  the  crystal  structure. 

2.4  The  Schottky  Barrier 

The  metal-semiconductor  interface  spans  the  subject  of  surface  science  and  semi¬ 
conductor  theory,  as  techniques  from  both  fields  can  effectively  measure  interface 
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properties.  One  such  property,  the  Schottky  barrier,  is  a  potential  barrier  that  exists 
at  the  interface  between  a  metal  and  semiconductor  and  is  related,  roughly,  to  the 
difference  in  work  functions  between  materials.  This  barrier  is  responsible  for  con¬ 
trolling  current  conduction  and  capacitance  behavior  and  therefore  warrants  careful 
consideration. 

As  an  example,  the  interface  electronic  properties  of  many  metal/GaN  systems 
have  been  directly  measured  and  studied  using  photoemission  spectroscopy  and  other 
techniques.  Specifically,  Cs,  Mg,  Al,  Ti,  Ni,  Au,  Pd,  and  Pt  on  GaN  surfaces  have 
been  reported  [14-26]  and,  more  recently,  Xiao  et  al.  presented  an  original  study  of 
rare-earth  metals  (Gd)  on  GaN  [27]. 

Alternatively,  electronic  measurement  of  interface  properties  via  current-voltage 
(I-V)  and  capacitance-voltage  (C-V)  characteristics  also  facilitate  Schottky  barrier 
height  calculation.  However,  the  benefit  of  direct  surface  barrier  height  measure¬ 
ment  is  avoidance  of  complications  associated  with  other  experimental  techniques. 
For  example,  defects  at  the  metal-semiconductor  interface  can  lead  to  complications 
in  I-V  measurements  of  the  surface  barrier  height  [28,29].  Similarly,  C-V  determina¬ 
tion  of  the  barrier  height  can  be  complicated  by  these  defects,  which  alter  the  space 
charge  region  and  affect  the  measured  flat  band  voltage  [28,29].  Accordingly,  the 
sections  that  follow  will  discuss  the  theoretical  models,  underlying  assumptions,  and 
experimental  techniques  associated  with  direct  spectroscopic  measurement  and  indi¬ 
rect  electronic  measurement  of  the  Schottky  barrier  height  of  a  metal-semiconductor 
interface.  The  final  section  will  discuss  the  Schottky  barrier  relevance  for  radiation 
detector  designers. 

2.4.1  Schottky-Mott  Model. 

The  Schottky-Mott  model,  which  is  the  most  simplistic  physical  model  and  as¬ 
sumes  intimate  and  abrupt  contact  as  shown  in  Figure  10,  predicts  a  barrier  height 
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at  the  metal  to  n-type  semiconductor  interface  of 


*&B,n  0ra  Xs  5 


(50) 


where  (pm  and  Xs  are  the  metal  work  function  and  semiconductor  electron  affinity, 
respectively.  Closer  inspection  of  Figure  10(b)  shows  that  an  alternate  expression 
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(a)  Energy  band  diagram  of  an  isolated  metal  (b)  Energy  band  diagram  of  a  metal- 
adjacent  to  an  n-type  semiconductor  under  semiconductor  contact  in  thermal  equilibrium, 
thermal  non-equilibrium  condition. 


Figure  10.  Ideal  Schottky  barrier  formation. 


for  <&B>n  is 

$B,n  =  Eg  —  (Ep  —  Ey)  ,  (51) 

where  Eg,  Ep,  and  Ey  are  the  band  gap,  surface  Fermi,  and  surface  valence  band 
maximum  energies,  respectively,  of  the  semiconductor.  Photoemission  spectroscopy 
permits  direct  measurement  of  these  energy  levels,  the  surface  valence  band  bending, 
and,  therefore,  the  Schottky  barrier  height  as  metal  overlayers  are  deposited  on  a 
semiconductor  surface. 

In  practice,  the  simple  theoretical  expression  for  barrier  height  as  given  by  equa¬ 
tion  (50)  is  never  realized  experimentally.  Specifically,  an  unavoidable  interface  layer 
between  the  metal  and  semiconductor  surfaces  as  well  as  the  presence  of  semiconduc¬ 
tor  surface  states  alter  the  measured  Schottky  barrier  from  the  ideal  value  calculated 
by  the  Schottky-Mott  model.  Furthermore,  on  an  atomic  scale,  the  transition  is 
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not  abrupt  but  gradual,  with  local  fields  varying  from  atom  to  atom.  One  model 
used  to  explain  practical  interfaces  is  the  metal  induced  gap  states  (MIGS)  model 
which  claims  that  defects  due  to  metal  deposition  and  the  presence  of  semiconductor 
surface  states  result  in  interface  traps  that  trap  charge  and  change  the  ideal  barrier 
height. 

2.4.2  Metal  Induced  Gap  States  Model. 

When  seeking  a  more  complete  theoretical  model  to  compare  with  the  measured 
results,  both  chemical  bonding  and  physical  concepts  must  be  considered.  Exten¬ 
sive  measurements  using  photoemission  spectroscopy  indicate  that  for  most  covalent 
semiconductors,  and  specifically  III-V  compounds,  the  Schottky  barrier  formation  is 
due  mainly  to  defects  generated  near  the  metal-semiconductor  interface  by  deposi¬ 
tion  of  the  metal  [30].  On  several  compound  semiconductors  such  as  GaAs,  GaSb, 
and  InP,  the  surface  Fermi-level  positions  obtained  from  various  metals  are  pinned 
at  an  energy  level  quite  independent  of  the  metal  [31].  Thus,  the  Fermi  level  pinning 
explains  the  fact  that  the  barrier  height  in  covalent  semiconductors  is  essentially 
independent  of  the  metal  work  function. 

For  ionic  semiconductors,  the  barrier  height  generally  depends  strongly  on  the 
metal  and  a  correlation  has  been  found  between  interface  behavior  and  the  electroneg¬ 
ativity  difference  between  the  metal  and  semiconductor  materials  [32,33].  Specifi¬ 
cally,  the  interface  bonds  at  metal-semiconductor  contacts  are  partly  ionic,  as  are 
all  heteropolar  bonds  [33].  Generalizing  Pauling’s  concept  [34]  that  the  ionicity  of 
covalent  single  bonds  is  described  by  the  electronegativity  difference  between  the 
atoms,  the  charge  transfer  at  a  metal-semiconductor  interface  is  characterized  by  the 
electronegativity  difference  between  the  materials. 

Using  simple  physical  arguments  as  a  basis,  Heine  [35]  pointed  out  that  wave 
functions  of  the  metal  electrons  tail  into  the  semiconductor  in  the  energy  range  where 
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the  metal  conduction  band  overlaps  the  semiconductor  band  gap.  He  concluded  that 
this  continuum  of  metal  induced  gap  states  represents  the  fundamental  mechanism 
that  determines  Schottky  barrier  height.  The  schematic  energy  diagram  in  Figure  11 
shows  a  practical  metal-semiconductor  interface  with  an  interface  layer  of  width  5 
between  the  surfaces. 


I 


Figure  11.  Schematic  energy  diagram  of  a  practical  metal-semiconductor  interface 
with  a  continuum  of  interface  states. 


The  MIGS  arise  from  nearby  bands  and  are  predominantly  donor-like  (charge 
neutral  when  occupied,  positively  charged  when  unoccupied),  closer  to  the  valence 
band  maximum  and  mostly  acceptor-like  (charge  neutral  when  unoccupied,  nega¬ 
tively  charged  when  occupied),  nearer  to  the  conduction  band  minimum.  The  energy 
level  (fibp  above  the  valence  band  maximum  at  which  the  dominant  behavior  changes 
is  called  the  branch  point.  Consequently,  when  is  above,  coincides  with,  or  is 
below  the  Fermi  level,  the  net  charge  in  these  interface  states  is  positive,  vanishes,  or 
negative,  respectively.  Lastly,  from  the  schematic  diagram,  the  zero-charge  transfer 
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barrier  height  4>zc  is  calculated  as 


4>zc  =  Eg  -  ( 4>bp  ~  Ev) .  (52) 

When  the  density  of  interface  traps  Dlt  in  the  region  5  is  large  ( Dlt  — >  oo),  the 
Schottky  barrier  height  approaches  the  value  given  by  equation  (52).  Inspection  of 
the  equation  indicates  that  barrier  height  is  independent  of  the  metal  work  function 
and  is  determined  entirely  by  the  interface  properties  of  the  semiconductor.  In  this 
case,  the  Fermi  level  at  the  interface  is  said  to  be  pinned,  by  the  surface  states  at  a 
value  4>bp  above  the  valence  band  maximum  [36].  In  the  other  limiting  case  (Dlt  — >  0), 
the  Schottky  barrier  height  approaches  the  value  given  by  the  Schottky-Mott  model 
in  equation  (50). 

For  the  case  0  <  Du  <  oo,  the  Schottky  barrier  height  is  linearly  corrected  from 
(52)  by  factors  that  are  dependent  upon  material  properties  of  both  the  metal  and 
the  semiconductor.  Specifically,  in  combining  the  physical  argument  of  MIGS  and 
the  chemical  concept  of  partly  ionic  interface  bonds,  Monch  [33]  predicted  Schottky 
barrier  heights  &B,n  to  be  larger  or  smaller  than  the  zero-charge  transfer  barrier 
height  (j>zc  if  the  electronegativity  difference  between  the  metal  and  semiconductor 
(Xm  —  Xs)  was  positive  or  negative,  respectively,  as  given  by 

®B,n  =  4>zc  +  Sx(Xm  —  Xs ) .  (53) 

The  relationship  in  (53)  simply  states  that  the  Schottky  barrier  varies  linearly 
as  a  function  of  the  electronegativity  difference  where  Sx  is  the  slope  parameter  as 
determined  by 

4^-l»0.1(eoo-l)2.  (54) 

Here  Ax  is  calculated  from  the  MIGS  model  as  either  0.86  eV/Miedema-unit  or 
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1.79  eV/Pauling-unit  when  using  the  Miedema  or  Pauling  values  for  electronegativ¬ 
ity,  and  Coo  is  the  optical  dielectric  constant  of  the  bulk  valence  electrons  of  the 
semiconductor  [33] .  Thus,  a  theoretical  value  for  <f> s,n  can  be  determined  from  equa¬ 
tion  (53)  using  material  properties  and  by  calculating  the  branch  point  energy  above 
the  valence  band  maximum  in  order  to  solve  (52). 

At  first  glance,  the  relationship  between  barrier  height,  the  slope  parameter  Sx, 
and  the  semiconductor  might  appear  unconvincing.  Full  theoretical  development 
of  (54)  requires  a  rather  rigorous  analysis  of  the  charge  transfer  dynamics  at  the 
metal-semiconductor  interface,  which  can  be  found  elsewhere  [36,37].  While  this 
analysis  elucidates  the  connection  between  the  physical  relationship  and  the  govern¬ 
ing  parameters,  the  association  between  them  may  be  better  understood  graphically, 
as  is  done  within  the  reported  literature.  As  an  example,  Figure  12  displays  the 
relationship  between  the  electronegativity  difference  Xm  —  Xs  at  the  hypothetical 
metal  to  semiconductor  interface  between  various  low  and  high  work  function  metals 
and  n-type  silicon.  Just  as  with  silicon  in  the  example,  these  figures  and  theoretical 
relationships  are  determined  by  experimentally  measuring  barrier  heights  between 
various  metals  for  a  given  semiconductor  surface.  Careful  attention  must  be  given 
to  semiconductor  growth,  sample  and  surface  preparation,  metal  deposition  method, 
and  barrier  height  measurement,  in  order  to  ensure  consistent  experimental  con¬ 
ditions;  thus,  generating  a  theoretical  MIGS  relationship  via  data  extrapolation  is 
a  significant  undertaking.  However,  these  theoretical  relationships  are  valuable  in 
determining  whether  or  not  MIGS  play  a  role  at  a  particular  metal-semiconductor 
interface. 

Whether  or  not  metal  induced  gap  states  play  a  role  at  the  metal  to  GaN  thin 
film  interface  is  a  subject  of  debate.  Kurtin,  McGill,  and  Mead  [38]  suggested  that 
the  Schottky  barrier  height  on  GaN  should  depend  directly  on  the  work  function  or 
electronegativity  difference  between  the  metal  electrode  and  GaN.  Foresi  and  Mous- 
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Figure  12.  Barrier  heights  of  various  metal  to  n-Si  contacts  as  a  function  of  electroneg¬ 
ativity  difference.  The  theoretical  MIGS  line  (solid)  is  drawn  using  a  least  squares  fit 
of  the  data  with  <I?ZC  =  0.792  and  Sx  —  0.113.  Adapted  from  [1]. 

takas  [39]  observed  this  direct  correlation  experimentally,  while  Gno  et  al.  [40]  and 
Mori  et  al.  [41]  observed  only  a  weak  dependence  of  the  Schottky  barrier  height  on  the 
metal  work  function  for  n-type  GaN  and  p-type  GaN,  respectively.  The  1998  review 
of  metal-GaN  contact  technology  by  Liu  and  Lau  [42]  reported  that,  for  a  variety  of 
contact  metals  with  both  low  and  high  work  functions,  Schottky  barrier  heights  at 
the  metal-GaN  interface  varied  with  metal  work  function,  within  the  experimental 
scatter.  Subsequent  work  by  Rickert  et  al.  [28]  supported  a  modified  Schottky-Mott 
model  at  the  metal-GaN  interface  for  Ni,  Pd,  and  Al,  yet  more  “complex”  behavior 
when  Au,  Ti,  and  Pt  were  used  as  the  contact  metals.  Additional  experiments  by 
Barinov  et  al.  [22,23]  reported  Schottky  barrier  heights  at  the  Au-GaN  interface  that 
exceeded  both  work  function  difference  (Schottky-Mott)  and  electronegativity  differ¬ 
ence  (metal  induced  gap  states)  models.  Thus,  a  point  worth  emphasizing  is  that 
regardless  of  the  particular  metal-GaN  interface  studied,  experimentally  measured 
barrier  heights  vary  considerably. 
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The  final  section  of  this  chapter  discusses  a  method  of  Schottky  barrier  height 
calculation,  which  was  not  conducted  in  this  research,  but  is  common  in  the  literature. 
This  section  is  presented  for  several  reasons.  First,  the  Schottky  barrier  heights 
at  the  metal  to  rare  earth  doped  GaN  interfaces,  determined  using  photoemission 
spectroscopy, will  be  compared  with  previously  reported  results  that  were  derived 
from  both  spectroscopic  data  and  device  I-V  and  C-V  characteristics.  For  this  reason, 
a  cursory  understanding  of  current  transport  in  Schottky  barrier  devices  is  necessary. 
Additionally,  the  application  of  the  results  to  devices,  and  particularly  to  radiation 
detectors,  requires  at  least  a  brief  discussion.  Lastly,  the  topic  and  its  application 
represent  areas  for  potential  future  research  in  the  field  of  rare  earth  doped  GaN  thin 
films. 

2.4.3  Current  Transport  in  Schottky  Barrier  Devices. 

Current  transport  across  a  Schottky  barrier  results  from  thermionic  emission  of 
majority  carriers  from  the  semiconductor  over  the  electrostatic  potential  barrier  into 
the  metal,  or  from  carrier  tunneling  through  the  barrier,  and  measurement  of  device 
I-V  properties  allows  calculation  of  this  barrier. 

In  thermal  equilibrium,  and  considering  an  n-type  semiconductor,  the  emission  of 
electrons  from  the  semiconductor  into  the  metal  is  balanced  by  the  flow  of  electrons 
from  the  metal  into  the  semiconductor.  The  concentration  of  electrons  nth  that 
possess  sufficient  energy  to  traverse  the  barrier  <f> Bn  is  [43] 

nth  =  Nc  exp  ,  (55) 

L  kBT  J 

where  Nq  is  the  effective  density  of  states  in  the  conduction  band,  kB  is  Boltzmann’s 
constant,  and  T  is  the  temperature. 

At  thermal  equilibrium,  the  current  density  conditions  described  above  are  ex- 
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pressed  as 


Jm^s  |  |  J s^m  I  OC  nth 

Jm— »s|  I1/*— >m\  G\Nc  GXp 


&B,n 

k^T 


(56) 


where  C\  is  a  proportionality  constant,  and  Jm^s  and  Js^m  are  the  current  densities 
from  the  metal  to  the  semiconductor  and  from  the  semiconductor  to  the  metal, 
respectively. 

An  applied  forward  bias  Vp  at  the  contact  reduces  the  electrostatic  potential 
difference  across  the  barrier  and  increases  the  electron  density  across  the  interface 
to 


nth  =  Nc  exp 


(®B,n  —  QVf) 

w 


(57) 


The  current  density  on  the  semiconductor  side  of  the  interface  Js^m  is  altered  by  the 
same  amount,  but  the  electron  current  density  remains  unchanged  from  the  metal 
to  the  semiconductor  because  the  barrier  remains  at  its  equilibrium  value  from  the 
viewpoint  of  electrons  on  the  metal  side  of  the  interface.  Therefore,  under  forward 
bias,  the  net  current  density  is 


J  Js—>m  Jm^rs 
=  C\Nc  exp 


( &B,n  —  QVf) 
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(58) 


The  same  argument  for  the  reverse  bias  condition  produces  an  identical  expression  to 
(58),  except  that  Vp  is  replaced  by  —Vr.  It  should  be  noted  that,  in  addition  to  the 
electron  (majority  carrier)  current,  a  hole  (minority  carrier)  current  exists  from  the 
metal  to  the  semiconductor  due  to  the  minority  carrier  injection.  However,  typical 
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minority  carrier  currents  are  orders  of  magnitude  smaller  than  those  of  minority 
carriers,  which  characterizes  Schottky  diodes  as  single  carrier  conduction,  or  unipolar, 
devices  [43]. 

The  coefficient  C\Nc  in  equation  (58)  can  be  shown  to  be  equal  to  A*T2  where 
A*  is  the  effective  Richardson  constant,  in  units  of  A/K2-cm2,  and  depends  directly 
on  the  electron  effective  mass  in  the  semiconductor  material.  Therefore,  the  current- 
voltage  characteristic  of  a  metal-semiconductor  contact  is 


Js  (eqV/kBT  -  1)  , 

(59) 

y[*  ^b,u/^bT 

(60) 

where  Js  is  the  saturation  current.  In  (59)  V  can  be  either  Vp  or  —  Vr,  depending 
on  bias  conditions. 

Experimental  calculation  of  Schottky  barrier  heights  via  equations  (59)  and  (60) 
is  done  as  follows.  Current  density  versus  applied  forward  voltage  is  plotted  and 
extrapolated  to  V  —  0,  yielding  a  value  for  Js,  from  which  the  barrier  height  is 
obtained  via  (60)  and  knowledge  of  the  effective  Richardson  constant.  This  device 
parameter  has  critical  implications  with  regard  to  radiation  detector  devices. 

Schottky  contacts  serve  an  important  function  as  rectifying  contacts  in  general 
semiconductor  device  design,  but  an  often  critical  function  for  semiconductor-based 
radiation  detectors.  Despite  efforts  to  limit  current  flow  due  to  leakage,  and  max¬ 
imize  the  charge  collection  volume,  the  generated  signal  due  to  radiation-induced 
currents  can  prove  difficult  to  detect  above  noise  due  to  leakage  currents.  In  such 
situations,  Schottky  contacts  are  desirable  to  reduce  the  noise  associated  with  leakage 
currents  and  maximize  the  signal  of  the  radiation-induced  currents.  Quantitatively, 
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the  contact  resistance  pc  is  determined  by  the  surface  barrier  height  [28]  as 


Pc 


qA*T 


exp 


$B,n 

kAT 


(61) 


It  is  evident  from  (61)  that  small  reductions  in  the  surface  barrier  height  will  have 
a  large  effect  in  reducing  the  contact  resistance  at  the  Schottky  contact,  which  ad¬ 
versely  affects  the  signal-to-noise  ratio  and  performance  of  a  diode  radiation  detector. 
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III.  Experimental  Procedures 


This  chapter  discusses  general  experimental  procedures  used  during  research.  Ad¬ 
ditional  details  pertaining  to  the  individual  experiments  can  be  found  in  Chapters  IV, 
V,  and  VI.  The  ultrahigh  vacuum  conditions  needed  for  photoemission  experiments 
required  that  all  experiments  occur  at  either  the  Center  for  Advanced  Microstruc¬ 
tures  and  Devices  (CAMD)  at  Louisiana  State  University  in  Baton  Rouge,  LA  or  at 
the  Nebraska  Center  for  Materials  and  Nanoscience  at  the  LIniversity  of  Nebraska 
(Lincoln)  campus.  Altogether,  13  experiments  (8  resonant  photoemission,  4  Schot- 
tky  barrier  measurement,  1  surface  Debye  temperature)  were  conducted  during  14 
total  weeks  spanning  November  2009-December  2010. 

3.1  Sample  Growth  and  Characterization 

The  rare  earth  doped  gallium  nitride  thin  films  were  prepared  at  the  LIniversity  of 
Puerto  Rico  and  the  LIniversity  of  Nebraska  (Lincoln)  using  plasma  assisted  molec¬ 
ular  beam  epitaxy  [1].  Molecular  beam  epitaxy  is  a  relatively  slow  deposition  rate 
process  (<  1000  nm/hour),  shown  schematically  in  Figure  13,  and  takes  place  in 
high  or  ultra  high  vacuum.  Vacuum  conditions  limit  the  introduction  of  unwanted 
impurities  during  the  growth  process  and  maximize  the  mean  free  paths  of  the  indi¬ 
vidual  atoms  in  order  to  minimize  their  interactions  between  the  respective  effusion 
cells  and  the  substrate  wafer.  During  film  grown,  ultra-pure  solid-form  elements  are 
heated  in  separate  effusion  cells  to  their  appropriate  sublimation  temperatures.  In 
their  gaseous  forms  and  within  the  nitrogen-rich  environment  indicated  by  the  figure, 
the  gaseous  elements  condense  epitaxially  on  the  substrate  wafer.  During  operation, 
reflection  high  energy  electron  diffraction  (RHEED)  monitors  crystal  layer  growth, 
and  computer-controlled  shutters  in  front  of  each  effusion  cell  allow  precise  control 
of  layer  thickness.  Details  pertaining  to  the  specific  growth  parameters  for  GaN:RE 
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Figure  13.  Typical  molecular  beam  epitaxy  schematic. 

thin  film  deposition  will  be  discussed  in  Chapters  IV,  V,  and  VI. 

Another  technique  which  was  used  to  prepare  samples  is  pulsed  laser  deposition 
[2],  depicted  in  Figure  14.  During  the  growth  process,  a  laser  beam  strikes  the  surface 
of  the  source  (target)  material,  and  the  deposited  energy  from  the  laser  evaporates 
the  target’s  surface.  The  short  laser  pulses,  typically  less  than  50  ns  in  duration, 
cause  the  surface  temperature  of  the  target  to  rise  rapidly,  while  the  bottom  of  the 
target  remains  cooler.  The  nonequilibrium  heating  produces  a  flash  of  evaporants 
that  deposit  on  the  substrate,  producing  a  film  with  composition  identical  to  that  of 
the  target  surface.  This  ability  to  replicate  the  composition  of  the  source  in  the  film 
is  perhaps  the  greatest  benefit  of  PLD.  Whereas  thermal  evaporation  techniques, 
such  as  molecular  beam  epitaxy,  produce  vapor  compositions  dependent  upon  the 
vapor  pressures  of  the  elements  in  the  target  material,  the  laser-induced  expulsion  of 
material  from  pulsed  laser  deposition  yields  a  plume  of  material  with  stoichiometry 
identical  to  the  target.  Lastly,  although  the  schematic  in  Figure  14  shows  only  a 
single  target  material,  experimental  configurations  capable  of  positioning  multiple 
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Laser  beam 


targets  within  the  laser  beam  exist. 

Following  growth,  film  thickness  was  measured  with  a  surface  profilometer,  and 
X-ray  diffraction  (XRD)  analysis  verified  orientation,  crystal  structure,  and  film 
purity.  XRD  was  conducted  either  using  a  Siemens  D5000  or  Rigaku  9- 2D  X-ray 
diffractometer,  using  a  copper  target  X-ray  tube  that  produced  radiation  at  the  Cu 
Ka  wavelength  (A  =  1.5406  A).  The  sample  and  detector  were  rotated  with  respect 
to  the  incident  beam  angle  6  and  spectra  were  collected. 

Figure  15  shows  the  X-ray  diffraction  pattern  of  a  GaN:Yb  film  grown  on  Si(lll) 
substrates.  The  Bragg  formalism  predicts  that  diffraction  maxima  between  parallel 
planes  of  ions  spaced  a  distance  d  apart  occur  when  incident  light  of  wavelength  A: 

(1)  reflects  from  a  crystal  plane  where  the  angle  of  incidence  6  equals  the  angle  of 
reflection; 

(2)  results  in  constructive  interference  of  reflected  rays  from  successive  planes. 
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Figure  15.  The  X-ray  diffraction  (XRD)  pattern  of  a  GaN:Yb  film  grown  on  Si(lll) 
substrates.  The  films  show  c-axis  orientation  of  the  texture  growth  and  a  high  level 
of  crystallinity.  The  presence  of  any  secondary  phases  or  spurious  peaks  has  not  been 
observed  [3]. 


The  condition  for  a  sharp  peak  in  the  intensity  of  scattered  radiation  from  a  crystal 

is 

n\  =  2dhki  sin(0) ,  (62) 


where  dhki  is  the  distance  between  planes  of  Miller  indices  (hkl) 


dhkl  — 


'h2  k 2  l2' 
~^+b2+d2 


-1/2 


(63) 


and  a,  b,  and  c  are  the  lattice  dimensions.  The  diffraction  maxima  in  Figure  15  occur 
for  the  GaN:Yb(0002)  and  GaN:Yb(0004)  planes  at  26  =  34.66°  and  29  =  73.14°, 
respectively.  From  (63),  d0oo2  =  c/2  and  d0oo4  =  c/4.  Substituting  these  values  into 
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(62)  and  letting  n  —  1  gives 


(0002): 

\  sm  6  J 
=  5.172  A, 

(0004):  c=(—\ 
\  sm  6  J 

=  5.172  A, 


/  1.5406  \ 

\sin(34.66/2) ) 


/  2(1.5406)  \ 
\sin(73.14/2)  J 


which  is  slightly  larger,  but  in  good  agreement,  with  the  undoped  GaN  c-axis  length 
of  5.166  A.  The  same  analysis  of  other  GaN:RE  samples  produced  similar  results. 

3.2  Sample  Surface  Preparation 

The  physical  dimensions  of  the  sample  holder  (~  8  mm  x  8  mm)  required  some 
thin  films  to  be  cut  using  a  South  Bay  Technologies  Model  650  low-speed  diamond 
wheel  saw.  The  GaN: RE  thin  film  samples  were  then  rinsed  in  ethanol  and  air- 
dryed.  Next,  electrical  contact  was  established  between  the  thin  film  sample,  tan¬ 
talum  mounting  plate,  and  molybdenum  sample  holder  using  several  spot  welds. 
Figure  16  shows  a  typical  configuration  of  the  above  mentioned  contacts.  The  sam¬ 
ple  holder  was  then  loaded  into  a  preparation  chamber  and  vacuum  conditions  to 
~  10~9  Torr  were  established,  which  typically  required  6-10  hours.  Next,  a  manipu¬ 
lator  arm  transferred  the  sample  holder  from  the  preparation  chamber  to  the  sample 
chamber  and  inserted  the  holder  into  the  copper  receiving  block.  Finally,  vacuum 
conditions  to  ~  10~10  Torr  were  established,  which  typically  required  less  than  one 
hour. 

Atomically  clean  sample  surfaces  were  obtained  by  several  sputtering  and  anneal¬ 
ing  cycles.  Typically,  sputtering  occurred  via  Ar+  ion  bombardment  at  ±45°  with 
respect  to  the  sample  normal  for  20  minutes  in  each  orientation.  Next,  resistive  heat- 
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(a)  Top  view  of  sample  holder. 


]  Molybdenum  sample  holder 


Tantalum  plate  /  strips 


]  Sample  film 


(b)  Side  view  schematic  of  sample  holder. 

Figure  16.  Sample  mounting  technique  used  at  the  Center  for  Advanced  Microstruc¬ 
tures  and  Devices. 


ing  to  a  temperature  of  approximately  500  K  via  a  tungsten  filament  located  behind 
the  sample  holder,  followed  by  gradual  return  to  room  temperature,  annealed  the 
sample  surface. 

3.3  Photoemission  Spectroscopy  Considerations 

Many  of  the  considerations  discussed  here  are  specific  to  the  experimental  setup 
at  the  Center  for  Advanced  Microstructures  and  Devices,  but  most  apply  to  the  XPS 
experiment  conducted  at  the  University  of  Nebraska  (Lincoln),  as  well. 

Every  photoemission  experiment  begins  with  a  photon  source.  At  the  CAMD 
synchrotron,  electron  packets  are  accelerated  through  an  underground  linear  accel¬ 
erator  into  a  current  storage  ring  inside  the  experimental  hall.  Within  the  storage 
ring,  magnetic  fields  impart  a  Lorentz  force  on  the  electrons,  which  are  subsequently 
turned  and  accelerated.  This  acceleration  results  in  photon  emission  perpendicular 
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to  the  direction  of  electron  propagation.  Next,  the  photons  travel  through  a  series  of 
optics,  which  produce  monochromatic,  linearly  polarized  light  that  is  used  to  irradi¬ 
ate  a  sample.  Figure  17  shows  a  cross  section  sketch  of  the  3-meter  toroidal  grating 
monochromator  (3m  TGM)  experimental  setup  used  for  all  experiments  conducted 
at  CAMD. 


Figure  17.  Schematic  of  the  3-meter  Toroidal  Grating  Monochromator  (3m  TGM) 
beamline  at  the  Center  for  Advanced  Microstructures  and  Devices  [4]. 

Sample  irradiation  results  in  photoelectron  emission  with  a  range  of  energies 
and  directions,  and  there  are  several  components  required  to  efficiently  filter  these 
electron  kinetic  energies,  as  depicted  schematically  in  Figure  18.  The  electron  optics, 
which  are  typically  a  set  of  electrostatic  or  magnetic  lenses,  allow  a  portion  of  the 
photoclectrons  at  the  correct  acceptance  angle  and  of  a  given  energy,  referred  to  as  the 
pass  energy,  to  arrive  at  the  analyzer  entrance  slit.  More  specifically,  photoclectrons 
of  the  specified  pass  energy  are  measured  by  setting  voltages  for  the  lens  system 
that  both  focus  electrons  onto  the  entrance  and  that  retard  the  electron  velocity 
so  that  their  kinetic  energy  after  passing  through  the  transfer  lenses  matches  the 
pass  energy  of  the  hemispherical  analyzer.  Therefore,  in  order  to  record  an  electron 
kinetic  energy  spectrum  over  a  range  of  energies,  it  is  necessary  to  scan  the  voltages 
applied  to  the  transfer  lenses,  referred  to  as  a  lens  function. 

Photoelectrons  that  pass  through  the  entrance  slit  travel  through  the  hemispheri¬ 
cal  electron  energy  analyzer,  subject  to  the  applied  1/r  field  between  the  two  concen- 
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trie  hemispheres,  and  depart  through  the  analyzer  exit  slit.  The  departing  electrons 
are  then  multiplied  via  an  electron  multiplier  and  counted  via  electronics  and  data 
collection  software. 

The  energy  resolution  of  the  analyzer  A E  is  dependent  on  the  slit  width  W,  mean 
radius  R,  acceptance  angle  cp  and  pass  energy  Ep  as  [5] 

(W  o?\ 

AE=E^2R  +  Y)'  <64> 

The  3m  TGM  configuration  is  equipped  with  a  50  mm  hemispherical  electron 
energy  analyzer  [4]  with  an  acceptance  angle  of  ±1°;  pass  energy  was  typically  set 
to  15  V.  The  operational  range  is  15-200  eV,  and  the  combined  resolution  of  the 
photoemission  configuration,  to  include  the  monochromator,  all  optics,  and  the  elec¬ 
tron  analyzer  is  about  70  meV  for  photon  energies  up  to  140  eV.  At  higher  photon 
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energies,  the  combined  beamline/analyzer  resolution  is  about  150  meV  [4], 

3.4  Photoemission  Spectra  Conditioning 

3.4.1  Conversion  to  Binding  Energy. 

Photoemission  spectroscopy  requires  establishing  a  meaningful  reference  level  to 
set  the  energy  scale.  Although  the  kinetic  energies  of  photoelectrons  are  measured, 
spectra  are  typically  displayed  in  binding  energies  relative  to  the  Fermi  edge  of  the 
sample.  In  a  semiconductor,  the  Fermi  level  is  not  the  best  reference  level,  as  it 
shifts  around  in  the  bandgap  at  the  surface  due  to  changes  in  band  bending  and 
doping  (n  or  p)  type.  However,  the  Fermi  level  of  a  metal  is  not  susceptible  to  these 
variations.  Thus,  the  accepted  practice  [6-9]  is  to  use  a  clean  metal  foil  to  determine 
the  conversion  from  the  measured  kinetic  energy  to  the  binding  energy,  as  depicted 
graphically  in  Figures  19-21. 


Figure  19.  Raw  spectra  of  intensity  versus  kinetic  energy  for  GaN:Yb  (black)  and 
tantalum  foil  (red)  using  90  eV  photons. 
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Figure  19  shows  the  raw  spectra  of  intensity  versus  kinetic  energy  for  GaN:Yb 
and  a  clean  tantalum  foil.  The  spectrum  for  tantalum  indicates  the  presence  of  a 
density  of  states  nearby  85  eV  that  resembles  the  Fermi-Dirac  distribution 

=  e(E-EF)/kBT  +  l  '  (65) 

Inspection  of  (65)  shows  that  f(Ep)  =  1/2.  Thus,  fitting  (65)  against  the  experi¬ 
mental  data  at  the  appropriate  temperature  and  such  that  f(E)  =  1/2  at  E  —  Ep 
provides  a  good  estimate  of  the  kinetic  energy  of  the  Fermi  level,  as  shown  in  Figure 
20  at  T  =  300  K.  The  kinetic  energy  of  the  Fermi  level  Ep  permits  shifting  from  a 


Figure  20.  Kinetic  energy  of  spectrometer  Fermi  level  using  90  eV  photons.  The  Fermi 
level  is  established  at  the  intersection  between  the  Fermi-Dirac  distribution  (65)  at 
T  —  300  K  (black)  and  the  experimentally  observed  tantalum  density  of  states  (red). 


kinetic  energy  to  a  binding  energy  scale  as 


Ep  —  -Ekin  —  Ep  , 


(66) 
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such  that  Eb  =  0  at  the  Fermi  level  with  all  other  binding  energies  expressed  as 
negative  values.  The  converted  GaN:Yb  spectrum  of  intensity  versus  binding  en¬ 
ergy  is  shown  in  Figure  21.  The  valence  band  maximum,  which  is  determined  by 


Figure  21.  Raw  spectrum  of  intensity  versus  binding  energy  for  GaN:Yb  using  90  eV 
photons.  The  valence  band  maximum  (VBM)  is  located  by  linear  extrapolation  of  the 
high  kinetic  energy  edge  of  the  spectrum  at  the  intersection  of  zero  background. 

linear  extrapolation  of  the  high  kinetic  energy  edge  of  the  spectrum,  is  located  at 
approximately  2.7  eV  below  the  Fermi  level  for  the  GaN:Yb  sample  shown.  Also, 
the  dominant  spectral  feature  located  at  a  binding  energy  of  approximately  —20.5 
eV  is  the  Ga  3 d  shallow  core  electronic  level.  Therefore,  the  separation  between  the 
valence  band  maximum  and  the  Ga  3 d  shallow  core  level  in  GaN:Yb  is  consistent 
with  the  accepted  value  of  17.7  eV  for  wurtzite  GaN  [8-10]. 

The  detailed  method  of  fitting  the  Fermi-Dirac  distribution  to  the  experimental 
data  is  shown  for  completeness.  In  practice,  visual  interpolation  of  the  Fermi  level  is 
sufficient  and  was  the  method  used.  Additionally,  the  kinetic  energy  corresponding 
to  the  Fermi  level  is  a  unique  photon  energy  dependent  value.  Thus,  in  theory,  the 
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Fermi  level  would  need  to  be  determined  for  every  photon  energy  used  during  an 
experiment.  In  practice,  however,  the  3m  TGM  beamline  historical  files  contain  a 
calibration  database  which  provides  the  Fermi  level  kinetic  energy  for  various  pho¬ 
ton  energies  throughout  the  operating  range,  such  that  values  not  included  in  the 
database  are  easily  interpolated.  During  this  research,  the  calibration  database  was 
both  verified  and  updated,  in  order  to  ensure  accuracy. 

3.4.2  Removal  of  Background  and  Secondary  Electron  Signals. 

After  converting  spectra  to  a  binding  energy  scale,  the  background  was  adjusted 
to  zero  and  secondary  electron  influence  was  removed.  Electrons  with  binding  ener¬ 
gies  above  EF  in  a  semiconductor  should  not  be  collected  during  photoemission,  so 
photoemission  intensities  at  values  greater  than  EF  were  subtracted  from  the  spectra. 
The  secondary  electron  signal  Is  influences  PES  spectra  by  increasing  exponentially 
with  decreasing  kinetic  energy,  as 

JS(KE)  =  Ae-(B'KE) ,  (67) 

where  A  and  B  are  fitting  parameters.  Therefore,  a  parametric  fit  for  (67)  can  be 
determined  by  using  several  intensity  values  at  locations  away  from  primary  spec¬ 
tral  features,  and  the  resulting  equation  can  be  subtracted  from  each  spectrum.  In 
practice,  this  was  accomplished  using  the  Origin  Baseline  Toolkit  with  6-10  values, 
dependent  upon  the  GaN:RE  sample  spectra,  and  the  resulting  baseline  was  sub¬ 
tracted  from  all  spectra  associated  with  the  specific  sample.  Figure  22  illustrates 
this  technique. 
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Figure  22.  Removal  of  secondary  electron  influence  using  Origin  Baseline  Toolkit. 
Several  intensity  values  (squares)  are  used  to  fit  a  baseline  (red)  via  (67),  which  is 
subtracted  from  the  sample  spectra. 

3.4.3  Spectra  Normalization. 

The  next  step  of  the  spectral  conditioning  process  requires  spectra  to  be  ad¬ 
justed  for  intensity  variations  caused  by  the  experimental  conditions.  The  electron 
ring  current  at  CAMD  experiences  an  exponential  decay  during  normal  operation, 
which  requires  regular  (typically  3  times  per  day)  electron  re-injections.  Hence,  the 
absolute  intensities  of  both  the  ring  current  and  the  photon  flux  at  the  3m  TGM 
beamline  endstation  vary  greatly  during  experiments.  However,  as  both  the  resonant 
photoemission  and  Schottky  barrier  experiments  required  comparison  of  the  relative 
intensities  of  many  spectra,  each  individual  spectrum  was  normalized  to  points  away 
from  the  major  spectral  features. 

Using  Figure  22  as  an  example,  the  saddle  between  the  valence  band  and  the  Ga 
3 d  shallow  core  level  was  identified  as  the  reference  level  for  that  spectrum.  Next,  the 
average  absolute  intensity  of  three  data  points  within  0.3  eV  was  used  to  calculate 
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the  intensity  of  the  spectrum  at  the  reference  point  as  Jrp;  this  process  was  repeated 
for  each  individual  spectrum.  The  largest  value  of  intensity  at  the  reference  point 
for  all  spectra  was  then  designated  as  /rP)  max.  Finally,  the  absolute  intensity  of  each 
individual  spectrum  J0  was  scaled  to  Is  as 

Is  =  R- 1 o  ;  R  =  /rprmax  (68) 

-*rp 

where  R  is  the  scaling  factor. 

The  technique  presented  here  for  Figure  22  applied  to  the  GaN:Yb  sample  only. 
Each  rare  earth  doped  GaN  thin  him  had  unique  spectral  features,  so  the  reference 
points  used  to  normalize  spectra  were  slightly  different  for  the  different  rare  earth 
dopants  studied,  as  would  be  expected.  Lastly,  the  intensities  presented  in  Chapters 
IV,  V,  and  VI  are  relative  intensities,  unless  stated  otherwise. 

3.4.4  Spectral  Feature  Deconvolution. 

The  final  step  of  PES  spectra  conditioning  is  to  deconvolute  the  features  of  each 
spectrum  into  their  subpeaks.  Photoemission  from  the  core  level  of  an  atom  produces 
a  Lorentzian  lineshape,  where  the  lifetime  of  the  state  determines  linewidth  [11],  In  a 
metal,  the  intrinsic  process  of  electron-hole  pair  production  via  photoemission  leads 
to  the  selection  of  an  asymmetric  Doniach-Sunjic  [12]  or  Mahan  [13]  lineshape  to 
represent  the  core  level.  However,  in  a  semiconductor  (insulator),  the  asymmetric 
lineshapes  are  difficult  to  observe  as  differential  charging  adds  a  symmetric  contribu¬ 
tion  to  linewidth  [11],  Moreover,  the  influence  of  Heisenberg’s  uncertainty  principle 
to  the  lifetime  of  a  state,  combined  with  the  effects  of  electron-phonon  coupling  fur¬ 
ther  broaden  the  core  level  and  yield  a  more  symmetric  lineshape.  Therefore,  the 
collective  effects  of  these  processes  on  GaN: RE  core  levels  indicated  a  least  squares 
fitting  of  Gaussian  line  shapes  as  the  logical  choice  to  represent  spectral  feature 
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subpeaks. 

Figures  23(a)  and  23(b)  show  the  Ga  3d  shallow  core  level  and  valence  band 
deconvolution  for  two  different  rare  earth  doped  GaN  samples.  The  full  width  at 


(a)  Deconvolution  of  Ga  3 d  shallow  core  level  (b)  Deconvolution  of  GaN:Er  valence  band 
from  GaN:Gd  photoemission  spectrum.  spectrum  using  multiple  Gaussian  subpeaks. 

Figure  23.  Deconvolution  of  valence  band  and  shallow  core  level  lineshapes  from 
example  GaN:RE  photoemission  spectra. 


half  peak  maximum  for  all  Gaussians  exceeded  the  overall  system  resolution  of  0.5  eV 
for  all  experiments.  Further  justification  for  the  number,  location,  and  identification 
(surface,  core,  alloy)  of  valence  band  and  shallow  core  level  subpeaks  for  the  different 
rare  earth  doped  GaN  thin  hints  will  be  discussed  in  Chapters  V  and  VI. 

Discussions  presented  here  and  in  Chapters  IV  and  V  provide  all  details  pertaining 
to  data  analysis  for  the  respective  experiments.  However,  a  few  special  considerations 
for  Schottky  barrier  height  measurement  require  additional  explanation. 
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3.5  Schottky  Barrier  Considerations 


Recall  from  Chapter  II  that  the  Schottky  barrier  height  at  the  metal-semiconductor 
interface  may  be  determined  experimentally  as 

®B,n  =  Eg  —  (Ep  —  Ey)  .  (69) 

During  metal  deposition,  alignment  of  the  metal  and  semiconductor  Fermi  levels  re¬ 
sults  in  an  upward  valence  band  bending,  as  shown  in  Figure  10  of  Chapter  II.  Band 
bending  may  be  expected  to  manifest  in  the  photoemission  spectra  as  a  shift  of  the 
valence  band  maximum  toward  the  Fermi  level  with  increased  metal  deposition,  from 
which  <f>B,n  can  be  calculated  via  (69).  However,  the  difficulty  in  interpreting  mixed 
metal-semiconductor  photoemission  spectra  lies  in  the  fact  that  as  metal  is  deposited 
onto  the  semiconductor  surface,  the  semiconductor  valence  band  edge  becomes  ob¬ 
scured  by  the  metal  overlayer  signal.  Surface  Fermi  level  movement  and  valence  band 
bending  were  determined  by  monitoring  the  Ga  3 d  core  level  shift  as  a  function  of  Au 
coverage.  Figure  24(a)  shows  the  Ga  3 d  core  level  and  the  valence  band  region  spec¬ 
tra  on  the  surface  for  a  clean  GaN:Gd  surface.  The  energy  between  these  two  levels, 
Evbm-Cj  is  considered  a  bulk  property  of  the  material  and  is,  therefore,  independent 
of  metal  coverage  [6,8,9].  When  metal  is  evaporated  on  the  sample  surface,  as  shown 
in  Figure  24(b),  the  edge  of  the  semiconductor  valence  band  maximum  is  obscured 
by  the  signal  originating  from  the  metal  overlayer.  Whereas  the  Fermi  level  of  the 
metal,  which  coincides  with  the  surface  Fermi  level  of  the  semiconductor,  can  now 
be  observed  (Figure  24(b)  inset),  the  valence  band  maximum  of  the  semiconductor 
is  no  longer  visible.  However,  owing  to  the  fact  that  Evbm-c  is  constant,  the  Ga  3d 
core  shift  results  in  an  equivalent  valence  band  bending  at  the  semiconductor  surface, 
from  which  the  Schottky  barrier  height  is  calculated  via  (69). 
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(a)  Photoemission  spectrum,  in  binding  energy,  for  a  clean  GaN:Gd  thin 
film  showing  the  Ga  3 d  core  level  peak  and  the  valence  band  region.  A 
least  squares  fit  is  used  to  determine  the  valence  band  maximum  (VBM) . 
The  value  of  Evbm-c,  shown  by  the  double  arrow,  is  determined  from 
these  features. 


(b)  Photoemission  spectrum  for  GaN:Gd  thin  film  with  ~  3  AAu  cov¬ 
erage.  The  Fermi  edge  of  the  metal  (inset)  indicates  the  Fermi  level 
position  at  the  semiconductor  surface.  The  core  level  shift,  as  com¬ 
pared  with  the  clean  spectrum,  determines  the  semiconductor  valence 
band  surface  band  bending. 


Figure  24.  Semiconductor  surface  valence  band  bending  via  core  level  shift. 
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IV.  The  Effective  Surface  Debye  Temperature  of  Yb:GaN 


The  contents  of  this  chapter  were  published  in  reference  [1].  The  effective  Debye 
temperatures  of  ytterbium  and  gallium  in  Yb:GaN  thin  films  have  been  obtained  us¬ 
ing  X-ray  photoemission  spectroscopy.  The  vibrational  motion  normal  to  the  surface 
results  in  a  dimunition  of  photoemission  intensities  from  which  we  have  estimated  the 
effective  Debye  temperatures  of  221  ±  30  K  and  308  ±  30  K  for  Yb  and  Ga,  respec¬ 
tively.  The  difference  between  the  measured  values  for  Yb  and  Ga  suggests  that  the 
Debye  temperatures  are  influenced  by  the  local  environment.  The  smaller  effective 
surface  Debye  temperature  for  Yb  correlates  to  a  soft,  strained  surface,  possibly  due 
to  an  increased  Yb-N  bond  length  as  compared  to  the  Ga-N  bond  length. 

4.1  Introduction 

During  the  past  decade,  rare  earth  doped  semiconductors  have  generated  consid¬ 
erable  attention  for  their  application  in  new  optoelectronic  devices  [2-5].  The  favor¬ 
able  thermal,  chemical,  and  electronic  properties  of  wide  band  gap,  Ill-nitride  semi¬ 
conductors  suggest  device  feasibility  using  lanthanide  doped  AIN  and  GaN.  More¬ 
over,  the  tunable  bandgaps  of  these  Ill-nitride  alloys  offer  device  applications  across 
the  visible  spectrum  through  the  ultra-violet  range,  to  include  optically  stimulated 
lasers  [6] ,  p  —  n  junction  light  emitting  diodes  [7],  and  thin  film  electroluminescence 
of  phosphors  [8-10]. 

Phonon  mediated  exciton  decay  can  lead  to  line  broadening  and  spectral  defor¬ 
mation.  While  phonons  can  affect  excited  state  lifetimes,  so  too  can  the  carrier 
concentration  and  local  structural  deformation.  A  signature  of  phonon  mediated 
processes  would  be  a  low  effective  Debye  temperature.  The  true  surface  Debye  tem¬ 
perature,  containing  the  in-plane  and  anharmonic  motions,  is  difficult  to  measure  in 
most  surface  spectroscopies  [11].  However,  the  effective  surface  Debye  temperature 
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is  readily  obtained  using  X-ray  photoemission  spectroscopy  (XPS)  and  ultraviolet 
photoemission,  low  energy  electron  diffraction,  inverse  photoemission  spectroscopy, 
and  other  surface  sensitive  techniques  [12-22],  The  experimentally  derived  effective 
surface  Debye  temperature  is  dominated  by  the  dynamic  motion  of  vibrational  modes 
normal  to  the  surface  and  tends  to  be  independent  of  surface  orientation  [11,19], 
Increased  thermal  vibrations  dictate  that  the  intensity  of  an  emitted  or  scattered 
electron  beam  decays  exponentially  with  increasing  temperature  as  [12-22] 

I  —  Iq  exp(— 2ID ) ,  (70) 


where  W  is  the  Debye-Waller  factor  and  is  given  by 

2W  =  |  A/c|2  (u0)2  ,  (71) 


where  A k  is  the  wave  vector  transfer  and  («o)2  is  the  mean  square  displacement  of 
the  atoms.  Within  the  Debye  model  of  thermal  vibrations,  in  the  case  of  isotropic 
vibrations,  W  is  described  as 


3h2{Ak)2T 
rn  ks  02d 


(72) 


where  h(Ak)  is  the  electron  momentum  transfer,  T  is  the  temperature  of  the  sample 
(in  Kelvin),  m  is  the  mass  of  the  scattering  center,  ks  is  the  Boltzmann  constant, 
and  @£>  is  the  effective  surface  Debye  temperature.  In  the  case  of  photoemission 
spectroscopy,  the  momentum  transfer  is  equal  to  the  momentum  of  the  emitted  elec¬ 
tron  [12-22] ,  and  the  scattering  center  mass  is  the  mass  of  the  specific  element  from 
which  the  emitted  photoelectron  originates  [17],  This  chapter  reports  the  effective 
surface  Debye  temperatures  of  Ga  and  Yb  via  XPS  spectra  taken  from  Yb:GaN  thin 
films  and  correlates  the  results  with  thin  him  strain. 
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4.2  Experimental 


Thin  films  of  Yb^Gai-^N  (50-300  nm)  were  fabricated  on  Si(lll)  substrates  by 
RF  plasma  (EPI  620)  assisted  molecular  beam  epitaxy  (MBE).  The  growth  param¬ 
eters  for  the  deposition  of  Yb  doped  (in  situ )  GaN  thin  films  were  base  pressure  of 
~  10~n  Torr,  nitrogen  flux  of  0.75-1.0  SCCM,  RF  power  of  500  W,  substrate  tem¬ 
perature  of  850-900  °C,  Ga  cell  temperature  of  850  °C,  and  Yb  cell  temperature 
of  500-850  °C.  Thickness  of  the  films  was  measured  with  a  surface  profilometer. 
The  orientation,  crystal  structure,  and  phase  purity  of  the  films  were  established 
by  Cu  Ka  (A=1.5406  A  radiation  X-ray  diffraction,  as  shown  in  Figure  25,  using  a 
Siemens  D5000  X-ray  diffractometer. 
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Figure  25.  The  X-ray  diffraction  (XRD)  pattern  of  a  GaN:Yb  film  grown  on  Si(lll) 
substrates.  The  films  show  c-axis  orientation  of  the  texture  growth  and  a  high  level 
of  crystallinity.  The  presence  of  any  secondary  phases  or  spurious  peaks  has  not  been 
observed. 


Surface  adsorbates  were  removed  by  several  preparatory  sputtering  and  thermal 
annealing  cycles  in  ultra  high  vacuum  (UHV).  The  resulting  surfaces  were  charac¬ 
terized  by  core  level  XPS  and  found  to  be  free  of  oxygen.  The  XPS  studies  were 
performed  in  an  UHV  chamber  with  a  hemispherical  electron  analyzer  (Thermo  VG 
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Scientific  VG100AX).  The  combined  resolution  of  the  system  was  about  1.1  eV  using 
a  Mg  Ka  line  (1253.6  eV)  from  a  fixed  anode  X-ray  source.  Surface  temperature  was 
controlled  by  a  combination  of  resistive  heating  and  cooling  with  liquid  nitrogen. 
The  temperature  was  monitored  with  a  W-5%Re/W-26%Re  thermocouple  with  an 
accuracy  of  5  K.  The  photoemission  spectra  were  taken  with  a  45°  incidence  angle 
and  normal  emission,  with  binding  energies  referenced  to  the  Fermi  level  of  a  clean 
Ta  foil. 


4.3  Results  and  Discussion 

As  anticipated,  the  core  level  photoemission  intensities  decrease  with  increasing 
temperature,  as  shown  in  Figures  26(a)  and  26(b).  Figure  27  shows  the  logarithm 
of  the  core  peak  photoemission  intensities  for  Ga  3 p  and  Yb  4 d  as  a  function  of 
temperature.  The  kinetic  energy  of  the  outgoing  photoelectron  is  calculated  as 


Akin  =  hv  -  \Eb\  -  4>a  , 


(73) 


where  hv  is  the  incident  photon  energy,  \Eb\  is  the  binding  energy  of  the  emitted 
photoelectron,  and  (j)a  is  the  electron  analyzer  work  function  (4.5  eV).  Thus,  A k 
values  in  Figure  27  are  calculated  from 

(A kf  =  ,  (74) 


where  me  is  the  electron  mass.  From  (70)  and  (72),  the  slope  S  for  each  XPS  data 
set  in  Figure  27  determines  ©£>  from 


3H2(Ak)2 

m  kp,  (~)  2d 


(75) 
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(a)  Ga  3 p  (b)  Yb  4 d 


Figure  26.  Temperature  dependent  X-ray  photoemission  spectra  of  the  indicated  core 
levels  after  background  subtraction. 


The  fact  that  the  effective  Debye  temperature  is  different  for  ytterbium  and 
gallium  in  a  Yb:GaN  thin  film  indicates  that  this  effective  Debye  temperature  is 
influenced  by  the  local  environment,  not  exclusively  by  the  overall  rigidity  of  the 
thin  film  lattice.  The  difference  between  the  measured  effective  Debye  temperatures 
for  Ga  and  Yb  provides  an  indication  of  the  surface  rigidity  and  lattice  strain.  In 
a  simple  picture,  examination  of  equations  (71)  and  (72)  predicts  an  inverse  rela¬ 
tionship  between  atomic  displacement  and  surface  Debye  temperature  such  that  an 
increased  0 o  predicts  lower  atomic  displacement,  or  a  “rigid”  surface.  It  follows 
that  a  decreased  0£>  is  indicative  of  a  “soft”  surface.  Thus,  our  respective  values 
of  221  ±  30  K  and  308  ±  30  K  for  the  effective  Yb  and  Ga  Debye  temperatures  in 
Yb:GaN  reflect  softening  from  the  introduction  of  Yb  atoms  into  the  GaN  wurtzite 
crystal  structure. 

If  our  measured  effective  Debye  temperatures  for  Yb  and  Ga  were  identical,  within 
experimental  uncertainty,  it  would  be  indicative  of  a  substitutional  occupation  of  a 
Ga  site  by  a  Yb  ion.  Although  the  measured  effective  Debye  temperatures  were 
not  identical,  the  experimental  uncertainty  indicates  that  the  Yb  and  Ga  Debye 
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Figure  27.  Logarithm  of  the  core  peak  photoemission  intensities  for  Yb  and  Ga  as  a 
function  of  temperature  following  background  (/bg)  subtraction  and  normalization  to 
the  peak  intensity  at  the  lowest  temperature  Ig.  The  experimental  values  were  fitted 
(solid  line)  with  the  Debye- Waller  factor. 


temperatures  could  differ  by  as  little  as  27  K.  This  similarity  between  the  Yb  and 
Ga  Debye  temperatures  suggests  that  Yb  doping  does  occur  by  Ga  site  substitution. 
Experimental  data  and  crystal-field  calculations  have  confirmed  that  rare  earth  ions 
in  a  GaN  host  occupy  relaxed  substitutional  Ga  sites  [2,3].  To  explain  the  slight 
difference  in  the  Debye  temperatures,  strain  should  be  considered. 

From  the  standpoint  of  bond  length,  shorter  bonds  have  higher  binding  energies 
than  longer  ones,  which  implies  that  decreasing  the  bond  distance  reflects  stiffening 
and  increasing  bond  distance  reflects  softening.  Thus,  the  slightly  lower  calculated 
effective  Debye  temperature  for  Yb  in  Yb:GaN  reveals  a  softened  surface,  and  we 
contend  that  this  softening  arises  from  a  lengthening  of  the  Yb-N  bond  distance,  as 
compared  with  the  Ga-N  distance,  which  would  manifest  as  thin  film  strain.  The 
decrease  in  Debye  temperature  should  be  reflected  in  a  decrease  in  exciton  lifetime 
as  this  means  that  new  phonon  decay  channels  are  opened.  This  in  turn  should  lead 
to  an  increase  in  fluorescence  peak  widths. 

From  Figure  25,  the  c-axis  length  of  our  Yb:GaN  thin  films  was  found  to  be 
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approximately  5.172  A,  which  is  slightly  larger  than  the  widely  reported  and  accepted 
c-axis  length  of  undoped  GaN  (5.166  A).  The  XRD  data  arc  limited  because  the  films 
are  so  very  nearly  crystalline,  that  is  to  say  highly  textured.  There  will  be  changes  to 
the  a-  and  6-axis  lengths  in  response  to  contraction  or  transverse  strain,  as  dictated 
by  Poisson’s  ratio.  However,  the  implementation  of  Poisson’s  ratio  depends  on  the 
details  of  the  stress-strain  tensor,  which  cannot  be  determined  accurately  from  our 
experimental  data.  Thus,  the  values  for  the  a-axis  and  6-axis  cannot  be  accurately 
established  from  our  data,  but  would  require  additional  measurements  like  grazing 
incidence  X-ray  scattering.  However,  the  small  increase  in  c-axis  lattice  constant, 
and  hence  thin  film  strain,  can  be  explained  due  to  the  fact  that  the  Yb3+  ion  has  a 
larger  ionic  radius  than  the  Ga3+  ion  by  at  least  40%. 

The  tetrahedral  arrangement  of  four  nearest  neighbors  in  wurtzite  GaN  yields 
a  Ga3+  ionic  radius  of  0.47  A  according  to  [23],  and  others.  Regarding  Yb  doping 
in  GaN,  references  have  not  been  found  for  the  Yb3+  ionic  radius  in  the  same  con¬ 
figuration.  However,  the  available  data  indicate  ionic  radii  for  Yb3+  ranging  from 
0.87-1.04  A  for  a  coordination  number  of  VI-1X.  Similarly,  the  range  of  ionic  radii  for 
Ga3+  are  0.47-0.62  A  for  a  coordination  number  of  IV- VI.  Therefore,  in  the  absence 
of  a  Yb3+  data  point  with  coordination  =  4,  then  the  most  conservative  comparison 
using  the  largest  possible  documented  ionic  radius  for  Ga3+  (0.62  A)  and  the  small¬ 
est  possible  documented  ionic  radius  for  Yb3+  (0.87  A)  further  supports  the  claim  of 
doping-induced  thin  him  strain  due  to  increased  Yb-N  bond  length,  compared  with 
the  Ga-N  bond  length. 

4.4  Conclusions 

X-ray  photoemission  spectroscopy  has  been  used  to  determine  the  effective  Debye 
temperatures  of  gallium  and  ytterbium  in  Yb:GaN  thin  films.  The  similarity  between 
their  measured  values  suggests  that  substitutional  occupation  of  a  Ga  site  by  a  Yb 
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ion  occurs,  which  supports  experimental  data  and  theoretical  calculations.  Careful 
comparison  of  the  Ga  and  Yb  values  indicates  that  the  slightly  smaller  effective  Debye 
temperature  of  the  Yb  atom  reflects  a  surface  that  is  softened  and  strained,  possibly 
due  to  an  increased  Yb-N  bond  length  as  compared  to  the  Ga-N  bond  length. 
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V.  Resonant  Photoemission  of  Rare  Earth  Doped  GaN 

Thin  Films 


The  contents  of  this  chapter  were  accepted  for  publication  in  reference  [1].  The 
4 d  — >■  4/  Fano  resonances  for  various  rare  earth  doped  GaN  thin  films  (RE  =  Gd,  Er, 
Yb)  were  investigated  using  synchrotron  photoemission  spectroscopy.  The  resonant 
photoemission  Fano  profiles  show  that  the  major  Gd  and  Er  rare  earth  4/  weight  is 
at  about  5-6  eV  below  the  valence  band  maximum,  similar  to  the  4/  weights  in  the 
valence  band  of  many  other  rare  earth  doped  semiconductors.  For  Yb,  there  is  very 
little  resonant  enhancement  of  the  valence  band  of  Yb  doped  GaN,  consistent  with 
a  largely  4/14  occupancy. 

5.1  Introduction 

During  the  past  decade,  rare  earth  doped  semiconductors  have  generated  consid¬ 
erable  attention  for  their  application  in  new  optoelectronic  devices  [2-5].  Ill-nitride 
semiconductors,  such  as  AIN,  GaN,  and  InN  offer  tunable  bandgaps  and  favorable 
thermal,  chemical,  and  electronic  properties,  which  facilitate  various  device  applica¬ 
tions  [6-10]  from  the  ultraviolet  through  the  visible  spectrum  to  the  infrared  range. 
Moreover,  thin  film  electroluminescent  phosphors  with  red,  blue,  and  green  emis¬ 
sions  [8-16]  imply  the  promise  of  full  color  (white)  light  capability.  Rare  earth  doping 
GaN  might  have  a  number  of  advantages:  there  is  the  promise  that  Eu  or  Er  doping 
will  improve  the  light  output.  Luminescence  due  to  the  Er  intra-4/-shell  transition 
from  the  4Ii3/2  excited  state  to  the  4Iis/2  ground  state  is  known  to  be  particularly 
intense  and,  above  all,  efficient.  Due  to  their  highly  localized  4/  electrons  [17,18],  the 
direct  /  —  /  interactions  between  the  neighboring  rare  earth  atoms  are  very  weak  and 
nonexistent  in  a  weakly  doped  semiconductor  host.  In  addition,  rare  earth  doping  is 
seen  to  increase  the  gold  electrode  Schottky  barrier  heights  significantly  [19]  thereby 
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decreasing  leakage  currents  in  particle  type  detector  devices. 

While  numerous  resonant  photoemission  studies  of  rare  earth  metals  have  been 
reported,  the  resonant  processes  resulting  from  photon  interaction  with  Ill-nitride 
semiconductors  are  significantly  less  understood.  Plucinski  et  al.  [20]  reported  a  reso¬ 
nant  photoemission  process  at  the  Ga  3 p  absorption  threshold  in  GaN  and  compared 
the  results  to  those  reported  for  GaP  and  GaAs  [21,22],  Lastly,  Maruyama  et  al.  [23] 
reported  on  Eu  doped  GaN  using  X-ray  photoelectron  spectroscopy  (XPS)  and  res¬ 
onant  photoemission  spectroscopy  (RPES)  and  concluded,  via  4 d  — >  4/  resonant 
photoemission  measurements,  that  the  transition  from  trivalent  to  divalent  Eu  ions 
occurred  near  the  surface  of  GaN.  The  key  value  of  resonant  photoemission,  how¬ 
ever,  is  to  probe  what  valence  bands  of  the  semiconductors  have  strong  4/  and/or 
rare  earth  weight  [24-26].  Of  course,  such  photon  energy  dependent  studies  must 
be  disentangled  from  bulk  band  structure  effects,  so  there  is  considerable  value  in 
studying  the  resonant  photoemission  process  of  a  semiconductor  with  several  differ¬ 
ent  rare  earths;  each  is  likely  to  dope  the  semiconductor  in  a  similar  fashion,  but  the 
resonant  enhancement  of  the  valence  band  will  occur  at  different  photon  energies. 

Although  likely  to  locally  strain  the  lattice,  the  4/  rare  earths  will  tend  to  adopt 
substitutional  sites  for  Ga  [2,3,27]  in  GaN  while  significantly  altering  magnetic  [28- 
35]  and  optical  properties  [8-15],  and  it  is  therefore  of  considerable  interest  to  know 
whether  even  low  concentrations  of  a  rare  earth  in  the  GaN  host  can  alter  the  surface 
electronic  structure.  This  is  likely,  as  although  rare  earths  are  isoclectronic  with 
Ga3+,  they  may  be  associated  with  other  defects  [36,37].  With  these  considerations 
in  mind,  we  have  engaged  in  investigations  of  the  surface  electronic  structure  and 
interface  properties  of  the  RE.xGai_xN  (RE  =  Gd,  Er,  Yb)  semiconductors. 


5.2  Experimental 


The  REiTGai_a;N  thin  films  (50-300  nm)  were  fabricated  on  Si(lll)  (RE  =  Gd, 
Yb)  and  sapphire  A1203  (RE  =  Er)  substrates  by  RF  plasma  (EPI  620)  assisted 
molecular  beam  epitaxy  (MBE).  The  growth  parameters  for  the  deposition  of  RE- 
doped  (in  situ)  GaN  thin  films  were  base  pressure  of  ~  10_n  Torr,  nitrogen  flux 
of  0.75-1.0  SCCM  (Gd,  Yb)  and  2.0  SCCM  (Er),  RF  power  of  500  W,  substrate 
temperature  of  850-900  °C,  Ga  cell  temperature  of  850  °C,  and  RE  cell  temperatures 
of  1050-1100  °C  (Gd),  1000-1100  °C  (Er),  and  500-850  °C  (Yb).  The  thickness  of 
the  films  was  measured  with  a  surface  profilometer  and  atomic  force  microscopy. 

The  orientation,  crystal  structure,  and  phase  purity  of  the  films  were  established 
by  Cu  Ka  (A  =  1.5406  A)  radiation  X-ray  diffraction  using  a  Siemens  D5000  X-ray 
diffractometer.  The  X-ray  diffraction  (XRD)  patterns  of  Gd,  Er,  and  Yb  doped  GaN 
films  show  c-axis  orientation  and  a  high  degree  of  crystallinity.  The  presence  of  any 
secondary  phases  or  spurious  peaks  has  not  been  observed,  as  described  elsewhere 
[19].  Slight  shifts  in  diffraction  peak  positions  towards  lower  Bragg  angles  has  been 
observed  with  Er  doped  GaN  grown  on  Al203(0001)  substrates  and  RE^Ga^^N  thin 
films  (50-300  nm)  fabricated  on  Si(lll)  (RE  =  Gd,  Yb),  which  is  indicative  of  some 
lattice  expansion,  as  is  expected  [27].  The  c-axis  length  of  GaN:Yb  was  found  to  be 
5.172  A  [27],  which  is  very  close  to  the  widely  reported  and  accepted  c-axis  length 
(5.166  A)  of  undoped  GaN. 

The  elemental  compositions  of  the  rare  earth  doped  GaN  thin  films  grown  under 
different  conditions  were  characterized  by  energy  dispersive  spectroscopy  (EDS)  and 
a  VG  Microtech  XPS  attached  to  the  MBE  growth  system  (VG  Microtech).  The 
measured  concentrations  were  found  to  be  at  1-2%,  as  confirmed  from  the  Ga  2 p3/2, 
Er  Ad,  Gd  4 d,  Yb  4 d,  and  N  Is  core  level  XPS  intensities  using  an  A1  Ka  (1486.8 
eV)  X-ray  source.  The  typical  values  for  Er  concentrations  were  found  to  be  ~  5%, 
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higher  than  the  EDS-  and  XPS-derived  Gd  and  Yb  concentrations.  In  the  rare  earth 
doped  GaN  samples,  surface  segregation  cannot  be  excluded  and  may  well  be  likely, 
at  least  in  the  selvedge  region  of  the  surface. 

The  photoemission  experiments  were  conducted  on  the  3m  TGM  beamline  [38] 
at  the  Center  for  Advanced  Microstructures  and  Devices  at  Louisiana  State  Uni¬ 
versity  [39-41],  The  beamline  is  equipped  with  a  photoemission  endstation  with  a 
50  mm  hemispherical  electron  energy  analyzer,  with  a  resolution  of  about  70  meV, 
as  described  elsewhere  [38,42],  Photoemission  spectra  were  taken  with  a  45°  inci¬ 
dence  angle  and  the  photoclectrons  collected  along  the  sample  normal.  All  spectra 
presented  are  normalized  to  the  photon  flux,  and  the  secondary  electron  background 
has  been  subtracted.  The  position  of  the  Fermi  level  was  established  using  a  clean  Ta 
foil  as  reference.  All  binding  energies  reported  here  are  with  respect  to  this  common 
Fermi  level  in  terms  of  E  —  Ep,  so  that  occupied  state  binding  energies  are  negative. 
Energy  distribution  curves  (EDCs)  were  obtained  by  fixing  the  photon  energy  hv 
and  sweeping  electron  kinetic  energy  Akin,  thus  measuring  binding  energies.  Con¬ 
stant  initial  state  spectra  were  obtained  by  simultaneously  sweeping  hv  and  in,  so 
as  to  hold  binding  energy  fixed. 

Atomically  clean  GaN: RE  surfaces  were  obtained  by  several  preparatory  cycles 
of  Ar+  ion  sputtering  and  annealing,  as  described  elsewhere  [19,27].  This  will  create 
a  number  of  point  defects,  but  photoemission  is  generally  insensitive  to  such  defects. 
The  photoemission  spectra  from  the  clean  sample  surfaces  indicated  that  the  surfaces 
were  free  of  contaminants. 
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5.3  Results  and  Discussion 


5.3.1  The  4/  Contributions  to  the  Valence  Band  of  Doped  GaN. 

The  effects  of  band  hybridization  and  dopant  induced  strain  between  the  various 
rare  earth  dopants  and  the  GaN  surfaces  are  not  identical.  The  valence  band  spectra 
for  Gd  and  Er  doped  GaN  are  very  similar  as  seen  in  Figures  28  and  29,  while  for  Yb 
doped  GaN  (Figure  30),  the  valence  band  is  significantly  broader  extending  from  —2 
eV  to  nearly  —14  eV.  The  spectra  do  change  with  with  photon  energy  and  significant 
enhancements  of  some  of  the  valence  band  features  are  seen  at  some  photon  energies. 
Figures  28,  29,  and  30  (RE  =  Gd,  Er,  Yb)  show  the  valence  band  photoemission 
spectra  at  various  photon  energies  in  the  vicinity  of  the  respective  rare  earth  dopant 
4 d  absorption  thresholds.  For  the  rare  earth  dopants  studied,  resonance  results  from 
a  signal  overlap  between  direct  emission  of  photoelectrons  from  the  4/  state 

4d104fN  +  hv  — >  4d104/Ar_1  +  e“  ,  (76) 

and  Auger-like  electrons  emitted  in  a  super  Coster-Kronig  process  [43] 

4d104/Ar  +  hv  — >  [AdUfN+1]*  — >  4d104/7V_1  +  e~  ,  (77) 

where  [  ]*  denotes  an  excited  state.  The  final  states  for  both  the  direct  and  recom¬ 
bination  processes  are  identical. 

For  Gd  doped  GaN  (GaN:Gd),  Figure  28  shows  that  the  intensity  of  the  spectral 
features  near  8-9  eV  below  the  Fermi  level  first  increase  slightly  at  ~  144  eV,  followed 
by  a  rapid  and  considerable  enhancement,  with  maximum  intensity  at  a  photon 
energy  of  ~  148  eV.  Figure  29  shows  a  similar,  albeit  not  identical,  response  for 
GaN:Er.  As  with  the  GaN:Gd  thin  film,  the  spectral  feature  at  an  approximate 
binding  energy  of  —9  eV  in  GaN:Er  resonates  in  the  range  of  the  Er  4 d  absorption 
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Figure  28.  Energy  distribution  curves  of  GaN:Gd  for  various  incident  photon  energies, 
as  shown  for  each  spectrum.  The  strong  surface  and  bulk  4/  Gd  components  of  the 
GaN:Gd  valence  band  are  illustrated.  Resonating  components  C  (—8.1  eV)  and  D 
(  —  9.6  eV),  used  for  the  constant  initial  state  curves  in  Figure  31,  are  shown  at  the 
bottom. 
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Figure  29.  Energy  distribution  curves  of  GaN:Er  for  various  incident  photon  energies, 
as  shown  for  each  spectrum.  The  resonating  components  C  (—8.9  eV)  and  D  (  —  10.7 
eV),  used  for  the  constant  initial  state  curve  in  Figure  32,  are  shown  at  the  bottom. 
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Figure  30.  Energy  distribution  curves  of  GaN:Yb  for  various  incident  photon  energies, 
as  shown  for  each  spectrum.  Resonating  components  C  (—6.2  eV)  and  D  (  —  7.6  eV), 
used  for  the  constant  initial  state  curve  in  Figure  33,  are  shown  at  the  bottom. 


threshold.  However,  Figure  29  shows  that  the  features  near  the  bottom  of  the  valence 
band  increase  in  intensity  considerably  first  at  ~  165  eV  and  then  again  at  ~  176 
eV  for  GaN:Er.  Lastly,  Figure  30  (GaN:Yb)  reveals  that  the  intensity  of  spectral 
features  at  binding  energy  ~  7  eV  increase  slightly  at  ~  178  eV,  and  slightly  again 
at  ~  190  eV. 

The  intensity  of  the  valence  band  features  for  which  the  direct  and  super  Coster- 


Kronig  processes  overlap  is  described  as  a  function  of  photon  energy  hv  by  the  Fano 


lineshape  [43,44]  as 


N{hv)  = 


(e  +  q? 

e2  +  1  ’ 


hu  —  huj 


where  huj  is  the  photon  energy  equal  to  the  core  level  absorption  threshold  binding 
energy,  q  is  a  fitting  parameter  for  the  core  level,  and  Y  =  A(h,Uj)/2  where  A (hvj) 
is  the  width  (FWHM)  of  the  core  level. 
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5.3.2  The  4/  Contributions  to  the  Valence  Band  of  GaN:Gd. 


The  GaN:Gd  valence  band  components,  at  binding  energies  in  the  region  of  —8 
to  —10  eV,  represent  the  valence  band  features  that  strongly  resonate  at  photon 
energies  in  the  vicinity  of  147  eV,  as  plotted  in  Figure  31.  These  intensity  reso¬ 
nances  (Figure  31)  are  very  similar  to  those  observed  for  Gd203  [45]  and  Gd  doped 
Hf02  [24,45].  Not  only  are  the  resonant  photon  energies  similar,  but  the  features 
at  the  bottom  of  the  valence  band  contain  a  strong  feature  at  about  —8  eV  bind¬ 
ing  energy  and  a  weaker  shoulder  in  the  vicinity  of  —10  eV  binding  energy.  This 
suggests  a  surface  and  bulk  component  for  the  Gd  in  GaN,  in  spite  of  the  low  Gd 
concentrations.  Figure  31  displays  the  constant  initial  state  spectra  and  the  calcu¬ 
lated  Fano  profiles  for  GaN:Gd.  The  weighted  oscillator  strength,  gf ’,  relates  directly 
to  the  electric  dipole  transition  probability  and,  hence,  line  intensity.  Theoretical 
interpretations  of  resonant  photoemission  in  lanthanide  metals  at  the  4 d  absorption 
threshold  via  gf  calculations  have  been  reported  [46]  and  have  been  shown  to  agree 
well  with  experimental  results  [47].  The  resonance  peak,  at  roughly  147.5  eV  for  the 
GaN:Gd  Fano  line  profile  in  Figure  31,  illustrates  reasonable  agreement  with  the  Gd 
gf  calculations  of  Sugar  [46],  which  indicated  a  strong  absorption  peak  centered  at 
149.0  eV.  This  agreement  supports  further  the  assertion  that  the  GaN:Gd  spectral 
features  at  the  bottom  of  the  valence  band,  in  the  region  of  —8  to  —10  eV  in  Fig¬ 
ure  28,  are  of  strong  Gd  4/  weight  or  represent  bands  that  strongly  hybridize  with 
the  rare  earth.  This  is  expected  based  on  the  Gd  occupied  4/  level  placement  in 
GdN  [28,48], 

5.3.3  The  4/  Contributions  to  the  Valence  Band  of  GaN:Er. 

Similarly,  for  GaN:Er,  the  valence  band  features,  again  at  the  bottom  of  the  va¬ 
lence  band  at  binding  energies  of  roughly  —8.9  and  —10.7  eV,  strongly  resonate  at 
photon  energies  of  166  and  173  eV,  as  plotted  in  Figure  32.  Two  Fano  line  shapes 
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Figure  31.  Constant  initial  state  curves  and  Fano  fit  (solid  lines)  for  feature  C  (squares) 
at  —8.1  eV  and  feature  D  (triangles)  at  —9.6  eV  in  GaN:Gd.  The  fitting  component, 
q,  for  feature  C  was  calculated  as  1.91;  for  feature  D  it  was  2.79. 

were  used  to  fit  the  experimental  results  for  GaN:Er,  shown  in  Figure  32.  The 
GaN:Er  Fano  line  profile  in  Figure  32  also  shows  reasonable  agreement  with  the 
gf  calculations  [46].  Specifically,  Sugar  reported  the  separation  between  oscillator 
strength  maxima  at  167.2  eV  and  174.8  eV  to  be  7.6  eV  for  Er  [46].  Our  experimen¬ 
tally  observed  separation  is  roughly  7  to  7.5  eV  between  absorption  peaks  at  166  eV 
and  173  eV  for  GaN:Er.  The  ratio  of  the  Er  computational  peak  maxima  reported 
by  Sugar  is  ~  1  :  2.6  [46],  somewhat  larger  that  the  ratio  of  the  peak  maxima  at 
roughly  166  eV  and  173  eV  of  about  1.552/2.142  or  ~  1  :  1.9  for  GaN:Er,  as  plotted 
in  Figure  32.  Although  inspection  of  equation  (78)  predicts  that  when  hu  — >  huj, 
N(hv )  — >  q 2,  nonetheless,  the  resonant  photoemission  of  the  GaN:Er  thin  him  at  ab¬ 
sorption  thresholds  of  roughly  166  eV  and  173  eV  represents  reasonable  agreement 
with  the  Er  4 d  shallow  core  (167.2  eV,  174.8  eV)  binding  energies.  Thus,  unlike  the 
case  of  the  Gd  doped  GaN,  there  are  spin-orbit  interactions  that  separate  out  the 
4d3/2  and  4 d5/2  in  GaN:Er. 

While  there  are  indications  that  the  occupied  density  of  states  would  have  the 
spectral  weight  density  at  the  bottom  of  the  valence  band,  although  complexed  with 
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Figure  32.  Constant  initial  state  curve  (triangles)  and  Fano  fit  (solid  lines)  for  the 
valence  band  feature  in  the  vicinity  of  —9  eV  binding  energy,  for  GaN:Er.  The  fit¬ 
ting  component,  q,  in  the  region  of  the  Er  4 d5/2  absorption  threshold  (166  eV)  was 
calculated  as  1.55;  in  the  region  of  the  Er  4d3/2  absorption  threshold  (173  eV)  it  was 
2.14. 


a  vacancy  defect,  contributions  at  the  top  of  the  valence  band  might  also  be  expected 
[37].  In  fact,  the  /  removal  energies  calculated  for  Er  placed  in  GaN  are  seen  to  be 
very  similar  to  those  noted  for  Gd  in  GaN  [49] .  This  differs  from  the  relative  position 
of  the  Gd  and  Er  4/  levels  in  GaAs  [49].  Judging  by  ErAs  [50],  we  would  expect 
4/  contributions  at  both  the  bottom  and  the  top  of  the  valence  band  of  GaN.  The 
region  of  4-6  eV  binding  energy  should  be  dominated  by  the  4/7/2  (mainly  the  Tg 
and  5I7  multiplet  components  [50-52]),  while  the  bottom  of  the  valence  band  should 
be  dominated  by  the  4 f5/2  (mainly  the  3M10  and  3L9  multiplet  components  [50-52], 
although  there  are  many  others).  Strong  hybridization  of  the  4/  states,  possibly 
modified  by  intra-atomic  /  —  d  and  /  —  s  hybridizations  [50,  53-55]  can  lead  to  a 
much  more  delocalized  core  exciton,  and  a  decrease  in  the  resonant  photoemission 
intensities  [56]. 
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5.3.4  The  4/  Contributions  to  the  Valence  Band  of  GaN:Yb. 


To  test  the  possibility  of  intra-atomic  f  —  d  and  /  —  s  hybridizations  [50,  53- 
55],  or  extra-atomic  4/  hybridization  with  the  Ga  and  N,  we  also  looked  at  the 
resonant  enhancement  of  the  valence  band  of  GaN:Yb,  as  seen  in  Figures  30  and 
33.  Although  much  weaker  than  the  resonant  enhancements  observed  for  GaN:Gd 
(Figure  31)  and  GaN:Er  (Figure  32),  there  is  a  resonant  enhancement  of  the  valence 
band,  particularly  in  the  region  of  approximately  —6  to  —7  eV  binding  energy,  as 
plotted  in  Figure  33.  Although  Yb  metal  should  have  an  electronic  configuration 


Photon  Energy  (eV) 


Figure  33.  Constant  initial  state  curve  (squares)  and  Fano  fit  (solid  lines)  for  the 
valence  band  features  in  the  vicinity  of  —7  eV  binding  energy,  for  GaN:Yb.  The 
fitting  component,  q,  in  the  region  of  the  Yb  4d5/2  absorption  threshold  (179  eV)  was 
calculated  as  0.33;  in  the  region  of  the  Yb  4 d3/2  absorption  threshold  (189  eV)  it  was 
1.15. 


of  4d104/14,  such  that  an  exited  4 d  electron  has  no  unfilled  4/  state  to  occupy, 
partial  4/  occupancy  is  still  possible.  Previous  experimental  results  [57,  58]  using 
the  oxidation-induced  valence  change  of  Yb  (4/14  — »  4/13)  determined  that  the  data 
for  Yb203  were  well  described  by  Fano’s  theory  of  interaction  between  discrete  and 
continuum  states.  There  are  similar  expectations  for  Er  doped  GaSb  [59].  Despite 
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the  relative  weakness  of  the  observed  Fano  line  shapes  for  GaN:Yb,  we  note  that 
the  enhancements  in  spectral  features  in  the  valence  band  do  exhibit  photoemission 
resonances  at  photon  energies  in  the  vicinity  of  179  eV  and  189  eV.  These  absorption 
thresholds  agree  reasonably  well  with  the  Yb  4d.5/2  (182.0  eV)  and  the  Yb  4^3/2  (190.8 
eV)  shallow  core  level  binding  energies.  Because  it  is  not  just  the  valence  band 
maximum  that  exhibits  the  weak  photoemission  resonance  in  the  expected  location 
of  the  Yb  4/  multiplets  [51],  the  depletion  of  the  4/  occupancy  appears  related  to 
strong  hybridization  with  the  GaN  lattice. 

5.4  Summary  and  Conclusions 

Resonant  photoemission  for  GaN  thin  films  doped  with  various  rare  earths  (RE 
=  Gd,  Er,  Yb)  was  investigated  using  synchrotron  based  photon  energy  dependent 
photoemission  spectroscopy.  The  calculated  Fano  profiles  predict  the  experimentally 
observed  4 d  —>  4/  super  Coster-Kronig  processes  and  strong  rare  earth  4/  weights 
were  observed  about  5-6  eV  below  the  valence  band  maximum  for  Gd  and  Er  doped 
GaN.  For  Yb,  the  photoemission  resonance  is  much  weaker,  indicating  that  there 
is  only  a  very  little  depletion  of  the  4/14  occupancy,  but  strong  hybridization  with 
GaN  is  implicated.  There  is  also  evidence  from  the  selectivity  of  the  resonant  photoe¬ 
mission  enhancement  of  the  valence  band  that  there  is  also  strong  4/  hybridization 
with  the  GaN  valence  bands,  particularly  within  the  Er  4 /5/2  envelope.  The  results 
here  on  the  placement  of  the  occupied  Gd,  Er  and  Yb  4/  levels,  deep  within  the  va¬ 
lence  band,  suggests  that  the  intra-atomic  /-/  transitions  may  be  more  ‘blue’  than 
predicted  by  many  theoretical  models.  The  resonant  photoemission  results  strongly 
support  the  predicted  hybridization  with  the  host  lattice  [55],  and  the  expected  f  —  s 
hybridization  [50,53-55]. 
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VI.  Schottky  Barrier  Formation  at  the  Au  to  Rare  Earth 
Doped  GaN  Thin  Film  Interface 


The  contents  of  this  chapter  were  accepted  for  publication  in  reference  [1].  The 
Schottky  barriers  formed  at  the  interface  between  gold  and  various  rare  earth  doped 
GaN  thin  films  (RE  =  Yb,  Er,  Gd)  were  investigated  in  situ  using  synchrotron 
photoemission  spectroscopy.  The  resultant  Schottky  barrier  heights  were  measured 
as  1.68  ±  0.1  eV  (Yb:GaN),  1.64  ±  0.1  eV  (Er:GaN),  and  1.33  ±  0.1  eV  (Gd:GaN). 
We  find  compelling  evidence  that  thin  layers  of  gold  do  not  wet  and  uniformly  cover 
the  GaN  surface,  even  with  rare  earth  doping  of  the  GaN.  Furthermore,  the  trend  of 
the  Schottky  barrier  heights  follows  the  trend  of  the  rare  earth  metal  work  function. 

6.1  Introduction 

During  the  past  decade,  rare  earth  doped  semiconductors  have  generated  consid¬ 
erable  attention  for  their  application  in  new  optoelectronic  devices  [2-5].  The  favor¬ 
able  thermal,  chemical,  and  electronic  properties  of  wide  band  gap,  Ill-nitride  semi¬ 
conductors  suggest  device  feasibility  using  lanthanide-doped  AIN  and  GaN.  More¬ 
over,  the  tunable  bandgaps  of  these  Ill-nitride  alloys  offer  device  applications  across 
the  visible  spectrum  through  the  ultra-violet  range,  to  include  optically  stimulated 
lasing  [6]  and  p-n  junction  light  emitting  diodes  in  the  red  [7]  using  lanthanide-doped 
AIN  and  GaN,  as  well  as  in  the  blue.  Lastly,  the  production  of  thin  film  electrolu¬ 
minescent  phosphors  with  red,  blue,  and  green  emissions  [8-10]  offers  the  promise  of 
full  color  (white)  light  capability. 

The  large  band  gap  (  ~  3.45  eV)  of  GaN  minimizes  the  effects  of  thermal  or 
visible  (or  longer  wave  length)  light  charge  carrier  generation,  while  alloying  with 
a  rare  earth  nitride  should  decrease  the  band  gap.  As  a  general  rule,  the  rare- 
earth  monopnictides  have  band  gaps  of  0.7  to  1.0  eV  [11],  and  in  many  cases  are 
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suspected  to  be  semimetals,  but  if  alloyed  with  AIN  or  GaN  are  of  considerable 
interest  as  semiconductors.  As  a  device  material,  RErGai_xN  is  unlikely  in  principle 
to  result  in  significant  changes  to  the  barrier  heights  and  the  band  gap  of  GaN,  if 
the  RE  doping  level  is  low.  Although,  if  metal  induced  gaps  states  play  a  significant 
role  [12, 13],  even  a  small  amount  of  rare  earth  could  have  a  significant  effect  on  the 
Schottky  barrier  heights.  It  is  important  to  recognize  that  rare  earth  dopant  induced 
strain,  and  a  bulk  concentration  of  even  a  dilute  amount  of  rare  earth  atoms,  can 
significantly  alter  the  surface  chemistry  and  the  surface  enthalpy  leading  to  a  means 
for  adjusting  Schottky  barrier  heights  that  can  accompany  an  engineering  of  the  GaN 
optical  properties. 

Although  gold  is  generally  considered  unreactive  [14]  complications  abound.  Sur¬ 
face  alloying  can  occur  [15]  and  a  large  range  of  experimentally  measured  Schottky 
barrier  heights  has  been  reported  (0.76  -  1.40  eV)  at  the  Au  to  n-type  GaN  in¬ 
terface  [14-30],  using  photoemission  spectroscopy  (PES),  current- voltage  (I-V)  and 
capacitance  voltage  (C-V)  characteristics,  and  internal  photoemission  [31].  However, 
the  generally  accepted  value  is  about  1.08  eV  [32],  Kurtin,  McGill,  and  Mead  [33] 
suggested  that  the  Schottky  barrier  height  on  GaN  should  depend  directly  on  the 
work  function  or  electronegativity  difference  between  the  metal  electrode  and  GaN. 
Foresi  and  Moustakas  [34]  observed  this  direct  correlation  experimentally,  while  Guo 
et  al.  [35]  and  Mori  et  al.  [36]  observed  only  a  weak  dependence  of  the  Schottky 
barrier  height  on  the  metal  work  function  for  n-type  GaN  and  p-type  GaN,  re¬ 
spectively.  The  1998  review  of  metal-GaN  contact  technology  by  Liu  and  Lau  [20] 
reported  that,  for  a  variety  of  contact  metals  with  both  low  and  high  work  func¬ 
tions,  Schottky  barrier  heights  at  the  metal-GaN  interface  varied  with  metal  work 
function,  within  the  experimental  scatter.  Subsequent  work  by  Rickert  et  al.  [24] 
supported  a  modified  Schottky-Mott  model  at  the  metal-GaN  interface  for  Ni,  Pd, 
and  Al,  yet  more  ‘complex’  behavior  when  Au,  Ti,  and  Pt  were  used  as  the  contact 
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metals.  Additional  experiments  by  Barinov  et  al.  [22,  23]  reported  Schottky  bar¬ 
rier  heights  at  the  Au-GaN  interface  that  exceeded  both  work  function  difference 
(Schottky-Mott)  and  electronegativity  difference  (metal  induced  gap  states)  models. 
Thus,  a  point  worth  re-emphasizing  is  that  regardless  of  the  particular  metal-GaN 
interface  studied,  experimentally  measured  barrier  heights  can  vary  considerably. 

Using  photoemission  to  measure  the  surface  barrier  height  is  advantageous  be¬ 
cause  the  technique  is  both  extremely  surface  sensitive,  and  one  can  avoid  some  of 
the  complications  associated  with  other  experimental  techniques.  For  example,  when 
using  traditional  I-V  and  C-V  measurements,  defects  at  the  metal-semiconductor  in¬ 
terface  can  often  lead  to  overestimates  of  the  surface  barrier  height  [24,37]. 

Although  likely  to  locally  strain  the  lattice,  the  4/  rare  earths  will  tend  to  adopt 
substitutional  sites  for  Ga  [2,  3,  38]  in  GaN  while  significantly  altering  magnetic 
and  optical  properties  [11],  and  it  is  of  considerable  interest  to  know  whether  even 
low  concentrations  of  a  rare  earth  in  the  GaN  host  will  routinely  lead  to  high  or 
low  barrier  height  values  at  the  Au-RE:GaN  interface.  With  this  is  mind,  we  have 
engaged  in  investigations  of  the  interface  properties  of  RExGai_xN  (RE  =  Yb,  Er, 
Gd)  semiconductors  with  Au  metal  overlayer  deposition  under  UHV  conditions.  Our 
studies  of  the  Au  to  RE^Gai.^N  semiconductor  interface  properties  were  performed 
much  in  the  manner  of  other  studies  of  the  Au  to  GaN  interface  [19,21-24], 

6.2  Experimental 

The  experimental  details  were  identical  to  those  discussed  for  resonant  photoe¬ 
mission  in  Section  5.2,  with  the  following  exceptions: 

(1)  The  valence  band  and  Ga  3d  shallow  core  spectra  were  taken  with  a  pho¬ 
ton  energy  hu  =  90  eV.  Au  4/  spectra  were  taken  with  a  photon  energy  of 
hu  =  140  eV. 

(2)  Binding  energies  are  shown  as  absolute  values  |  E  —  Ef  |,  instead  of  negative 
values. 
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(3)  The  Au  deposition  was  made  by  thermal  evaporation  on  the  clean  RE:GaN  sur¬ 
faces  at  room  temperature.  To  prevent  contamination,  the  background  vacuum 
pressure  was  generally  kept  at  <  10~9  torr.  The  evaporation  rate  and  aver¬ 
age  coverage,  reported  here  in  A,  was  monitored  by  means  of  a  quartz  crystal 
thickness  monitor  located  in  the  evaporation  chamber.  A  low  deposition  rate 
of  ~  0.2  A/min  was  used  in  a  best  effort  to  enhance  uniform  film  growth  and 
ensure  accurate  determinations  of  the  average  Au  film  thickness. 

6.3  Results  and  Discussion 

6.3.1  The  Rare  Earth  Doped  GaN  Surfaces. 

The  relative  position  of  the  valence  band,  and  changes  due  to  the  band  bending, 
were  determined  by  monitoring  the  Ga  3d  core  level  shift  as  a  function  of  Au  coverage. 
The  measured  Ga  3d  peaks  in  Figures  34(a),  35(a),  and  35(b)  are  very  similar  to  the 
photoemission  spectra  reported  previously  [21-23,39,40].  The  Ga  3d  shallow  core 
level  peaks  are  broad,  with  maxima  at  17.7  eV  (Yb:GaN),  17.1  eV  (Er:GaN),  and  17.4 
eV  (Gd:GaN)  below  the  valence  band  maximum  (VBM),  in  acceptable  agreement 
with  the  expected  value  of  17.7  eV  [22,23,39]  for  wurtzite  GaN.  The  energy  difference 
between  the  valence  band  maximum  and  a  core  level,  Evbm-c,  is  considered  a  bulk 
property  of  the  material  and  is,  therefore,  independent  of  metal  coverage  [22-24], 
The  valence  band  maximum,  determined  by  extrapolation  of  the  high  kinetic  energy 
edge  of  the  clean  spectrum,  as  is  the  common  practice  [21-24,40,41],  is  at  2.7  eV 
below  Ep  for  the  Yb:GaN  surface.  Assuming  that  the  low  RE-doping  levels  (1-2 
%)  of  our  samples  leaves  the  3.45  eV  band  gap  for  GaN  relatively  unchanged,  and 
given  that  Ep  is  approximately  50  meV  below  the  conduction  band  minimum  in  the 
bulk  of  n-type  GaN  [21],  determines  the  upward  band  bending  at  the  clean  Yb:GaN 
surface  (bare  surface  barrier  height)  as  0.70  eV.  A  consequence  of  the  high  photon 
flux  environment  (>  1011  photons-/im-2-s-1)  required  by  photoelectron  spectroscopy 
are  surface  photovoltage  effects,  which  have  been  known  to  cause  band  flattening  [22], 
By  varying  the  sample  temperature  in  a  separate  experiment,  we  demonstrated  that 
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surface  photovoltage  was  negligible  for  our  samples.  Thus,  we  consider  our  measured 
bare  surface  barrier  height  to  be  in  good  agreement  with  values  reported  previously 
for  UHV  prepared  n-type  wurtize  GaN  surfaces  (0.50-1.40  eV)  [21,22,40-42], 

The  ionization  energy 

IE  =  Ey AC  —  Ey BM  =  (Ep  ~  Ey bm)  +  <t>a  ,  (79) 

given  the  analyzer  work  function  <fia  of  4.4  eV,  is  7.1  ±  0.1  eV  for  Yb:GaN,  which  is 
in  reasonable  agreement  with  prior  values  (6. 7-6. 9  eV)  [21,42]  for  GaN.  The  electron 
affinity 

X  =  4>a  ~  [Eg  —  (EF  —  EVBm)]  ,  (80) 

is  thus  3.7  ±  0.1  eV  for  Yb:GaN,  which  is  in  reasonable  agreement  (2.2-4. 1  eV) 
[21,40,42,43]  with  previously  reported  results  for  GaN.  Table  1  summarizes  the 
experimentally  measured  surface  electronic  properties  of  each  RE:  GaN  sample  and 
indicates  reasonable  agreement,  given  the  low  RE-doping  levels,  between  our  samples 
and  previously  published  GaN  values. 


Table  3.  Measured  properties  of  clean  RE:GaN  surfaces.  Experimental  uncertainties 
are  listed  only  when  explicitly  stated  within  the  indicated  references. 


Material 

EF  —  Ey  BM 

BSBH 

IE 

X 

Yb:GaN 

2.70  ±  0.1 

0.70  ±  0.1 

7.10  ±  0.1 

3.70  ±  0.1 

Er:GaN 

2.80  ±  0.1 

0.60  ±  0.1 

7.20  ±  0.1 

3.80  ±  0.1 

Gd:GaN 

2.65  ±  0.1 

0.75  ±  0.1 

7.05  ±  0.1 

3.65  ±  0.1 

GaN 

GaN 

2.60  [21] 

2.80  1  0.1  [22] 

0.75  ±  0.1  [21] 
0.50  ±  0.1  [22] 

6.90  ±  0.1  [21] 

3.50  ±  0.1  [21] 

GaN 

1.90  ±  0.2  [40] 

1.40  ±  0.2  [40] 

2.70  ±  0.2  [40] 

GaN 

2.70  [41] 

0.70  [41] 

2.20  -  4.10  [43] 

GaN 

2.50  [42] 

0.90  [42] 

6.70  [42] 

3.30  [42] 
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Binding  Energy  |E  -  Ep|  [eV] 

(a)  Yb:GaN 


90  88  86  84  82  80  14  12  10  8  6  4  2  0 


Binding  Energy  |E  -  Ep|  [eV] 

(b)  Au  4f 


Figure  34.  Deconvolution  of  Ga  3 d  core-level  spectra  (a)  and  evolution  of  valence 
band  and  Au  4/  bands  (b)  with  increasing  Au  coverage  on  Yb:GaN  thin  film.  The 
components  attributed  to  bulk  GaN  features  are  shown  with  solid  black  lines.  The 
lower  binding  energy  ‘surface’  component,  shown  as  the  orange  line,  is  removed  with 
increasing  Au  coverage  and  is  replaced  with  Au-GaN  and  Au-RE  alloy  features  (solid 
green  lines). 
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Binding  Energy  |E  -  Ep|  [eV] 

(a)  Er:GaN 


Binding  Energy  |E  -  Ep|  (eV) 

(b)  Gd:GaN 


Figure  35.  Deconvolution  of  Ga  3 d  core-level  spectra  and  valence  band  evolution  with 
increasing  Au  coverage  on  Er:GaN  (a)  and  Gd:GaN  (b)  thin  films.  The  components 
attributed  to  bulk  GaN  features  are  shown  with  solid  black  lines.  The  lower  binding 
energy  ‘surface’  component,  shown  as  the  orange  solid  line,  is  removed  with  increasing 
Au  coverage  and  is  replaced  with  Au-GaN  and  Au-RE  alloy  features  (solid  green  lines). 


Ill 


6.3.2  The  Growth  of  Gold  Overlayers  on  RE:GaN  Surfaces. 

Figure  34  shows  the  evolution  of  the  valence  band  and  Ga  3 d  spectra  for  a  Yb:GaN 
sample  with  increasing  An  coverage.  The  progressively  metallic  nature  of  the  over¬ 
layer  film  is  reflected  in  the  valence  band  features  by  a  density  of  states  above  the 
valence  band  maximum  at  a  Au  overlayer  him  thickness  of  about  6  A,  leading  to 
emission  at  or  near  the  Fermi  level  (Ep)  at  about  12  A  of  An.  The  valence  band 
of  the  gold  overlayer,  as  seen  in  the  photoemission  data,  is  dominated  by  the  An  5 d 
weighted  bands  [44],  as  is  expected,  with  the  An  5 d  levels  at  binding  energies  of  3.5 
and  6.1  eV,  in  reasonable  agreement  with  published  values  for  bulk  Au  [45]  and  gold 
deposited  on  III-V  semiconductors  [46]. 

We  also  measured  the  Au  4/7/2  and  Au  4/5/2  levels  as  a  function  of  Au  coverage, 
as  shown  in  Figure  34(b).  The  Au  4/  features  appear,  at  a  coverage  of  4  A,  at 
binding  energies  that  are  approximately  0.3  eV  greater  that  the  Au  4/7/2  and  Au 
4/5/2  levels  of  84.0  eV  and  87.7  eV  in  metallic  Au  [47].  With  increasing  Au  coverage 
up  to  16  A,  these  peaks  sharpen,  increase  in  intensity,  and  shift  toward  the  lower 
binding  energies  of  bulk  Au.  This  decrease  in  the  Au  4/7/2  and  Au  4 /5/2  levels,  with 
increasing  gold  coverage,  is  one  of  many  indicators  that  initial  gold  adsorption  is  not 
uniformly  wetting  the  surface  and  that  island  growth  is  likely. 

From  photoemission  work  with  Au  alloys,  it  is  known  that  Au  surface  distribution 
may  be  probed  by  following  the  details  of  the  5 d  bands  [48] .  When  Au  atoms  are  well 
dispersed,  as  in  a  dilute  Au  alloy,  5 d  bands  are  shifted  toward  higher  binding  energies, 
as  compared  to  the  bulk  values.  The  splitting  of  the  5 d  valence  band  features  should 
increase  from  the  value  for  atomic  Au  (1.5  eV)  to  the  bulk  value  (2.5  eV)  with 
increased  deposition.  In  our  studies,  we  find  the  Au  5 d  band  splitting  is  about  2.6 
eV  for  Yb:GaN  (Figure  34(b))  and  Gd:GaN  (Figure  35(b)),  almost  independent  of 
Au  overlayer  coverage,  indicative  of  rapid  metallic  Au  island  formation. 
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Attenuation  of  the  Ga  3 d  core  level  in  Figures  34(a),  35(a),  and  35(b)  (RE 
=  Yb,  Er,  Gd)  also  indicates  that  the  growth  mode  of  the  An  on  the  respective 
RE:GaN  surfaces  is  not  uniform  (Frank-van  der  Merwe  or  layer-by-layer  growth)  but 
rather  the  An  overlayer  growth  follows  Volmer- Weber  (island  growth)  or  Stranski- 
Krastanov  growth.  While  the  growth  mode  strongly  depends  on  the  particular  metal- 
semiconductor  pair  and  also  on  the  experimental  conditions  such  as  substrate  tem¬ 
perature  and  evaporation  rate  [49],  the  evidence  for  Volmer- Weber  (island  growth) 
or  Stranski-Krastanov  growth  in  describing  Au  on  the  respective  RE:GaN  surfaces 
is  robust.  If  the  gold  were  of  uniform  thickness,  the  Ga  3d  core  level  photoemission 
signal  would  be  absent  at  thicker  Au  overlayer  coverages  due  to  the  limited  mean 
free  path  of  approximately  5  A  in  Au  [50,51],  estimated  from  the  kinetic  energy  of 
a  Ga  3d  photoelectron  using  a  photon  energy  of  90  eV. 

Since  the  Ga  3d  core  levels  are  visible  in  the  spectra  up  to  16  A  Au  coverage,  thin 
spots  must  exist  in  the  Au  overlayer  him.  The  intensity  of  the  Ga  3d  core  level  before 
I0  and  after  Is  deposition  leads  to  a  change  in  photoemission  intensity,  described  by 
Volmer- Weber  growth  [24,49,52,53],  as 

y-  =  (1  -  0)  +  6e~tlx  ,  (81) 

where  t  is  the  him  thickness,  A  is  the  mean  free  path  of  the  electrons,  and  9  represents 
the  fractional  surface  coverage  reached  prior  to  island  growth  in  three  dimensions. 
Figure  36  shows  the  photoemission  intensity  ratios  of  the  Ga  3d  core  level  for  each 
of  the  RE:GaN  thin  films  studied.  As  expected,  the  ht  parameter  differs  for  each 
Au  to  RE:GaN  interface  studied,  which  rehects  the  strong  growth  mode  dependence 
upon  the  particular  metal-semiconductor  pair. 

The  data  indicate  that  Au  is  not  growing  in  a  layer-by-layer  manner  but  rather 
by  Volmer- Weber  (island  growth)  or  Stranski-Krastanov  growth.  Although  layer-by- 
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layer  growth  for  Au  on  GaN  has  been  reported  [19],  there  has  also  been  evidence  of 
island  formation  [19,24,54].  Here,  we  find  compelling  evidence  that  thin  layers  of 
gold  do  not  wet  and  cover  the  GaN  surface,  even  with  rare  earth  doping  of  the  GaN. 
We  would  expect  that  a  dilute  surface  coverage  of  rare  earth  atoms  would  nucleate 
more  uniform  Au  overlayer  films,  but  this  does  not  appear  to  be  the  case. 


Figure  36.  The  intensity  ratio  of  the  Ga  3 d  core  level  before  I0  and  after  Is  Au 
deposition  as  a  function  of  film  thickness  for  the  deposition  of  Au  on  RE:GaN,  as 
indicated.  The  expected  ratio  as  a  function  of  thickness  (dashed  lines)  is  plotted  using 
the  Volmer- Weber  growth  mode. 


6.3.3  Schottky  Barrier  Formation  at  RE:GaN  Surfaces. 

When  metal  is  evaporated  on  the  sample  surface,  the  edge  of  the  semiconductor 
valence  band  maximum  is  obscured  by  the  signal  originating  from  the  metal  overlayer. 
However,  owing  to  the  fact  that  AVbm-C  is  constant,  the  Ga  3d  core  shift  is  roughly 
equivalent  to  the  valence  band  shift  at  the  semiconductor  surface,  from  which  the 
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Schottky  barrier  height  $B,n  is  calculated  as 


*&B,n  —  Eg  —  (Ep  —  -Evbm)  •  (82) 

To  describe  the  surface  Fermi  level  movement  and  Schottky  barrier  formation  during 
Au  deposition,  Ga  3d  core  levels  were  deconvoluted  into  surface  and  bulk  components 
with  Gaussian  form.  There  exists  some  variation  concerning  the  number  of  fitting 
components  for  the  Ga  3d  lineshape.  Some  authors  choose  a  single,  dominant  bulk 
subpeak  and  one  surface  subpeak  [39,55]  to  represent  the  surface  to  bulk  core  level 
shift  that  affects  the  Ga  3d  lineshape,  whereas  other  researchers  deconvolute  the 
lineshape  using  two  bulk  subpeaks  and  one  surface  subpeak  [22,23,56].  We  calculated 
Schottky  barrier  heights  using  both  methods,  and  their  differences  proved  negligible 
(<  0.05  eV)  for  all  Au  to  RE:GaN  interfaces.  However,  we  present  only  our  results 
using  the  latter  method  to  deconvolute  the  Ga  3d  peak,  as  this  is  deemed  more 
reliable.  We  selected  this  method  owing  to  the  compelling  arguments  of  Barinov  et 
al.  [22, 23]  and  Lambrecht  et  al.  [56]  who  demonstrated  that  the  dominant  and  high 
binding  energy  components  of  the  peak  are  intrinsic  features  of  the  Ga  3d  semi-core 
levels  of  GaN  and  are  related  to  hybridization  effects  [20],  while  the  lower  binding 
energy  component  behaves  as  a  ‘surface’  component.  In  Figures  34(a),  35(a),  and 
35(b),  black  subpeaks  represent  bulk  components  and  orange  subpeaks  represent  the 
surface  components.  As  Au  coverage  increases,  Au  forms  a  surface  alloy  with  GaN 
and  the  RE  metals  on  the  sample  surface,  and  their  features  are  represented  as  green 
subpeaks.  These  features  supplant  the  surface  component  of  the  Ga  3d  core  levels. 

The  effect  of  Au  deposition  on  the  various  RE  doped  GaN  surfaces  is  not  identical, 
as  inferred  from  the  Ga  3d  shallow-core  spectra.  The  clean  Yb:GaN  Ga  3d  shallow- 
core  spectra  in  Figure  34(a)  and  the  clean  Gd:GaN  Ga  3d  shallow-core  spectra  in 
Figure  35(b)  required  three  fitting  components  to  represent  the  Ga  3d  lineshape. 
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During  Au  deposition  on  Yb:GaN,  the  single  surface  component  was  replaced  with 
two  alloy/Au  interface  components  that  were  separated,  at  lower  binding  energies, 
from  the  dominant  bulk  peak  by  1.3  eV  and  2.9  eV.  During  Au  deposition  on  Gd:GaN, 
the  single  surface  component  was  replaced  with  a  single  alloy/Au  interface  component 
that  was  separated,  at  a  lower  binding  energy,  from  the  dominant  bulk  peak  by  0.7 
eV. 

The  clean  Er:GaN  spectrum  in  Figure  35(a)  required  four  subpeaks  to  deconvolute 
the  Ga  3d  lincshape.  We  associate  the  fourth  subpeak,  at  binding  energy  24.8  eV, 
with  the  Er  5 p  core  level,  which  has  a  binding  energy  of  28.0  eV  in  bulk  Er.  During 
Au  deposition,  the  single  surface  component  was  replaced  with  one  alloy  component 
that  was  separated,  at  a  lower  binding  energy,  from  the  dominant  Ga  3d  bulk  peak 
by  0.6  eV.  Additional  fitting  components  were  needed  to  deconvolute  the  spectral 
features  that  evolved  around  the  Er  5 p  subpeak.  One  alloy  component  located, 
at  a  higher  binding  energy,  from  the  Er  5 p  subpeak  by  3.0  eV  is  present  in  all 
spectra  from  4  A-  16  A  Au  coverage.  The  intensity  of  this  Ga  3d  shallow-core 
spectra  component  increases  from  4  A-  12  A  Au  coverage,  but  is  attenuated  at 
16  A  Au  coverage,  possibly  the  result  of  surface  alloy  formation.  When  this  Ga 
3d  shallow-core  spectra  component  is  at  maximum  intensity  at  12  A  Au  coverage, 
precise  spectral  deconvolution  requires  two  additional  alloy  components.  These  two 
additional  alloy/Au  interface  components  are  separated  from  the  Er  5 p  subpeak  by 
lower  and  higher  binding  energies  of  1.7  eV  and  5.3  eV,  respectively. 

Fermi  level  movement,  surface  valence  band  bending,  and  Schottky  barrier  forma¬ 
tion  were  determined  by  monitoring  the  binding  energy  shift  of  the  dominant,  bulk 
Ga  3d  component  during  Au  deposition.  Figures  34(a),  35(a),  and  35(b)  show  these 
measured  shifts  to  be  0.98  eV  (Yb:GaN),  1.04  eV  (Er:GaN),  and  0.58  eV  (Gd:GaN). 
Figure  37  shows  the  surface  Fermi  level  movement  37(a)  and  resultant  Schottky  bar¬ 
rier  formation  37(b)  for  Gd:GaN,  which  was  calculated  as  1.33  ±0.1  eV  via  (82). 
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The  same  methods  and  calculations  described  above  yielded  Schottky  barrier  heights 
of  1.64  ±  0.1  eV  and  1.68  ±  0.1  eV  for  Er:GaN  and  Yb:GaN,  respectively.  These  are 
generally  far  larger  than  the  values  found  for  undoped  GaN,  as  measured  via  PES, 
I-V,  C-V,  and  internal  photoemission  (IPE),  and  summarized  in  Table  4. 


Table  4.  Comparison  between  experimentally  measured  and  theoretically  predicted 
(Schottky-Mott)  barrier  heights.  Experimental  uncertainties  are  listed  only  when 
explicitly  stated  within  the  indicated  references. 


Material 

Xs  [eV] 
via  Table  3 

$b,  „  [eV] 
(theory) 

$b,  „  [eV] 
(measured) 

Measurement 

technique 

Yb:GaN 

3.70  ±  0.1 

1.45  ±  0.1 

1.68  ±  0.1 

PES 

Er:GaN 

3.80  ±  0.1 

1.35  ±  0.1 

1.64  ±  0.1 

PES 

Gd:GaN 

3.65  ±  0.1 

1.50  ±  0.1 

1.33  ±  0.1 

PES 

GaN 

4.10  [43] 

1.05 

1.15  [19] 

PES 

GaN 

1.20  ±0.1  [21] 

PES 

GaN 

1.40  ±0.1  [22,23] 

PES 

GaN 

0.90  ±0.1  [24] 

PES 

GaN 

1.18  ±0.07  [26] 

I-V 

GaN 

0.94  [16] 

I-V 

GaN 

1.19  [17] 

I-V 

GaN 

1.15  [17] 

C-V 

GaN 

0.87  [18] 

I-V 

GaN 

0.98  [18] 

C-V 

GaN 

1.10  [14] 

I-V 

GaN 

1.22  [15] 

I-V 

GaN 

0.81  [29] 

I-V 

GaN 

0.76  [29] 

IPE 

GaN 

1.0  [30] 

IPE 

GaN 

1.11  [25] 

IPE 

GaN 

0.95  ±0.04  [27] 

IPE 

GaN 

0.97  ±0.05  [28] 

IPE 

Because  of  the  possible  interfacial  reactions  between  the  Au  and  RE: GaN  surfaces, 
the  magnitude  of  the  shift  in  the  Ga  3d  shallow-core  spectra  due  to  Au  charge 
donation  cannot  be  determined  quantitatively.  While  we  associate  the  Ga  3d  core 
shift  with  band  bending,  we  must  be  careful  to  state  that  other  causes  cannot  be 
completely  excluded  by  the  data  presented  here. 
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(a)  Gd:GaN  Fermi  level  movement 


0 

Au  Thickness  (A) 

(b)  Gd:GaN  Schottky  barrier  formation 


Figure  37.  Fermi  level  movement  (a)  and  Schottky  Barrier  formation  (b)  at  the  surface 
of  Gd:GaN  based  upon  the  binding  energy  of  the  bulk  component  of  a  Ga  3 d  core  level 
taken  at  a  photon  energy  of  90  eV.  The  filled  region  in  (b)  indicates  the  range  of  ideal 
barrier  height  values,  as  predicted  by  the  Schottky-Mott  model  using  our  measured  x 
values  from  Table  3. 
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The  Schottky-Mott  relationship  characterizes  contacts  for  n-type  semiconclnctors 
such  that  when  the  work  function  of  the  metal  contact  is  greater  than  the  work 
function  of  the  semiconductor,  i.e.  (f)m  >  c j)s ,  the  contact  is  rectifying.  The  work 
function  0S  for  GaN  has  been  calculated  from  photoemission  experiments  as  4.2  ± 
0.2  eV  [41],  whereas  the  work  function  0m  of  Au  is  5.15  eV  [57].  The  Schottky- 
Mott  theory  also  predicts  that  for  n-type  semiconductors,  with  (f)m  >  <j)s,  the  surface 
barrier  height  is  calculated  as  the  difference  between  the  contact  metal  work  function 
and  the  electron  affinity  of  the  semiconductor  sample  Xs  as 

*&B,n  =  0m  Xs  ■  (83) 

The  most  commonly  used  [18,34,39]  electron  affinity  for  undoped  GaN  was  reported 
as  4.1  eV  [43],  which  yields  a  theoretical  barrier  height  of  1.05  eV  between  Au 
and  GaN.  Our  measured  values  for  yy,  determined  from  the  clean  RE: GaN  bare 
surface  barrier  heights  (Table  3),  suggest  higher  theoretical  Schottky  barrier  heights 
(1.45  ±  0.1  eV  (Yb:GaN),  1.35  ±  0.1  eV  (Er:GaN),  1.50  ±  0.1  eV  (Gd:GaN)  than 
for  undoped  GaN  (Table  2).  Thus,  both  measured  and  theoretically  predicted  Au- 
RE:GaN  Schottky  barrier  heights  are  larger  than  those  for  the  Au-GaN  interface. 
Given  a  constant  electrode  metal  work  function,  either  a  decreased  semiconductor 
electron  affinity  (83)  or  decreased  semiconductor  work  function  (increased  contact 
potential)  due  to  rare  earth  doping  of  GaN  would  produce  higher  barrier  heights 
than  for  undoped  GaN,  but  more  likely,  both  occur. 

We  note  that,  in  general,  the  trend  of  the  Schottky  barrier  heights  follows  the 
trend  of  the  rare  earth  metal  work  function,  so  that  although  the  RE  ion  occupies  a 
Ga  site  in  the  GaN  [2,3,38],  perturbation  of  the  rare  earth  on  the  surface  electronic 
structure  of  GaN  is  possible.  The  work  function  of  Yb  is  2.60  eV  [58],  that  of  Er  is  2.97 
eV  [59],  while  Gd  is  3.10  eV  [57],  all  of  which  are  lower  than  that  of  GaN  (4.2  eV). 
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Regardless  of  the  interplay  between  electron  affinity  and  GaN  surface  work  func¬ 
tion,  the  resultant  Au-RE:GaN  Schottky  barrier  heights,  as  measured  and  noted 
above,  are  1.68  ±  0.1  eV  (Yb:GaN),  1.64  ±  0.1  eV  (Er:GaN),  and  1.33  ±  0.1  eV 
(Gd:GaN),  and  are  significantly  higher  than  observed  for  undoped  GaN.  Thus,  band 
bending  at  the  Au-RE:GaN  (RE  =  Yb,  Er,  Gd)  interface  is  larger  than  the  nor¬ 
mal  (Au-GaN).  This  implies  that  barrier  heights  might  be  engineered  to  optimize 
depletion  widths  and  charge  collection  volumes  for  sensor  device  applications. 

6.4  Conclusions 

Photoemission  studies  using  synchrotron  radiation  showed  that  the  Schottky  bar¬ 
rier  heights  between  Au  and  RE:GaN  thin  film  samples  were  measured  to  be  1.68 
±  0.1  eV  (Yb:GaN),  1.64  ±  0.1  eV  (Er:GaN),  and  1.33  ±  0.1  eV  (Gd:GaN).  This 
trend  of  the  Schottky  barrier  heights  follows  the  trend  of  the  rare  earth  metal  work 
function.  The  Au  overlayer  does  not  wet  and  cover  the  GaN  surface,  even  with  rare 
earth  doping  of  the  GaN.  But  in  spite  of  the  imperfections  of  the  Au-RE:GaN  inter¬ 
face,  the  resulting  Schottky  barrier  interfaces  might  lead  to  significant  improvement 
in  device  performance  in  sensor  applications. 
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VII.  Neutron  Detection  Using  Gadolinium 


7.1  Overview 

Owing  to  the  fact  that  neutron  detection  offers  the  greatest  likelihood  to  identify 
the  presence  of  special  nuclear  material,  this  chapter  will  examine  the  rare  earth  metal 
gadolinium  in  that  capacity.  The  thermal  neutron  capture  cross  section  of  15 'Gd 
(255,000  barns)  is  the  largest  of  all  known  natural  isotopes,  which  distinguishes  the 
material  as  a  logical  candidate  for  neutron  detection.  Recall  that  neutron  detection 
requires  the  conversion  of  neutrons  to  ionizing  radiation  that  generates  electron-hole 
pairs  via  energy  deposition  in  the  active  collection  volume  of  a  device.  Ignoring 
contact  types  for  the  following  discussions,  the  charge  produced  from  a  single  15 'Gd 
neutron  capture  event  will  be  analyzed  and  compared  to  the  intrinsic  noise  level  of 
typical  commercially  available  charge-to-voltage  preamplifiers.  Next,  material  and 
geometry  considerations  will  be  discussed  and  applied  to  multiple  capture  events 
due  to  a  hypothetical  neutron  flux.  The  single  and  multiple  event  calculations  will 
be  compared  to  boron  carbide  (B5C),  which  is  an  emerging  material  in  the  held  of 
solid  state  neutron  detectors.  Although  the  calculations  presented  in  the  subsequent 
sections  are  somewhat  idealized,  they  serve  to  identify  those  parameters  of  device 
design  that  might  be  suited  for  further  analysis  via  detailed  computational  model¬ 
ing  prior  to  device  fabrication.  Additionally,  the  fundamental  calculations  serve  to 
answer  overarching  questions  associated  with  an  engineering  study  of  this  type,  to 
include: 

(1)  Can  a  157Gd  neutron  capture  event  generate  a  signal  that  is  detectable  above 
the  intrinsic  noise  of  the  pulse  processing  electronics  and  the  leakage  current 
of  a  typical  diode? 

(2)  What  is  the  anticipated  frequency  response  of  the  diode  in  the  presence  of  a 
thermal  neutron  flux? 
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(3)  What  are  the  ranges  and  mean  free  paths  in  a  Gd-based  semiconductor  for  the 
charged  particles  and  photons  created  in  a  capture  event,  and  how  do  those 
distances  compare  with  the  depletion  width  of  the  diode? 

(4)  How  do  Gd-based  devices  perform  against  the  boron  carbide  benchmarks  for 
single  and  multiple  capture  events? 


7.2  Single  Neutron  Capture  Event 

7.2.1  Gd-157  Neutron  Capture  Scheme. 


19.1%  of  de-excitation  events 
lead  to  the  first  three  levels  of 
the  rotational  band  K  =  0 


0  keV 


b-  261  keV 
J-79keV_J 

a  «  30% 


y^+  7676.0  keV  (4.6%) 


6750.0 , 925.9  (5.55%) 

6757.9 , 918.2  (3.69%) 

2939.9 . 4736.1  (3.46% ) 

5659.8  ,  2016.2  (334%  ) 

....plus  22  more  with  abs.  int.  >  0.03% 


y  +  y_  =  7857.9  keV  (11.3%) 

67503  . 1107.6  (13.6%  ) 

5177.8  ,  2680.0  (2.61%) 

5676.8 , 2181.0  (2.34%) 

6672.3  ,  1185.6  (1.95%) 

_ plus  66  more  with  abs.  int.  >  0.03% 


•Internal  Conversion  Electrons 
182  keV  ) 

79keV  I  96.7%  of*  Energy 


•X-Rays 

~43  IteV,  K-Shell  1  94.2%of 
~7  keV,  L-Shell  j  X-ray  Energy 

•Coster-Kronig  Auger  Electrons 
~1  keV,  L-Shell 


y^+  y  =  7037 A  keV  (3.2%) 


6750.0  . 1187.2  (24.6% ) 

6420.3  .  15173  (3.15%) 

5972.5  ,  19650  (2.68%) 

5678.0 , 2259.4  (2.18%) 

— plus  38  more  with  abs.  int.  >  003% 


Figure  38.  De-excitation  schemes  of  158Gd*  following  neutron  capture  by  157Gd. 
Adapted  from  [1—3]. 


To  determine  the  amount  of  energy  available  for  electron-hole  pair  production 
in  a  Gd-based  semiconductor  device,  the  157Gd  neutron  capture  and  subsequent  de- 
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excitation  reactions  must  be  examined.  Neutron  capture  by  157Gd  results  in  an 
excited,  or  compound  158Gd*  nucleus  that  is  7937  keV  from  the  ground  state  energy 
of  158Gd.  The  compound  nucleus  can  decay  to  the  ground  state  by  several  paths, 
and  19.1%  of  these  de-excitations  are  well  understood  and  summarized  in  Figure  38. 

Decay  path  C  in  the  figure  is  the  most  straightforward,  occurs  3.2%  of  the  time, 
and  results  in  coincident  emission  of  two  gamma  rays  (71,72)  whose  energy  sum 
equals  7937  keV,  as  expected.  The  distribution  of  the  gamma  ray  energies  and  their 
probabilities  of  occurrence  are  listed  below  the  gamma  ray  energy  sum  at  the  lower 
left  corner  of  the  figure.  Decay  paths  A  and  B  occur  4.6%  and  11.3%  of  the  time, 
respectively,  and  result  in  a  158Gd*  nucleus  that  is  261  keV  and  79  keV,  respectively, 
from  the  ground  state  energy  of  158Gd.  The  decays  to  those  energy  levels  produce 
gamma  rays  whose  energy  sum  equals  (7937  keV  -  261  keV)  and  (7937  keV  -  79  keV) 
for  paths  A  and  B,  respectively,  and  their  probabilities  of  occurrence  are  listed  below 
their  gamma  ray  energy  sums.  Following  these  initial  decays,  and  shown  as  event  D 
in  the  figure,  the  158Gd*  nucleus  may  decay  to  the  ground  state  energy  of  158Gd  via 
emission  of  a  182  keV  or  79  keV  gamma  ray,  or  through  the  competing  process  of 
internal  conversion  electron  emission. 

Internal  conversion  (IC)  occurs  when  the  excited  nucleus  transfers  energy  to  an 
atomic  electron,  typically  in  the  K  (n  =  1)  or  L  (n  =  2)  electronic  shell,  where  n  is 
the  principle  quantum  number.  The  preference  for  the  competing  processes  of  IC 
and  gamma  emission  is  described  by  the  internal  conversion  coefficient  a  as 


Ne 

a  AT7 

where  Ne  and  N 7  are  the  number  of  internal  conversion  electrons  and  gamma  photons, 
respectively,  for  the  nuclear  transition.  This  coefficient  increases  with  the  cube  of  the 
atomic  number  and  decreases  with  E*\  thus,  internal  conversion  is  prevalent  in  the 
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decay  of  low-lying  excited  states  (small  E*)  of  heavy  nuclei  [4],  Lastly,  the  kinetic 
energy  Ee  of  the  ejected  atomic  electron  is  very  nearly  equal  to  the  difference  between 
the  excitation  energy  of  the  E*  nucleus  and  the  binding  energy  of  the  electronic  shell 
Eb  as 

Ee  =  E*-EB.  (85) 

Returning  to  Figure  38,  the  internal  conversion  coefficient  for  the  79  keV  transi¬ 
tion  is  expected  to  be  0.3.  This  predicts  that  30%  of  11.3%  of  events  from  transition 
B  — >  C  will  result  in  IC  emission,  or  3.4%  of  all  157Gd  neutron  capture  events.  It 
follows  that  30%  of  4.6%,  or  1.4%  of  events  from  transition  A  — >  B  (A E  =  182  keV) 
will  result  in  IC  emission.  As  IC  for  Gd  is  known  to  be  preferential  to  L-shcll  electron 
emission  [5],  and  the  average  L-shell  binding  energy  for  Gd  bas  been  theoretically 
calculated  and  experimentally  verified  to  be  approximately  7  keV  [6],  equation  (85) 
predicts  that  72  keV  electron  emission  will  dominate  capture  events  that  result  in  IC. 
The  vacancy  following  ejection  of  an  IC  electron  from  a  low-lying  electronic  orbital 
will  be  filled  by  an  electron  from  a  higher  orbital,  which  will  cause  the  higher  orbital’s 
vacancy  to  be  filled  by  an  electron  from  another  outer  orbital,  and  the  downward 
cascade  will  continue.  The  corresponding  energy  difference  between  electronic  or¬ 
bitals  results  in  a  series  of  characteristic  X-rays  of  approximately  43  keV  and  7  keV 
for  L  — >  K  and  M  — >  L  transitions,  respectively  [6].  The  energy  difference  may  also 
be  transferred  to  an  L-shell  electron,  in  which  case  an  Auger  electron  will  be  emitted 
with  a  kinetic  energy  equal  to  the  difference  between  the  transition  and  the  binding 
energy  of  the  electronic  state  from  which  it  originated,  or  approximately  35  keV. 

If  the  detailed  examination  of  the  157Gd(n,  7)158Gd  reaction  above  appeared  te¬ 
dious,  full  understanding  of  the  ensuing  reactions  facilitates  detector  design  dimen¬ 
sions  and  detection  equipment  settings  to  ensure  prudent  discrimination  when  analyz¬ 
ing  spectral  features.  The  salient  points  of  the  analysis  are  summarized  as  follows: 
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(1)  less  that  5%  (3.4%  +  1.4%)  of  all  157Gd  neutron  capture  events  produce 
charged  particles  (internal  conversion  and  Auger  electrons)  capable  of  gen¬ 
erating  electron-hole  pairs  in  a  Gd-based  semiconductor  device, 

(2)  relatively  low  energy  (keV-range)  photon  production  occurs  in 

(a)  decay  chain  D  via  emission  of  43  keV  77-shell  and  7  keV  L-shcll  charac¬ 
teristic  X-rays,  following  internal  conversion  electron  emission, 

(b)  decay  chain  D  where  182  keV  and  79  keV  gamma  rays  are  emitted  for 
3.2%  (70%  of  4.6%)  and  7.9%  (70%  of  11.3%),  respectively,  of  transitions, 

(3)  the  dominant  high  energy  photon  (MeV-range)  production  occurs  in  decay 
chain  B,  where  13.6%  of  11.3%,  or  1.5%  of  transitions  produce  coincident 
gamma  rays  of  6.7  MeV  and  1.1  MeV. 

Alternatively,  the  significant  reactions  discussed  above,  which  account  for  nearly 
all  of  the  158Gd*  de-excitations,  may  be  tabulated  as  shown  in  Table  5. 


Table  5.  Summary  of  ionizing  radiation  produced  from  neutron  capture  by  157Gd. 


Absolute 
Intensity  [%] 

Radiation 

Produced 

Energy(ies) 

Decay 

Chains 

3.4 

IC  electron 

72  keV 

B,D 

77-shell  X-ray 

43  keV 

L-shell  X-ray 

7  keV 

1.4 

IC  electron 

174  keV 

A,D 

77-shell  X-ray 

43  keV 

L-shell  X-ray 

7  keV 

7.9 

7-ray 

79  keV 

B 

3.2 

7-ray 

182  keV 

A 

1.5 

7-ray 

1.1,  6.7  MeV 

B 

KLL  Auger  electron 

35  keV 

LM M  Auger  electron 

6  keV 

7.2.2  Ionizing  Radiation  Charge  Production. 

The  quantity  of  electron-hole  pairs  produced  by  ionizing  energy  loss  events  is 
dominated  by  the  bandgap  of  the  semiconductor  in  which  they  were  created.  Con¬ 
sider  a  single  157Gd  neutron  capture  event  that  results  in  a  72  keV  internal  conversion 
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electron.  Assuming  a  cost  of  approximately  3 Eg  per  electron-hole  pair  [7],  and  that 
the  internal  conversion  electron  is  doing  most  of  the  work  Eion,  the  maximum  charge 
Qmax  that  can  be  created  from  a  single  event  is  given  as 


Qmax  =  X  (1.6  X  10-19)  [C]  .  (86) 

The  bandgap  for  Gd^Gai-^N  is  assumed  to  vary  linearly  as  a  function  of  molar 
concentration  x  between  GdN  (1.0  eV)  and  GaN  (3.4  eV)  as 

Eg  [Gd^Gai-zN]  =  xEg  [GdN]  +  (1  -  x)  Eg  [GaN] .  (87) 

Charge  generation  dependence  upon  semiconductor  bandgap  for  a  single  157Gd  neu¬ 
tron  capture  event,  via  (86)  and  (87),  is  depicted  in  Figure  39  for  GdxGai_iTN.  The 


Figure  39.  Qmax  creation  from  (red):  a  single  157Gd  neutron  capture  event  resulting 
in  79  keV  internal  conversion  electron  emission,  as  a  function  of  molar  concentration  x 
in  GdxGai_xN;  (black):  a  single  neutron  capture  event  in  B5C  resulting  in  0.84  MeV 
7Li  ion  production,  as  shown  in  (88). 


figure  indicates  that  the  charge  generated  varies  between  1.2  fC  for  GaN  and  4.2  fC 
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for  GdN.  As  a  benchmark  for  comparison,  the  boron  conversion  reactions  in  a  B5C 
device  were  considered 


94%  :  10B  +  11  — >  7Li  +  4He  +  7  ,  (88) 

(0.84)  (1.47)  (0.48) 

6%  :  10B  +  11  — >  7Li  +  4He  ,  (89) 

(1.02)  (1.77) 

where  the  branching  ratios  and  reaction  product  energies  (MeV)  are  shown  to  the 
left  and  beneath  the  reactions,  respectively.  Using  the  most  conservative  charged 
particle  energy  (0.84  MeV)  and  assuming  a  B5C  bandgap  of  0.9  eV  [8],  Figure  39 
shows  that  the  charge  generated  from  a  single  neutron  capture  event  is  more  than 
an  order  of  magnitude  greater  in  B5C  than  in  GdN. 

7.3  Detector  Circuit  Considerations 

The  charge  produced  from  the  interaction  of  a  single  particle  must  be  collected 
to  form  a  basic  electrical  signal.  This  collection  is  accomplished  typically  by  the 
imposition  of  an  electric  field  within  the  detector.  The  collection  time  tc  can  be 
approximated  as 

tc=~.  (90) 

v 

Here,  d  is  the  charge  collection  distance  and  v  is  the  drift  velocity  expressed  as  a 
function  of  charge  mobility  fi  [cm2  V”1  sec"-1]  and  the  applied  electric  field  £  as 


v  =  fi£ . 


(91) 
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Therefore,  the  transient  current  It  produced  by  a  single  event  radiation-induced 
charge  pulse  is  expressed  as  the  ratio  of  equations  (86)  and  (90)  as 

qEion\  f  /i£\ 

3  EgJ\dJ- 

Assuming  a  typical  diode  leakage  current  on  the  order  of  10'9  A,  a  rough  comparison 
between  the  transient  and  leakage  currents  can  be  made  for  a  Gd-based  device.  Using 
H  =  1  cm2  V'1  sec'1,  S  =  105  V  cm'1,  and  d  corresponding  to  a  depletion  width  of 
5 /irn  as  reasonable  order  of  magnitude  estimates,  yields  It  ~  0.2/xA,  which  is  two 
orders  of  magnitude  above  the  leakage  current.  As  a  first  approximation,  and  under 
ideal  conditions,  the  transient  current  is,  as  a  minimum,  marginally  better  than 
leakage. 

For  a  prototype  detector,  the  current  pulse  may  be  illustrated  as  shown  in  Figure 
40,  and  it  follows  from  the  figure  that  the  time-integrated  current  over  the  duration 


c 


Figure  40.  Example  of  detector  current  pulse  where  tc  represents  collection  time. 
Adapted  from  [5]. 
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of  the  pulse  collection  time  must  equal  the  collected  charge  Q,  or 

I(t)  dt  —  Q  .  (93) 

Whereas  Figure  40  shows  a  symmetric  pulse,  in  practice,  many  pulses  will  be  asym¬ 
metric,  owing  to  time  dependent  current  fluctuations,  ft  is  also  assumed  that  either 
the  radiation  interaction  rate  is  low  enough,  or  the  time  constant  of  the  external 
circuit  is  small  compared  to  the  charge  collection  time,  so  that  each  individual  in¬ 
teraction  gives  rise  to  a  distinguishable  current.  In  most  detection  applications,  the 
time  integral  of  each  transient  current,  or  total  charge  Q,  is  recorded,  and  the  de¬ 
posited  energy  in  the  detector  is  related  directly  to  Q.  In  applications  designed  to 
measure  the  energy  of  individual  events,  as  would  be  required  to  identify  an  individ¬ 
ual  internal  conversion  electron  following  neutron  capture  by  15 'Gd,  the  pulse  mode 
of  detector  operation  is  best  [5].  As  an  example,  the  pulse  height  spectra  for  a  B5C 
hetero junction  neutron  detector  in  the  presence  of  various  reactor  fluxes  [9]  is  shown 
in  Figure  41  and  clearly  discerns  the  dominant  spectral  features  predicted  by  the 
boron  neutron  capture  reactions  (88)  and  (89). 

In  practice,  generating  pulse  height  spectra  for  a  solid  state  neutron  detector 
as  shown  in  Figure  41  requires  the  detector  to  be  connected  to  pulse  processing 
electronic  equipment.  This  requirement  is  levied  owing  to  the  fact  that  the  generated 
charge  is  so  small  that  signal  pulses  require  an  intermediate  amplification  step.  The 
elements  of  a  typical  pulse  height  spectrometry  signal  chain  are  shown  in  Figure  42. 
The  voltage  supply  imposes  the  applied  electric  field  in  the  detector,  such  that  the 
radiation  induced  charge  may  be  collected,  converted  to  a  voltage,  amplified,  and 
pulse-shaped  after  departing  the  preamplifier-amplifier  configuration.  The  significant 
contribution  to  noise  occurs  at  the  beginning  of  the  signal  chain,  when  the  signal  is 
minimized.  Any  noise  contribution  at  this  point  undergoes  the  same  amplification  as 
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Figure  41.  Pulse  height  spectra  and  background  (inset)  from  a  B5C  heterojunction 
detector.  Each  spectra  was  collected  for  600  seconds  at  the  reactor  power  indicated. 
The  three  identifiable  peaks  are  attributable  to  the  boron  capture  products  shown  in 
reactions  (88)  and  (89).  The  1.02  MeV  peak  is  lost  in  the  continuum  between  the  0.84 
MeV  and  1.47  MeV  peaks  owing  to  its  branching  ratio  [9]. 


Figure  42.  Elements  of  a  typical  signal  chain  for  pulse  height  analysis.  Adapted 
from  [5]. 


Figure  43.  System  for  measuring  charge-sensitive  preamplifier  noise. 
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the  signal,  whereas  noise  generated  during  later  stages  in  the  signal  chain  is  typically 
much  smaller  than  the  signal.  It  is  therefore  of  significant  interest  to  know  whether 
the  charge  generated  by  a  157Gd  neutron  capture  event  exceeds  the  intrinsic  noise 
level  of  the  preamplifier  input  stage. 

The  intrinsic  noise  level  of  a  charge-sensitive  commercial  preamplifier,  typically 
expressed  as  a  function  of  input  capacitance,  is  determined  by  charge  injection  into 
the  preamplifier,  as  shown  in  Figure  43.  The  charge  can  be  injected  by  a  detector 
and  radiation  source  or  a  step  pulse  generator  connected  to  the  preamplifier  input 
through  a  capacitor,  sometimes  referred  to  as  a  charge  terminator.  The  preamplifier 
noise  can  be  determined  by  measuring  the  root-mean-square  (rms)  noise  voltage  Id-ms 
at  the  output  of  the  filter  amplifier  in  the  absence  of  any  pulses.  The  ratio  Vvms/Vp, 
referred  to  as  the  equivalent  noise  charge  [5],  and  reported  in  either  eV  or  number 
of  electrons,  is  used  to  determine  noise  as 

FWHM  (noise)  =  2.35  e  ^  [eV] ,  (94) 

Vp 

where  e  is  the  deposited  energy  required  to  create  one  charge  carrier.  Typical 
capacitance-dependent  performance  values  are  based  upon  the  silicon  equivalent 
(e  =  3.6  eV)  and  range  between  2-27  keV  for  commonly  used  ORTEC  systems  [10]. 
Conversion  to  the  corresponding  number  of  electrons  is  done  by  dividing  the  reported 
values  by  the  silicon  equivalent,  which  yields  a  range  of  555-7500  electrons,  or  an 
equivalent  noise  charge  of  approximately  0. 1-1.2  fC.  Keeping  in  mind  that  manufac¬ 
turer  performance  specifications  are  determined  under  ideal  conditions,  the  charge 
generated  from  a  15'Gd  neutron  capture  event,  as  shown  in  Figure  39,  approaches  the 
intrinsic  noise  level  of  many  commercially  available  charge-to-voltage  preamplifiers. 

On  the  basis  of  charge  collection,  the  governing  physics  implies  marginal  perfor¬ 
mance  of  Gd-based  semiconductor  neutron  detectors.  However,  if  the  preamplifier 
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gain  can  be  minimized  to  account  for  the  low  charge  production,  gadolinium  may 
warrant  consideration  as  a  detection  material.  Assuming  the  pulse  processing  elec¬ 
tronics  can  be  optimized  as  stated,  the  charge  collection  analysis  conducted  above 
reveals  nothing  about  the  anticipated  count  rate  or  collection  time  required  to  ob¬ 
tain  statistically  significant  data  in  a  pulse  height  spectrum.  In  order  to  predict  this 
detector  response,  one  must  account  for  material  and  geometry  considerations  such 
as  neutron  flux,  capture  cross  section,  detector  collection  volume,  and  conversion 
reaction  branching  ratios. 

7.4  Material  and  Geometry  Considerations 

Assuming  that  a  charge  pulse  resulting  from  a  15  'Gd  neutron  capture  event  could 
be  detected,  the  idealized  pulse  rate  fs  may  be  determined  as 

fa  =  tfvoi  Vd  K  [sec'1] .  (95) 

The  first  term,  volumetric  reaction  rate  Rvo i,  depends  directly  on  the  macroscopic 
capture  cross  section  Ec  and  the  thermal  neutron  flux  or 

Avoi  =  £c  iyth  ■  (96) 

The  macroscopic  capture  cross  section  is  calculated  as 

Y*c  =  N  ac  ,  (97) 

where  <jc  is  the  microscopic  neutron  capture  cross  section  in  barns  [10-24  cm2]  and 
N  is  the  concentration  of  nuclei  [cm'3],  calculated  as 

N  =  pNav  , 

Maw 


137 


via  the  material  density  p,  Avogadro’s  number  Nav,  and  the  molar  mass  Maw- 
The  second  and  third  terms  in  (95)  are  the  detector  collection  volume  Vd)  and  the 
conversion  reaction  branching  ratio  /C.  Collection  volume  is  calculated  as 

Vd  =  AcW,  (99) 

where  Ac  is  the  contact  area  and  W  is  the  depiction  region  of  the  solid  state  device. 
Assembling  the  components  from  (95) — (99) ,  the  idealized  pulse  rate  is  expressed  as 

fs  =  (°c  *th  ACW  K)  .  (100) 

Assuming  that  collection  volumes  are  identical  for  hypothetical  B5C  and  Gd- 
based  diodes,  and  that  both  are  subject  to  an  identical  neutron  flux,  fs  would  scale 
linearly  with  the  macroscopic  capture  cross  section  and  branching  ratio,  making  a 
comparison  between  the  two  materials  straightforward.  Using  parameters  summa¬ 
rized  in  Table  6,  and  assuming  100%  10B  and  15 'Gd  in  the  diodes,  the  ratio  of  charge 
pulse  rates  in  the  10B  neutron  capture  reaction  (88)  and  the  157Gd  neutron  capture 
reaction  in  Figure  38  is  /S(B5C)  /  /s(GdN)  =  3. 2/2. 3  =  1.4.  Therefore,  assuming 

Table  6.  Diode  parameters  for  charge  pulse  rate  comparison.  Parameters  common  to 
both  materials  are:  $th  —  104  n  cm“2s_1;  Ac  —  7r(0.1  cm)2/4;  W  —  4.2  pm 


Material 

p  (g  cm  3) 

M AW  (g  mol  x) 

<jc  (barn) 

1C 

b5c 

2.52 

55.2 

3840 

0.94 

GdN 

2.30 

171.2 

2.55  xlO5 

(0.113)  (0.30) 

the  pulse  processing  electronics  can  be  optimized  to  account  for  the  low  charge  pro¬ 
duction  from  a  15 'Gd  neutron  capture  event,  the  ratio  of  charge  pulse  rates  above 
predicts  that  a  pulse  height  spectrum  for  a  Gd-based  diode,  such  as  that  in  Figure 
41  for  B5C,  would  not  be  significantly  reduced  on  the  basis  of  the  poor  branching 
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ratio  for  internal  conversion  electron  in  15 'Gd,  nor  significantly  improved  owing  to 
its  favorable  neutron  capture  cross  section. 

With  questions  (1)  and  (2)  from  Section  7.1  addressed,  question  (3)  requires  con¬ 
sideration.  During  device  operation,  the  depletion  width,  and  hence  charge  collection 
volume,  are  controlled  by  the  bias  conditions.  Making  the  assumption  that  a  Gd- 
based  material  could  be  grown  on  a  substrate  with  a  similar  bandgap,  like  p-Si,  such 
that  an  abrupt  homojunction  diode  could  be  constructed,  W  is  determined  from 


W 


2es 

q 


(Na±Nd\ 

V  nand  ) 


1/2 


(Hi  -  V) 


(101) 


The  dielectric  constant  es  and  acceptor  NA  and  donor  Np  concentrations  are  ma¬ 
terial  properties  and  can,  to  some  extent,  be  influenced  during  device  fabrication. 
Additionally,  the  built  in  voltage  Vu  is  also  a  function  of  doping  concentration  as 


Vu 


(NAND\ 

V  nl  ) 


(102) 


where  nt  is  the  intrinsic  carrier  density.  For  the  sake  of  brevity,  full  development 
of  (101)  and  (102)  will  not  be  shown,  but  can  be  found  in  various  semiconductor 
textbooks  [11,12], 

The  requirement  to  optimize  the  depletion  width  in  a  semiconductor  device  be¬ 
comes  apparent  from  inspection  of  pulse  rate  (95)  and  transient  current  (92).  Specif¬ 
ically,  any  improvement  to  detection  efficiency  (fs)  has  the  opposite  effect  upon  the 
minimum  detectable  level  (A)  within  a  detector.  Furthermore,  the  effects  of  growth 
(doping)  and  device  operation  (bias)  allow  depletion  widths  to  be  engineered  via 
(101)  such  that  they  approach  the  mean  free  path  of  the  neutron  capture  reaction 
products.  Figure  44  shows  the  bias  dependent  depletion  width  in  a  hypothetical 
Gd-based  homojunction  neutron  detector,  using  estimates  of  Np  =  1019,  NA  =  1014, 
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nt  =  1010,  and  e  =  9.  As  an  example,  during  a  low  power  application  with  a  bias 


Figure  44.  Depletion  width  as  a  function  of  bias  for  a  hypothetical  Gd-based  (100% 
157Gd)  homojunction  neutron  detector  using  iVo  =  1019  and  Na  =  1014. 

of  —  1  V,  W  =  4.2  /ini.  Detailed  parameter  studies  would  reveal  the  full  influence 
of  depletion  width  upon  the  design  considerations  discussed  above,  but  any  analysis 
should,  as  a  minimum,  address  the  mean  free  paths  of  the  electrons  and  photons 
summarized  in  Table  5. 

Recall  that  a  narrow  beam  of  monoenergetic  photons,  with  an  incident  intensity 
Jo,  interacting  with  a  layer  of  material  with  thickness  x  and  density  p,  will  emerge 
with  intensity  /  given  by  the  exponential  attenuation  law 

I(x)  =  I0e~^px  ,  (103) 

where  p/p  is  the  photon  mass  attenuation  coefficient,  which  is  a  function  of  incident 
photon  energy  and  absorber  material  density.  Alternatively,  photons  can  be  char¬ 
acterized  by  their  mean  free  path  A,  which  is  the  average  distance  traveled  in  an 
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absorber  material  before  an  interaction  takes  place,  or 


/0°°  xe  'lxdx  1 

J0°°  e~fJjXdx  p 


(104) 


NIST  values  [13]  of  p/p  for  gadolinium  oxysulfide  at  photon  energies  of  182  keV, 
79  keV,  43  keV,  and  7  keV  from  Table  5  are  0.481,  2.613,  5.854,  and  130.6  cm2  g-1, 
respectively,  and  p  =  2.3  g-cm-3.  The  resulting  mean  free  paths  are  9  mm,  2  mm, 
700  /mi,  and  33  pm  for  the  respective  photon  energies  above.  The  possibility  of  a 
7  keV  photon  interaction  with  the  material  cannot  be  ruled  out  from  these  calcu¬ 
lations.  However,  the  charge  generated  from  a  single  event  interaction  would  be 
~  0.1  fC,  well  below  the  intrinsic  noise  level  of  a  commercial  preamplifier.  There¬ 
fore,  the  effects  of  photon  interaction  within  the  depletion  widths  of  Figure  44  of  the 
hypothetical  diode  are  negligible. 

Based  on  limited  availability  of  data,  electron  interactions  will  be  addressed 
more  qualitatively.  Theoretical  electron  mean  free  path  calculations  in  solids  [14, 15] 
are  limited  to  electron  kinetic  energies  up  to  ~  10  keV,  and  range-energy  plots  for 
more  energetic  electrons  are  available  for  silicon  [5],  but  are  not  obtainable  for  non- 
conventional  semiconductor  devices.  However,  in  silicon,  the  internal  conversion 
electrons  created  from  neutron  capture  by  157Gd  have  ranges  of  approximately  40- 
160  pm  [5].  Electron  stopping  power,  and  hence  range,  results  from  collisional  and  ra¬ 
diative  electron  energy  losses,  and  scales  linearly  with  electron  density  n.  Therefore, 
the  ratio  of  nsi/ncd  ~  1/5  provides  a  very  rough  scaling  factor  for  approximating  the 
internal  conversion  electron  range  in  a  Gd-based  semiconductor  device.  Based  upon 
this  scaling  factor,  the  range  of  internal  conversion  electrons  can  be  approximated  as 
5-20  pm. 

The  most  effective  method  to  determine  electron  ranges  is  by  Monte  Carlo  sim¬ 
ulations,  which  are  computational  methods  that  rely  on  random  sampling.  In  fact, 


141 


the  estimate  of  electron  range  using  fundamental  physics  and  reasonable  approxima¬ 
tions  is  in  fair  agreement  with  CASINO  and  GEANT4  simulations  [16,17]-  CASINO, 
which  predicts  electron  trajectories  in  solid  materials,  was  used  to  determine  the  in¬ 
ternal  conversion  electron  penetration  depth  in  a  Gd203  on  SiC  diode,  and  found  a 
mean  penetration,  for  full  energy  deposition,  of  ~  26  fim.  [16].  However,  out  of  the 
105  incident  electrons  with  79  keV  energy,  only  21  fully  deposited  their  energy  at 
all  ranges  between  8-42  /jrn  [16] .  While  many  of  the  remaining  electrons  would  have 
deposited  a  portion  of  their  energy  in  the  diode  material,  the  generated  charge  would 
almost  certainly  fall  below  the  intrinsic  noise  level  of  a  commercial  preamplifier. 

GEANT  is  an  abundant  set  of  physics  models  designed  to  handle  the  interactions 
of  particles  with  matter  across  a  very  wide  energy  range.  Like  CASINO,  GEANT 
uses  Monte  Carlo  simulations  to  replicate  charged  particle  interactions  with  matter, 
but  adds  the  flexibility  to  simulate  neutral  particles.  Figure  45  is  a  schematic  of 
a  GEANT  simulation  [17]  in  which  106  thermal  neutrons  were  incident  on  a  Gd- 
doped  Hf02  thin  film  grown  on  a  thick  n-silicon  wafer.  Natural  gadolinium  was  used 
(~  16%157Gd),  and  the  doping  level  of  the  Gd  in  Hf02  was  simulated  at  15%.  The 
depletion  width  at  1  V  reverse  bias  was  the  full  width  of  the  Gd  layer  and  3  /m i 
in  silicon.  Figure  46(a)  shows  that  roughly  0.2%  (or  about  2000)  of  the  incident 
neutrons  produced  electrons.  The  low  production  was  attributed  to  the  doping  level 
in  the  Hf02  [17],  and  the  large  range  of  electrons,  relative  to  the  depletion  width, 
is  another  likely  cause.  The  Auger  and  internal  conversion  electron  peaks  predicted 
by  theory  were  observed  in  the  simulation,  as  shown  in  Figure  46(a).  Figure  46(b) 
shows  the  energy  deposited  in  the  thin  film  and  the  number  of  electrons  that  fully 
deposited  their  energy.  As  an  example,  there  were  just  over  100  events  that  occurred 
where  30  keV  was  deposited  in  the  active  region  of  the  diode  (solid  black  line),  and 
of  those  events,  most  were  full  depositions  (red  dashed  line).  Thus,  a  few  events 
occurred  where  more  energetic  electrons  deposited  a  portion  of  their  energy  in  the 
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Figure  45.  Schematic  of  GEANT4  simulation  of  thermal  neutrons  incident  on  Hf02=Gd 
on  silicon  heterojunction.  The  doping  concentration  was  15%  natural  Gd. 


(a)  The  energy  with  which  each  electron  is  ini¬ 
tially  formed  following  neutron  capture  in  nat¬ 
ural  gadolinium. 


(b)  Solid  black:  Total  energy  deposited  by  elec¬ 
trons  in  the  active  volume  of  the  diode.  Dashed 
red:  Full  energy  deposition  by  electrons  in  the 
active  detector. 


Figure  46.  GEANT4  simulation  output  [17]. 
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active  region.  The  dominant  deposition  at  30  keV  in  Figure  46(b)  corresponds  to  a 
generated  charge  of  0.4  fC,  which,  unfortunately,  is  below  the  intrinsic  noise  level  of 
many  commercial  preamplifiers. 

The  CASINO  and  GEANT4  simulations  [16, 17]  reinforce  the  challenges  of  low 
charge  production  from  neutron  capture  by  157Gd  and  large  electron  range  relative 
to  typical  diode  dimensions.  However,  there  exists  an  opportunity  to  optimize  the 
depletion  width  via  doping  concentrations  during  him  growth  and  bias  conditions 
during  operation,  which  will  be  discussed  later. 

The  disparity  in  capture  cross  section  between  the  two  materials  also  affects  the 
quantity  of  neutron  interactions.  Recall  that  in  a  simplistic  one-dimensional  analysis, 
if  a  monoenergetic  narrow  beam  neutron  experiment  were  conducted,  the  entering 
neutron  flux  <h0  would  attenuate  exponentially  [5,18],  via  the  neutron  capture  pro¬ 
cess,  with  absorber  thickness  x  as 

<h(a;)  =  <h0  e~T‘cX  .  (105) 

Alternatively,  (105)  permits  the  fraction  of  neutrons  deposited  in  material  T  to  be 
expressed  as 

I  =  1  -  =  1  -  e~^cX  .  (106) 

A  graph  of  (106)  is  shown  in  Figure  47  for  B5C  (100  %  10B)  and  GdN  (100  %  157Gd) 
homojunctions  and  clearly  indicates  that  the  massive  cross  section  of  157Gd  most 
significantly  affects  device  design  for  the  parameter  of  him  thickness.  Simply  stated, 
an  equivalent  number  of  neutron  capture  events  would  occur  for  much  thinner  Gd- 
based  films  compared  to  B,5C  films.  As  an  example,  the  fraction  of  neutrons  absorbed 
in  a  10  /irn  B5C  him  would  required  a  157Gd  him  thickness  of  0.5  /mi. 

In  theory,  him  thickness  appears  to  be  a  signiheant  design  advantage  for  Gd- 
based  devices  over  their  boron  counterparts.  However,  in  practice,  and  as  discussed 
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Figure  47.  Fraction  of  absorbed  neutrons  as  a  function  of  film  thickness  for  B5C  and 
GdN.  A  one  dimensional  monoenergetic  narrow  beam  of  thermal  neutrons  is  assumed. 


earlier,  the  range  of  internal  conversion  electrons  created  from  neutron  capture  by 
15 'Gd  is  likely  about  20-25  prn,  and  is  larger  than  the  range  of  the  heavy  charged 
particles  created  from  neutron  capture  by  B5C.  Therefore,  thicker  depletion  widths 
would  be  required  to  achieve  full  energy  deposition  of  the  neutron-capture-generated 
internal  conversion  electrons,  which  might  also  limit  the  gamma-blindness  of  a  de¬ 
tector.  Furthermore,  to  generate  a  depiction  width  of  25  /irn  would  require  a  bias  of 
—62.1V,  using  the  parameters  discussed  for  Figure  44.  Thus,  low  power  operation 
would  be  inhibited. 

Figure  48  shows  the  dependence  of  mean  free  path  on  photon  energy  for  natural 

Gd,  determined  from  photon  energy  dependent  p/p  values  from  the  NIST  database 

[13]  and  using  p  =  7.9  g  cm-3  for  Gd;  GaN  has  p  =  6.1  g  cm-3.  The  figure  indicates 
that  a  diode  depleted  to  the  full  25  pm  range  of  internal  conversion  electrons  created 

from  neutron  capture  by  157Gd  will  be  gamma  blind  to  photons  with  hv  >17  keV. 

The  possibility  of  depletion  region  interactions  between  photons  with  energies  below 

this  threshold  and  a  diode  with  a  comparable  material  density  cannot  be  ruled  out, 
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Figure  48.  Photon  energy  dependent  mean  free  path  A,  determined  from  NIST- 
tabulated  [13]  values  of  p/p  for  natural  Gd. 

based  upon  the  anticipated  reaction  dependent  photon  energies  discussed  in  Section 
7.2.2.  The  likelihood  also  exists  for  higher  energy  photons  to  be  Compton  scattered 
below  the  gamma  discrimination  energy.  However,  a  17  keV  energy  deposition  in 
the  active  region  of  an  ideal  Gd^Gai.^N  p  —  n  homojunction  diode  corresponds  to  a 
generated  charge  of  0.4(x  =  0)-0.9(x  =  1)  fC,  which  is  below  the  intrinsic  noise  level 
of  many  commercial  preamplifiers.  Thus,  it  is  expected  that  a  diode  depleted  to  the 
25  pm  range  of  79  keV  internal  conversion  electrons  will  be  either  gamma  blind  or 
charge  insensitive  to  gamma  ray  interactions  within  the  depletion  region. 

In  the  end,  the  calculations  presented  here  are  made  possible  under  the  critical 
assumption  that  the  charge  generated  in  a  Gd-based  semiconductor  by  an  internal 
conversion  electron  created  from  157Gd  neutron  capture  exceeds  the  intrinsic  noise 
level  of  the  preamplifier  required  for  pulse  height  analysis.  Based  upon  previous 
research  [1,16]  and  topical  review  of  solid  state  neutron  detection  [19],  this  criti¬ 
cal  assumption  may  be  invalid  and  gadolinium  may,  in  fact,  be  “romanticized,”  as 
predicted  by  many  researchers.  This  is  not  to  say  that  development  of  novel  solid 
state  neutron  detectors  using  gadolinium  should  be  abandoned,  but  until  the  pulse 
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processing  electronics  necessary  to  detect  a  neutron  induced  signal  are  improved, 
the  endeavor  may  prove  vexing.  On  the  other  hand,  I  find  the  words  of  John  D. 
Rockefeller  compelling: 


I  do  not  think  there  is  any  other  quality  so  essential  to  success  of  any 
kind  as  the  quality  of  perseverance.  It  overcomes  almost  everything,  even 
nature. 
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VIII.  Conclusions 


8.1  Summary  of  Findings 

During  the  past  decade,  rare  earth  doped  semiconductors  have  generated  con¬ 
siderable  attention  for  their  application  in  new  optoelectronic  devices  [1-4],  but  in¬ 
complete  understanding  of  the  physical  and  electronic  properties  of  these  materials 
classifies  rare  earth  doped  GaN  heterostructures  as  non-convent ional,  as  compared 
to  silicon-based  devices.  To  determine  the  extent  to  which  even  low  concentrations 
of  a  rare  earth  in  a  GaN  host  can  alter  device  properties,  investigations  of  the  sur¬ 
face  electronic  structure  and  interface  properties  of  RExGa1_.rN  (RE  =  Gd,  Er,  Yb) 
semiconductors  are  necessary. 

With  these  considerations  in  mind,  X-ray  photoemission  spectroscopy  was  used  to 
determine  the  effective  Debye  temperatures  of  gallium  and  ytterbium  in  GaN:Yb  thin 
films.  The  similarity  between  their  measured  values  suggested  that  substitutional 
occupation  of  a  Ga  site  by  a  Yb  ion  occurs,  which  supports  experimental  data  and 
theoretical  calculations.  Careful  comparison  of  the  Ga  and  Yb  values  indicated  a 
slightly  smaller  effective  Debye  temperature  of  the  Yb  atom,  which  reflects  a  surface 
that  is  softened  and  strained,  possibly  due  to  an  increased  Yb-N  bond  length  as 
compared  to  the  Ga-N  bond  length. 

Next,  resonant  photoemission  for  GaN  thin  films  doped  with  various  rare  earths 
(RE  =  Gd,  Er,  Yb)  was  investigated  using  synchrotron  based  photon  energy  de¬ 
pendent  photoemission  spectroscopy,  in  order  to  probe  what  valence  bands  of  the 
semiconductors  have  strong  4/  and/or  rare  earth  weight.  The  calculated  Fano  pro¬ 
files  predicted  the  experimentally  observed  4 d  4/  super  Coster-Kronig  processes. 
Additionally,  strong  rare  earth  4/  weights  were  observed  about  5-6  eV  below  the 
valence  band  maximum  for  Gd  and  Er  doped  GaN.  For  Yb,  the  photoemission  res¬ 
onance  was  much  weaker,  indicating  that  there  was  only  a  very  little  depletion  of 
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the  4/14  occupancy,  but  strong  hybridization  with  GaN  was  implicated.  There  was 
also  evidence,  from  the  selectivity  of  the  resonant  photoemission  enhancement  of 
the  valence  band,  for  strong  4/  hybridization  with  the  GaN  valence  bands,  particu¬ 
larly  within  the  Er  4/5/2  envelope.  Furthermore,  the  resonant  photoemission  results 
strongly  supported  the  predicted  hybridization  with  the  host  lattice  [5],  and  the  ex¬ 
pected  f  —  s  hybridization  [5-8].  The  results  on  the  placement  of  the  occupied  Gd,  El¬ 
and  Yb  4/  levels,  deep  within  the  valence  band,  suggested  that  the  intra-atomic  /-/ 
transitions  may  be  more  ‘blue’  than  predicted  by  many  theoretical  models,  which 
has  implications  for  both  the  physics  and  the  device  fabrication  and  engineering 
communities. 

Device  fabrication  and  neutron  detector  applications  were  explored  more  carefully 
by  measuring  the  Schottky  barrier  heights  at  the  gold  to  rare  earth  doped  GaN 
interface  for  the  Gd,  Er,  and  Yb  dopants  via  ultraviolet  photoemission  spectroscopy. 
The  An  overlayer  did  not  wet  and  cover  the  GaN  surface,  even  with  rare  earth 
doping  of  the  GaN.  The  photoemission  studies  showed  that  the  Schottky  barrier 
heights  between  Au  and  RE:GaN  thin  him  samples  were  measured  to  be  1.68  ±  0.1 
eV  (Yb:GaN),  1.64  ±  0.1  eV  (Er:GaN),  and  1.33  ±  0.1  eV  (Gd:GaN),  which  was 
25-55%  larger  than  those  reported  at  the  gold  to  undoped  GaN  interface  [9-25]. 
This  trend  of  the  Schottky  barrier  heights  followed  the  trend  of  the  rare  earth  metal 
work  function.  These  results  were  uniform  for  all  of  the  rare  earths  studied,  and 
imply  further  reduction  in  leakage  currents  as  compared  with  undoped  GaN  devices. 
Moreover,  the  larger  barrier  height  at  the  metal  to  rare  earth  doped  GaN  interface 
means  that  radiation-induced  charge  collection  volumes  are  increased,  and  thus  more 
efficient,  over  what  would  be  the  case  in  pure  GaN  devices. 

Last,  the  utility  of  gadolinium  as  a  neutron  detection  material  in  a  hypothetical 
solid  form  direct  conversion  homojunction  was  examined  using  fundamental  nuclear 
and  semiconductor  physics,  and  compared  to  a  boron  carbide  benchmark.  In  the 
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single  neutron  detection  limit,  the  charge  produced  in  a  Gd-based  semiconductor 
by  an  internal  conversion  electron  created  from  15 'Gd  neutron  capture  was  found 
to  approach  the  intrinsic  noise  level  of  commercial  preamplifiers  required  for  pulse 
height  analysis.  In  general,  the  fraction  of  captured  neutrons,  and  hence  efficiency, 
was  maximized  with  Gd- layer  thickness,  and  was  found  to  exceed  that  of  10B  layers 
by  at  least  an  order  of  magnitude.  However,  the  range  of  internal  conversion  elec¬ 
trons  created  from  neutron  capture  by  15 'Gd  was  estimated  to  be  about  20-25  /irn, 
which  is  larger  than  the  range  of  the  heavy  charged  particles  created  from  neu¬ 
tron  capture  by  B,5C.  Therefore,  thick  Gd  films  (~  30-40  //in)  would  be  required  to 
achieve  full  energy  deposition  of  the  neutron-capture-generated  internal  conversion 
electrons,  which  might  also  limit  the  gamma-blindness  of  a  detector.  Moreover,  full 
depletion  of  a  Gd-based  homojunction  to  the  range  of  the  72  keV  internal  conversion 
electrons  required  large  bias  of  V  =  —62.1V,  thus  inhibiting  low  power  operation. 
For  gadolinium  to  warrant  serious  consideration  as  a  neutron  detection  material  will 
require  maturation  of  preamplifier  technology,  such  that  the  gain  can  be  minimized 
to  account  for  inherently  low  charge  production. 

8.2  Future  Work 

There  are  several  natural  experimental  and  computational  extensions  of  the  cur¬ 
rent  hirelings  toward  future  research.  Researchers  who  are  interested  in  semiconductor 
device  characteristics  could  fabricate  multiple  Schottky  diodes  on  a  single  GaN:RE 
thin  him.  The  Schottky  barrier  could  then  be  calculated  from  device  current- voltage 
characteristics,  as  discussed  in  Chapter  II.  Ideally,  metal  contacts  would  be  sputtered 
onto  thin  films  to  produce  multiple  diodes  from  a  single  sample,  and  I-V  character¬ 
istics  could  be  obtained  either  from  a  conventional  source  measurement  unit  (via 
soldered  wires)  or  directly  through  a  probe  station  coupled  with  a  semiconductor 
characterization  system.  To  summarize  the  method,  the  current-voltage  characteris- 
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tic  of  a  metal-semiconductor  contact  is 


Js  {eqV/kBT  -  1)  , 

(107) 

^*2n2e-4,s,„/fcsT 

(108) 

where  Js  is  the  saturation  current  and  A*  is  the  effective  Richardson  constant,  in 
units  of  A/K2-cm2,  which  depends  directly  on  the  electron  effective  mass  in  the 
semiconductor  material. 

Experimental  calculation  of  Schottky  barrier  heights  via  equations  (107)  and 
(108)  is  done  as  follows  and  as  shown  in  Figure  49.  For  values  of  V  such  that 


V 

Figure  49.  Logarithm  of  forward  current  density  versus  applied  voltage  for  NTE555 
PIN  Schottky  Diode  (points)  and  least  squares  fit  (line)  in  the  linear  region  of  the 
diode  characteristic. 


qV  >  3 ksT,  it  follows  that  J  ~  JseqV^BT .  The  logarithm  of  current  density 


In  J  =  In  Js  + 


qV 
~^bT  ’ 


(109) 


versus  applied  forward  voltage  is  plotted  and  extrapolated  to  V  —  0,  yielding  a  value 
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for  Js,  from  which  the  barrier  height  is  obtained  via  (108)  and  knowledge  of  the  effec¬ 
tive  Richardson  constant.  Figure  49  shows  this  procedure  for  a  commercial  NTE555 
Schottky  barrier  PIN  diode,  which  yielded  a  calculated  saturation  current  density  of 
1.21  nA/cm2.  Assuming  a  circular  diode  contact  area  of  0.55  mm  [26],  and  A*  =  110 
for  n-Si  [27],  (108)  yields  a  Schottky  barrier  height  of  0.94  eV.  Experimental  results 
of  various  metal-n-Si  systems  are  known  [28]  to  obey  the  empirical  relationship 

$s,n  =  0.27</>m  -  0.52  .  (110) 

Therefore,  the  calculated  Schottky  barrier  height  of  0.94  eV  indicates  that  a  high 
work  function  contact  metal  (f)m  =  5.43  eV  was  used  in  the  NTE555  diode.  Although 
the  manufacturer’s  specifications  do  not  identify  the  contact  metal,  common  high 
work  function  metals  such  as  Au  and  Pt  have  work  function  values  that  range  from 
5.10-5.42  eV  and  5.12-5.93  eV,  respectively. 

When  the  experimental  data  are  plotted  against  the  theoretical  current  density 
(107)  using  &B,n  =  0.94  eV  and  the  other  parameters  specified  above,  there  is  poor 
agreement  between  the  results,  as  shown  in  Figure  50.  The  relatively  weak  correlation 
between  data  and  theory  is  attributable  to  the  fact  that  (107)  assumes  an  ideality 
factor  rj  of  1.0.  In  practice  [12,13,29,30],  a  more  accurate  expression  for  the  current 
density  at  metal-semiconductor  interfaces  and  p  —  n  junctions  [27,28]  is 

J  =  JS  (eqV/rikBT  -  1)  ,  (111) 


where  1  <  rj  <  2. 

For  typical  p—n  junctions,  when  ideal  diffusion  current  dominates,  r/  =  1,  whereas 
when  recombination  current  dominates,  rj  =  2.  Accordingly,  when  both  currents  are 
comparable,  rj  has  a  value  between  1-2.  For  typical  Schottky  diodes  with  doping 
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V 

Figure  50.  Logarithm  of  forward  current  density  versus  applied  voltage  for  NTE555 
PIN  Schottky  Diode  (points)  and  theoretical  current  density  (lines)  for  a  Schottky 
contact  with  &B,n  —  0.94  eV  and  ideality  factors  rj  —  1  (yellow)  and  rj  —  1.55  (purple). 

concentrations  <  1017  cm-3  operated  at  room  temperature,  the  dominant  current 
transport  mechanism  is  thermionic  emission  of  majority  carriers  from  the  semicon¬ 
ductor,  over  the  potential  barrier,  and  into  the  metal  [27].  In  this  case,  where 
tunneling  current  and  depletion  layer  recombination  are  minimized,  J  is  determined 
almost  entirely  by  thermionic  emission  and  r)  is  close  to  unity.  However,  as  doping 
concentrations  increase,  the  probabilistic  process  of  tunneling  increases,  owing  to  the 
increased  occupancy  of  available  states  in  the  conduction  band  of  the  semiconductor. 
In  this  case,  tunneling  currents  increase  and  rj  can  be  much  larger  than  unity.  The 
relationship  between  tunneling  current  and  doping  concentration  for  Au-Si  Schottky 
diodes  [28]  shows  that  rj  «  1.55  when  Njj  ps  1019  cm-3.  Although  the  complete  design 
details  of  the  NTE555  Schottky  barrier  PIN  diode  are  proprietary,  the  diode  cannot 
be  classified  solely  as  either  a  p  —  n  junction  or  a  Schottky  diode.  Therefore,  it’s 
likely  that  the  calculated  value  of  r)  —  1.55  for  the  NTE555  is  due  to  a  combination 
of  p  —  n  junction  recombination  and  Schottky  barrier  tunneling  currents. 

Regardless  of  the  success  or  failure  to  reverse  engineer  design  characteristics  of 
the  NTE55  Schottky  barrier  PIN  diode,  the  process  to  calculate  the  Schottky  barrier 
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via  current-voltage  characteristics  is  valid  and  may  be  applied  to  thin  films  used 
during  this  research.  The  calculated  barrier  heights  from  I-V  characteristics  could 
be  compared  to  the  values  obtained  via  photoemission  spectroscopy.  If  successful 
calculation  of  Schottky  barrier  heights  via  I-V  characteristics  were  accomplished  for 
all  of  the  rare  earth  dopants  studied  during  this  research,  the  results  may  merit 
publication. 

Researchers  who  are  motivated  by  device  modeling  could  examine  device  param¬ 
eters  and  conduct  sensitivity  analyses  for  a  Gd-based  diode  or  other  device  structure 
in  a  thermal  neutron  environment.  This  research  endeavor  would  most  likely  re¬ 
quire  integration  of  multiple  computational  tools,  to  include  the  GEometry  ANd 
Tracking  (GEANT)  toolkit,  which  can  be  used  to  accurately  simulate  the  passage  of 
both  charged  and  neutral  particles  through  matter,  as  well  as  a  finite  element  device 
simulation  code,  such  as  Synopsis  TCAD,  to  model  the  semiconductor  device  struc¬ 
ture  and  characteristics.  At  least  one  researcher  [31]  has  already  conducted  similar 
research  for  B5C  heterostructures,  which  would  provide  an  interesting  comparison. 

One  type  of  heterojunction  application  that  may  warrant  consideration  is  the 
quantum  well  [28],  depicted  in  Figure  51.  Quantum  wells  use  the  differences  between 
conduction  band  A Eq  and  valence  band  A Eq  levels  between  two  hetero junctions  or 
three  layers  of  materials  in  order  to  form  barriers  for  charge  carriers.  In  Figure  51, 
the  middle  GdN  layer,  of  relatively  low  bandgap  ( Eg  ~  1  eV),  is  inserted  between 
GaN:Gd  layers,  of  relatively  large  bandgap  ( Eg  ~  3  eV),  such  that  the  middle  layer 
has  the  lowest  Eq  for  an  electron  well  and  the  highest  Ey  for  a  hole  well.  Whereas 
electrons  in  a  bulk  semiconductor  are  free  to  move  in  three  dimensions,  a  quantum 
well  confines  electrons  or  holes  in  a  two-dimensional  system.  In  this  2-D  system,  the 
solutions  for  the  electron  and  hole  wavefunctions  inside  the  quantum  well  obey  the 
Schrodinger  equation,  and  the  energy  eigenvalues  are  quantized  at  distinct  levels, 
shown  as  E\  and  E2  in  the  figure.  A  critical  parameter  for  quantum  well  design 
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Figure  51.  Single  quantum  well  schematic. 


is  the  well  width  Lx.  If  the  well  width  is  smaller  than  the  de  Broglie  wavelength 
of  the  electron  or  hole,  electron  and  hole  wavefunctions  between  adjacent  quantum 
well  structures  begin  to  overlap,  and  distinct  energy  minibands  are  formed  in  a 
superlattice  structure,  as  shown  in  Figure  52. 
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Figure  52.  Superlattice  structure  schematic. 


If  multiple  quantum  wells  can  be  grown  and  fabricated  into  a  heterostructure 
superlattice,  such  that  the  generated  charge  from  157Gd  neutron  conversion  can  be 
collected  efficiently,  the  challenge  of  low  charge  generation  may  be  abated.  At  this 
point,  a  quantum  well  or  superlattice  approach  to  neutron  detection  using  Gd  is  a 
concept  only,  as  good  heteroepitaxy  growth  with  the  different  materials  is  relatively 
unproven.  However,  under  the  assumption  that  structures  could  be  fabricated,  de- 
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tailed  device  modeling  may  provide  a  reasonable  estimate  of  the  charge  collection 
efficiency  of  these  devices. 

An  extension  of  the  combined  achievements  of  both  research  efforts  discussed  in 
the  previous  paragraphs  could  be  measurement  of  the  radiation  response  of  a  Gd- 
based  diode  in  a  thermal  neutron  environment,  such  as  the  AF1T  thermal  pile  or 
The  Ohio  State  University  research  reactor,  and  comparison  with  the  computational 
model. 

Lastly,  researchers  who  are  interested  in  the  surface  electronic  structure  and  in¬ 
terface  properties  of  thin  him  semiconductors  might  use  photoemission  spectroscopy 
to  duplicate  many  of  the  experiments  conducted  during  this  research  effort,  but  with 
an  actinide  rare  earth,  such  as  uranium.  Plans  to  pursue  this  endeavor  already  exist. 
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